Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,158 G→C 5.8% intergenic (‑141/+30) yedQ ← / ← yodC putative membrane‑anchored diguanylate cyclase/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,028,1580GC 5.8% 692.0 / 27.6 190intergenic (‑141/+30)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (84/95);  new base C (7/4);  total (91/99)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

Only 100 of 193 total aligned reads displayed.

GTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATC  >  NC_000913/2028009‑2028306
                                                                                                                                                     |                                                                                                                                                      
  gctgcACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtgtg                                                                                                                                                      >  2:3725171/1‑150 (MQ=255)
    tgcACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                                                                    <  1:3300549/150‑1 (MQ=255)
      cACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                                                                  <  1:2465563/150‑1 (MQ=255)
      cACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                                                                  >  1:2182264/1‑150 (MQ=255)
       aCCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCg                                                                                                                                                 <  1:3681993/150‑1 (MQ=255)
       aCCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCg                                                                                                                                                 <  1:1778656/150‑1 (MQ=255)
        cccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                                                                                >  1:2070021/1‑150 (MQ=255)
        cccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCAGAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGTTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                                                                                <  1:3238806/150‑1 (MQ=255)
         ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACATCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGc                                                                                                                                               <  2:1192776/150‑1 (MQ=255)
         ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc                                                                                                                                               <  1:795329/150‑1 (MQ=255)
         ccTGACGACTCACAACCTGCTTTTTTATGATTCTGGCATAACAGGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGc                                                                                                                                               <  2:2442552/150‑1 (MQ=255)
          cTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGcc                                                                                                                                              <  1:3733564/150‑1 (MQ=255)
            gacgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCTATAATTTTATATTGGTGGCATGATAAGACGCGACTATAATTCATTGCAGGAATTAGCACACCTGGCCGGATGTGATCCGGCCTg                                                                                                                                            <  1:3546811/150‑1 (MQ=255)
            gacgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg                                                                                                                                            <  2:2379202/150‑1 (MQ=255)
               gacTCCCACCCTGCTTTTTGATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGCTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCCGGCCGGATGTGATCCGGCCTGATg                                                                                                                                         <  2:2296314/150‑1 (MQ=255)
                   cccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑Taa                                                                                                                                   <  1:3945345/150‑1 (MQ=255)
                      aaCCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGa                                                                                                                                <  1:179071/150‑1 (MQ=255)
                             ttttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCtt                                                                                                                         >  1:125867/1‑150 (MQ=255)
                              tttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCttt                                                                                                                        <  2:1106919/150‑1 (MQ=255)
                              tttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCttt                                                                                                                        >  1:3723700/1‑150 (MQ=255)
                                tttATGATTGTGGCATATGCTGGCGTTAGCGTACCGCATAACAGCGATATTTTTATAATGGTGGCATGATAAGACGCGACTATCATTGATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTca                                                                                                                      <  2:2540912/150‑1 (MQ=255)
                                  tATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTcaca                                                                                                                    <  1:866724/150‑1 (MQ=255)
                                    tGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACAct                                                                                                                  >  2:1282092/1‑150 (MQ=255)
                                     gATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACActt                                                                                                                 <  2:2168510/150‑1 (MQ=255)
                                     gATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACActt                                                                                                                 >  2:2044488/1‑150 (MQ=255)
                                        tCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACActtctt                                                                                                              <  2:302097/150‑1 (MQ=255)
                                         cTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACttcttc                                                                                                             <  1:2860045/150‑1 (MQ=255)
                                          tGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTcc                                                                                                            <  1:3793486/150‑1 (MQ=255)
                                          tGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTcc                                                                                                            <  1:2876274/150‑1 (MQ=255)
                                           ggCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCg                                                                                                           <  2:1920039/150‑1 (MQ=255)
                                            gCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGc                                                                                                          >  1:1125690/1‑150 (MQ=255)
                                             cATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCa                                                                                                         >  1:604741/1‑150 (MQ=255)
                                             cATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCa                                                                                                         >  2:437228/1‑150 (MQ=255)
                                                aaCCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGaa                                                                                                      <  1:3133153/150‑1 (MQ=255)
                                                  ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACg                                                                                                    >  2:1280314/1‑150 (MQ=255)
                                                      gCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGActc                                                                                                >  1:2185307/1‑150 (MQ=255)
                                                       cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                                                                               <  2:1540097/150‑1 (MQ=255)
                                                       cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                                                                               >  1:2826521/1‑150 (MQ=255)
                                                        gATTGCGTACCGCAAAACAGCGATTATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGActccc                                                                                              <  2:3340473/150‑1 (MQ=255)
                                                         aTAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGGTGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccc                                                                                             >  2:459374/1‑150 (MQ=255)
                                                         aTAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccc                                                                                             >  1:4028382/1‑150 (MQ=255)
                                                         aTAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccc                                                                                             >  1:2679298/1‑150 (MQ=255)
                                                         aTAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccc                                                                                             >  1:549087/1‑150 (MQ=255)
                                                              gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACtcccctcccc                                                                                        >  1:495685/1‑150 (MQ=255)
                                                              gTACCGCAAAACAGCGATAATTTAATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCAGATGA‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACtccccacccc                                                                                        >  1:4130114/1‑150 (MQ=255)
                                                                  cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGAACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                                                                    >  2:3169694/1‑150 (MQ=255)
                                                                    cAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc                                                                                  >  1:3948910/1‑150 (MQ=255)
                                                                      aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACaa                                                                                <  2:2297474/150‑1 (MQ=255)
                                                                      aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACaa                                                                                <  2:674097/150‑1 (MQ=255)
                                                                             gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGt                                                                         >  1:156209/1‑150 (MQ=255)
                                                                             gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGt                                                                         >  2:3616793/1‑150 (MQ=255)
                                                                                 aTTTTATAATGGTGGCATTATAAGACGCGACTATCATTCATTGCAGGAATTAGCACACCAGGCCGGATGTGATCCGGCCTTATTT‑‑TAATTATACCCTTTCACACTTCTTCCGCACAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                     <  1:1294512/150‑1 (MQ=255)
                                                                                 aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCTGGCCGGTTGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                     <  1:2947464/150‑1 (MQ=255)
                                                                                 aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                     <  1:667190/150‑1 (MQ=255)
                                                                                 aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                     <  1:4028724/150‑1 (MQ=255)
                                                                                 aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTGTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                     <  2:1594084/150‑1 (MQ=255)
                                                                                 aTTTTATAATGGGGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATTATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                     <  1:858444/150‑1 (MQ=255)
                                                                                   tttATAATGGTGGCATGATAAGACGCGACTATCTTTCATTGCAGGCATTAGCACACCTGGCCGGATGTGATCCGGCCTGTTGT‑‑TATTGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                                                                   <  1:3814867/150‑1 (MQ=255)
                                                                                   tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                                                                   >  1:2553644/1‑150 (MQ=255)
                                                                                       taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa                                                               >  2:1489505/1‑150 (MQ=255)
                                                                                        aaTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg                                                              >  1:2904376/1‑150 (MQ=255)
                                                                                            gtggCATTATAAGACGCTACTATAATTCATTGCAGGCATTAGCACACCAGGCCGGATGTTATCCGGCCTTATTT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTc                                                          <  1:2054794/150‑1 (MQ=255)
                                                                                             tggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTcc                                                         >  2:615106/1‑150 (MQ=255)
                                                                                               gCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCg                                                       >  1:1385100/1‑150 (MQ=255)
                                                                                                 aTGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt                                                     >  1:1985004/1‑150 (MQ=255)
                                                                                                  tGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCtt                                                    <  2:354895/150‑1 (MQ=255)
                                                                                                  tGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCtt                                                    <  2:950857/150‑1 (MQ=255)
                                                                                                    aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa                                                  <  2:380395/150‑1 (MQ=255)
                                                                                                    aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa                                                  <  2:3723700/150‑1 (MQ=255)
                                                                                                      aaGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcc                                                <  2:1332659/150‑1 (MQ=255)
                                                                                                      aaGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcc                                                <  2:125867/150‑1 (MQ=255)
                                                                                                         aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGGGATCCGGCCAGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCGTGACCCCg                                             >  2:771660/1‑150 (MQ=255)
                                                                                                             gACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAAc                                         >  2:1515137/1‑150 (MQ=255)
                                                                                                              aCTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAAcc                                        <  2:3497771/150‑1 (MQ=255)
                                                                                                                tATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGt                                      <  2:1684747/150‑1 (MQ=255)
                                                                                                                  tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                    >  2:1386239/1‑150 (MQ=255)
                                                                                                                  tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                    <  1:2768779/150‑1 (MQ=255)
                                                                                                                  tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                    >  1:3498623/1‑150 (MQ=255)
                                                                                                                  tcattcatTGCAGGAATTAGCACACCAGGCCGGATGTGATCGGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                    <  1:1874451/150‑1 (MQ=255)
                                                                                                                    attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATa                                  >  2:2664813/1‑150 (MQ=255)
                                                                                                                     ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                 >  1:1181812/1‑150 (MQ=255)
                                                                                                                     ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                 <  1:2305493/150‑1 (MQ=255)
                                                                                                                     ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                 <  2:4050057/150‑1 (MQ=255)
                                                                                                                     ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                 <  1:3364302/150‑1 (MQ=255)
                                                                                                                       cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCa                               <  1:3017765/150‑1 (MQ=255)
                                                                                                                       cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCa                               <  1:3541567/150‑1 (MQ=255)
                                                                                                                       cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCa                               >  2:1346615/1‑150 (MQ=255)
                                                                                                                       cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCa                               >  2:3885029/1‑150 (MQ=255)
                                                                                                                        attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCaa                              >  2:2065218/1‑150 (MQ=255)
                                                                                                                         ttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAAc                             >  1:2386431/1‑150 (MQ=255)
                                                                                                                               gCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACt                       >  2:87350/1‑150 (MQ=255)
                                                                                                                                cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTc                      <  2:4088333/150‑1 (MQ=255)
                                                                                                                                 aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTCACCGTCATACCAACGACACGCa                     >  2:2785918/1‑150 (MQ=255)
                                                                                                                                 aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCa                     >  1:1652614/1‑150 (MQ=255)
                                                                                                                                 aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCa                     >  1:3236363/1‑150 (MQ=255)
                                                                                                                                         atcAGGCCGGATCACATCCGGCCTGGTGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCAGAcc             >  2:3858720/3‑150 (MQ=255)
                                                                                                                                          tcAGGCCGGATCACATCCGGCCTGGTGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCg            >  2:2999018/2‑150 (MQ=255)
                                                                                                                                          tcAGGCCGGATCACATCCGGCCTGGTGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCg            >  2:882012/2‑150 (MQ=255)
                                                                                                                                            aGGCCGGATCACATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCt          >  1:996199/1‑150 (MQ=255)
                                                                                                                                             ggCCGGATCACATCCGGCCTGGTGTGCTAATGATACCCTGTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGc           >  1:503908/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATC  >  NC_000913/2028009‑2028306

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: