Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,264 A→C 8.7% intergenic (+34/‑183) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2640AC 8.7% 247.9 / 8.6 107intergenic (+34/‑183)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (55/41);  new base C (3/7);  total (58/49)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 108 total aligned reads displayed.

GTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTA  >  NC_000913/4060116‑4060409
                                                                                                                                                    |                                                                                                                                                 
gTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa                                                                                                                                                   >  2:554382/1‑149 (MQ=255)
    aTCGATGATGACGAACTGGTAGAAGTGACTCCTCCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAAAGATCGTCGCCGCGCCAACCGCGCCCCGAAAGACGATTAATTACGTTTTTTATTCATAAAAAACCAGCCAGCGATGGc                                                                                                                                               >  2:745812/1‑149 (MQ=255)
    aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACTGAAAACGATCGTCGCCGCTCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                               <  2:372953/149‑1 (MQ=255)
    aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                               >  1:679638/1‑149 (MQ=255)
    aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                               >  1:182091/1‑149 (MQ=255)
     tctatgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTAATAAAAAACCTGACAGCGATGGCa                                                                                                                                              >  1:4091/4‑149 (MQ=255)
      cGATGAAGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAg                                                                                                                                             <  2:256982/149‑1 (MQ=255)
       gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAgg                                                                                                                                            <  2:507823/149‑1 (MQ=255)
       gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                                                                            <  1:238665/149‑1 (MQ=255)
       gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATGTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                                                                            >  2:928876/1‑149 (MQ=255)
       gatgatGACGAACTGGTAGAAGTGACTCAGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAgg                                                                                                                                            <  1:454886/149‑1 (MQ=255)
         tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                                                                          >  1:980042/1‑149 (MQ=255)
         tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATAAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                                                                          >  2:789847/1‑149 (MQ=255)
         tgatgaAGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                                                                          >  2:729010/1‑149 (MQ=255)
           atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCAGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                                                                        <  2:8540/149‑1 (MQ=255)
            tgaCGAACTGGTAGAAGTGAATCCGACCTCTATCCGTATTCGTAAACGTCACCTGAAGGAAAACGATCGTCGGCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttt                                                                                                                                       <  2:881873/149‑1 (MQ=255)
             gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCACCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATTAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                                                                      >  1:666749/1‑149 (MQ=255)
             gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGGCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                                                                      <  2:72038/149‑1 (MQ=255)
             gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGGCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                                                                      <  2:342543/149‑1 (MQ=255)
               cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttat                                                                                                                                    <  2:378487/149‑3 (MQ=255)
                  aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGc                                                                                                                                 <  1:773105/149‑1 (MQ=255)
                   cTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGAGATGGCAGGTTTTTTTTTGcc                                                                                                                                >  1:654407/1‑149 (MQ=255)
                      ggaGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTgg                                                                                                                             >  2:343749/3‑148 (MQ=255)
                      gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGCCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGt                                                                                                                             >  1:1006034/1‑149 (MQ=255)
                      gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTgg                                                                                                                             >  2:843869/1‑148 (MQ=255)
                      gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGt                                                                                                                             >  1:1006895/1‑149 (MQ=255)
                       tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTc                                                                                                                            <  1:395703/149‑1 (MQ=255)
                          aaGTGACTCCGACCTCTATCCGTATTCGTATACGTCACCTGCCGGAAAACGATCGTCGCCGCGCCAACCGCGCGCCGAAAGACGATTAATTTCGTTATTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                                                                         >  2:813259/1‑149 (MQ=255)
                            gTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTAGCCGCGCCAACCGCGCACCGAAAGACGATTAATGTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAAGGCACGTTATTATTTGCCTGTCGTGAg                                                                                                                       >  1:743379/1‑149 (MQ=255)
                                 tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAGCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcac                                                                                                                  >  2:917311/1‑149 (MQ=255)
                                    gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTTTCGTGAGCGCACAcc                                                                                                               >  1:203800/1‑149 (MQ=255)
                                      ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgcg                                                                                                                     >  2:517901/1‑141 (MQ=255)
                                      ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgcg                                                                                                                     <  1:517901/141‑1 (MQ=255)
                                        tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCGACCGCGCCCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                                                                           >  2:79163/1‑149 (MQ=255)
                                         ctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                                                                          >  2:154214/1‑149 (MQ=255)
                                           aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATTGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                        >  2:849480/1‑149 (MQ=255)
                                           aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCCCACCAATGCGa                                                                                                        >  1:110846/1‑149 (MQ=255)
                                           aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                        >  2:264538/1‑149 (MQ=255)
                                           aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                        >  2:241263/1‑149 (MQ=255)
                                              cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                                                                     >  2:65931/1‑149 (MQ=255)
                                                 aTTAGTAAACGTAACCTGACGGAAAACGATGGTCGCCGCGACAACCGCGCACCGAAAGACGTATAATTTCGTTGTTTAGACATAAAAAACCTGCCAGCGAGGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                                                                  <  1:663320/149‑1 (MQ=255)
                                                  tgagTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                                                                 <  1:459898/146‑1 (MQ=255)
                                                      tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGCTGGCAGGTTTTGTTTTGCCTGTCGTGAGCGCACACAAATGCGAGTATCTCCCtt                                                                                             >  1:1009663/1‑149 (MQ=255)
                                                           gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                                                                        >  2:301200/1‑149 (MQ=255)
                                                            tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                                                                       <  1:365051/149‑1 (MQ=255)
                                                            tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                                                                       <  1:360752/149‑1 (MQ=255)
                                                                cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                                                                   >  2:1026038/1‑149 (MQ=255)
                                                                 tGACGGAAAACGTTAGTCGCCGTGCCAACCGGGGACCGTAAGTCTATTAATTTCGTTTTTTAGTCATAAAATACCTGCCAGCGTTGGCATGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGAGAGTATCTCCCTTATAAGTCtgtg                                                                                  <  2:517659/149‑1 (MQ=255)
                                                                 tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGTCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                                  <  2:756276/149‑1 (MQ=255)
                                                                 tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                                  <  1:568278/149‑1 (MQ=255)
                                                                 tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                                  >  1:12150/1‑149 (MQ=255)
                                                                  gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                 <  2:803289/149‑1 (MQ=255)
                                                                    cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGtt                                                                               <  1:660838/149‑1 (MQ=255)
                                                                     ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTTTCGTGAGCGCACACCAATGCTAGTATCTCCCTTATAAGTCTGTGGTTa                                                                              >  2:856053/1‑149 (MQ=255)
                                                                     ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACCCCAAGGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                              >  2:730330/1‑149 (MQ=255)
                                                                     ggAAAACGATCGTCGACGCGCCAACCGCGCACCGAAAGACGTTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATTGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                              <  1:677056/149‑1 (MQ=255)
                                                                            taTCGTAGCAGCGCAAAACGAGCACCGAAAGAAGATTAATTTAGTTTTTTAGTAATAAATAACCTGCCAGCGATGGCAGGTTTTTTTTTGACTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                                                                       <  1:986581/148‑1 (MQ=255)
                                                                            gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTGTGTGGTTACGTCAAt                                                                       >  2:224921/1‑149 (MQ=255)
                                                                              tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                                                                     >  2:277734/1‑149 (MQ=255)
                                                                                gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCTTGAGCGCACACCAATGCGATTATCTCCCTTATAAGTATGTGGTTACGTCAATagtg                                                                   >  2:957474/1‑147 (MQ=255)
                                                                                gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                                                                   >  2:296829/1‑149 (MQ=255)
                                                                                 tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATGTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                                                                  >  1:375777/1‑149 (MQ=255)
                                                                                   gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                                                                >  2:533492/1‑149 (MQ=255)
                                                                                    ccgcGCCCACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt                                                               >  1:133360/1‑149 (MQ=255)
                                                                                        gcCAACCGCGCACCGAAAGACGATTAATTTAGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt                                                           <  1:277742/149‑1 (MQ=255)
                                                                                         ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTTAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc                                                          <  2:797931/149‑1 (MQ=255)
                                                                                          cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca                                                         >  1:16662/1‑149 (MQ=255)
                                                                                           aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac                                                        >  1:275860/1‑149 (MQ=255)
                                                                                             ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGCTGGCAGGTTTTTTTTTGCCTGTCGGGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                      <  2:204262/149‑1 (MQ=255)
                                                                                               gcgaACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt                                                    <  1:581319/149‑1 (MQ=255)
                                                                                                   aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt                                                <  1:366790/149‑1 (MQ=255)
                                                                                                   aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt                                                <  1:104041/149‑1 (MQ=255)
                                                                                                   aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt                                                >  1:101648/1‑149 (MQ=255)
                                                                                                       aaaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACGAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTAACTGCaa                                            >  2:425762/1‑149 (MQ=255)
                                                                                                         aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATATAGAGCTTTGTCACAAATATCTGCaaat                                          >  2:831614/1‑148 (MQ=255)
                                                                                                          gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGt                                         >  2:686553/1‑149 (MQ=255)
                                                                                                            cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCa                                       >  2:885094/1‑149 (MQ=255)
                                                                                                            cGATTAATTTCGTATTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGATCGCACACCAATGCGAGTATCTCCCTTATAACTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCa                                       <  2:920616/149‑1 (MQ=255)
                                                                                                              aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTGTTTGCCTGTCGTGCGCGCACACCAATGCGAGGATCTCCCTTATAAGTCTGTGGTTACGTGAATAGAGCGCGTTGTCACAATTATCTTCAAAGTCATa                                     >  2:462191/1‑149 (MQ=255)
                                                                                                                     tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTGTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTTTCTGCCTTTTAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTGACAATTATGTGGAAAGTCATACGCCGtt                              >  2:492989/1‑149 (MQ=255)
                                                                                                                         tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAAtt                          <  1:65931/149‑1 (MQ=255)
                                                                                                                            ttattaaTACAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCATATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGct                       <  2:207851/143‑1 (MQ=255)
                                                                                                                              aGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGcttt                     <  1:798824/149‑1 (MQ=255)
                                                                                                                               gcaaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACAAATGCGAGTATCTCCCTTTTAAGTCTGTGGTTACGTCAAAAGAGAGCTTTGTCAGAATTATCTGCAAAGTCATACGCCTTTTAGTGctttg                    >  1:1017518/6‑148 (MQ=38)
                                                                                                                                tCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttct                   >  2:186623/1‑149 (MQ=255)
                                                                                                                                  aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAAGGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttcttt                 >  2:473363/1‑149 (MQ=255)
                                                                                                                                    aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGGGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCttatt               >  2:213887/1‑149 (MQ=255)
                                                                                                                                     aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt              <  1:533492/149‑1 (MQ=255)
                                                                                                                                     aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt              <  2:730129/149‑1 (MQ=255)
                                                                                                                                      aacaCCTCCCAGCGATGGCAGGGTTTTTTTTGCCTGTCGCTAGCGCACACCAATGCGAGTATCTCCCTAATAAGTCTGGGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg             <  2:666749/146‑1 (MQ=255)
                                                                                                                                      aaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg             <  1:928876/149‑1 (MQ=255)
                                                                                                                                       aaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg            >  2:552874/1‑149 (MQ=255)
                                                                                                                                        aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACATTTATCGGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc           <  1:745812/149‑1 (MQ=255)
                                                                                                                                        aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc           <  2:392829/149‑1 (MQ=255)
                                                                                                                                        aaCCTGCCAGCGATGGCAGGTTTTTTATTGCCTGTCGAGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc           <  1:849480/149‑1 (MQ=255)
                                                                                                                                             gCCAGCGATGGCAGGTTTTTTTTTGCCTGGCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaag      <  1:354481/149‑1 (MQ=255)
                                                                                                                                             gCCAGCGAGGGCAGGTTTTTTTTTGCATGTCGTGAGCGCACACCAATGCGAGTAGCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaag      <  1:902358/149‑1 (MQ=255)
                                                                                                                                              ccAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagc     <  1:1026038/149‑1 (MQ=255)
                                                                                                                                                aGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgt   <  1:241263/149‑1 (MQ=255)
                                                                                                                                                 gCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgta  <  1:264538/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTA  >  NC_000913/4060116‑4060409

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: