Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,160 G→C 5.7% intergenic (‑143/+28) 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,1600GC 5.7% 682.3 / 24.1 190intergenic (‑143/+28)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (84/94);  new base C (7/4);  total (91/99)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

GCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAG  >  NC_000913/2028011‑2028309
                                                                                                                                                     |                                                                                                                                                       
gctgcACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtgtg                                                                                                                                                         >  2:3725171/1‑150 (MQ=255)
  tgcACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                                                                       <  1:3300549/150‑1 (MQ=255)
    cACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                                                                     <  1:2465973/150‑1 (MQ=255)
    cACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                                                                     >  1:2182264/1‑150 (MQ=255)
     aCCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCg                                                                                                                                                    <  1:1778656/150‑1 (MQ=255)
      cccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                                                                                   >  1:2070021/1‑150 (MQ=255)
      cccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCAGAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGTTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                                                                                   <  1:3238806/150‑1 (MQ=255)
       ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc                                                                                                                                                  <  1:795329/150‑1 (MQ=255)
       ccTGACGACTCACAACCTGCTTTTTTATGATTCTGGCATAACAGGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGc                                                                                                                                                  <  2:2442552/150‑1 (MQ=255)
        cTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                                                                                 <  1:3765610/150‑1 (MQ=255)
          gacgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCTATAATTTTATATTGGTGGCATGATAAGACGCGACTATAATTCATTGCAGGAATTAGCACACCTGGCCGGATGTGATCCGGCCTg                                                                                                                                               <  1:3546811/150‑1 (MQ=255)
          gacgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg                                                                                                                                               <  2:2379202/150‑1 (MQ=255)
                    aaCCTGCTTTTTTATGATTCTGGCATACCCTGGCGATAGCGTACCGCATCACAGCGCTAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGa                                                                                                                                   <  2:3662145/150‑1 (MQ=255)
                    aaCCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGa                                                                                                                                   <  1:179071/150‑1 (MQ=255)
                           ttttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCtt                                                                                                                            >  1:125867/1‑150 (MQ=255)
                            tttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCttt                                                                                                                           <  2:1135501/150‑1 (MQ=255)
                              tttATGATTGTGGCATATGCTGGCGTTAGCGTACCGCATAACAGCGATATTTTTATAATGGTGGCATGATAAGACGCGACTATCATTGATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTca                                                                                                                         <  2:2540912/150‑1 (MQ=255)
                              tttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAAGGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCAGTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTca                                                                                                                         >  2:2174488/1‑150 (MQ=255)
                                  tGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACAct                                                                                                                     >  2:1282092/1‑150 (MQ=255)
                                    aTTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACActtc                                                                                                                   <  2:2935856/150‑1 (MQ=255)
                                       cTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACttcttc                                                                                                                <  1:2860045/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                                                                                                            >  2:437228/1‑150 (MQ=255)
                                           cATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCa                                                                                                            >  1:604741/1‑150 (MQ=255)
                                                ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACg                                                                                                       >  2:2771914/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)
                                                       aTAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGGTGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccc                                                                                                >  2:459374/1‑150 (MQ=255)
                                                       aTAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccc                                                                                                >  1:549087/1‑150 (MQ=255)
                                                       aTAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccc                                                                                                >  1:2679298/1‑150 (MQ=255)
                                                            gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACtcccctcccc                                                                                           >  1:495685/1‑150 (MQ=255)
                                                            gTACCGCAAAACAGCGATAATTTAATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCAGATGA‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACtccccacccc                                                                                           >  1:4130114/1‑150 (MQ=255)
                                                             tACCGCAAAACAGCGATAATTTTATAATGGTGGTATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTccccc                                                                                          >  1:859555/1‑150 (MQ=255)
                                                                cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                                                                       >  1:2331424/1‑150 (MQ=255)
                                                                cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                                                                       >  1:3803399/1‑150 (MQ=255)
                                                                  cAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc                                                                                     >  1:3948909/1‑150 (MQ=255)
                                                                    aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACaa                                                                                   <  2:674097/150‑1 (MQ=255)
                                                                    aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACaa                                                                                   <  2:2297474/150‑1 (MQ=255)
                                                                       cAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                                                                                <  2:3590663/150‑1 (MQ=255)
                                                                        aGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTc                                                                               >  1:4006668/1‑150 (MQ=255)
                                                                         gCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg                                                                              >  2:1006670/1‑150 (MQ=255)
                                                                          cGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                                                                             >  1:972988/1‑150 (MQ=255)
                                                                           gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGt                                                                            >  1:156209/1‑150 (MQ=255)
                                                                              aaTTTTATAATGGTGGCATTATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTc                                                                         <  1:3836586/150‑1 (MQ=255)
                                                                               aTTTTATAATGGTGGCATTATAAGACGCGACTATCATTCATTGCAGGAATTAGCACACCAGGCCGGATGTGATCCGGCCTTATTT‑‑TAATTATACCCTTTCACACTTCTTCCGCACAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                        <  1:1294512/150‑1 (MQ=255)
                                                                               aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCTGGCCGGTTGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                        <  1:2947464/150‑1 (MQ=255)
                                                                               aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                        <  1:4028724/150‑1 (MQ=255)
                                                                               aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                        <  1:667190/150‑1 (MQ=255)
                                                                               aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTGTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                        <  2:1594084/150‑1 (MQ=255)
                                                                               aTTTTATAATGGTGGAATGATAACACGCGACTATAATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                        <  1:3087520/150‑1 (MQ=255)
                                                                               aTTTTATAATGGGGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATTATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                        <  1:858444/150‑1 (MQ=255)
                                                                                   tataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                                                                    >  2:2160118/1‑150 (MQ=255)
                                                                                     taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa                                                                  <  2:2036112/150‑1 (MQ=255)
                                                                                     taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa                                                                  >  2:1489505/1‑150 (MQ=255)
                                                                                     taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa                                                                  <  2:3000429/150‑1 (MQ=255)
                                                                                      aaTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg                                                                 >  1:2904376/1‑150 (MQ=255)
                                                                                        tggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt                                                               <  1:3701682/150‑1 (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)
                                                                                             gCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCATTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCg                                                          <  2:1231798/150‑1 (MQ=255)
                                                                                               aTGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt                                                        >  1:1985004/1‑150 (MQ=255)
                                                                                                tGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCtt                                                       <  2:950857/150‑1 (MQ=255)
                                                                                                  aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa                                                     <  2:3322728/150‑1 (MQ=255)
                                                                                                  aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa                                                     <  2:380395/150‑1 (MQ=255)
                                                                                                    aaGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcc                                                   <  2:1332659/150‑1 (MQ=255)
                                                                                                    aaGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcc                                                   <  2:125893/150‑1 (MQ=255)
                                                                                                       aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGGGATCCGGCCAGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCGTGACCCCg                                                >  2:771660/1‑150 (MQ=255)
                                                                                                           gACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAAc                                            >  2:1515137/1‑150 (MQ=255)
                                                                                                              tATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGt                                         <  2:1684747/150‑1 (MQ=255)
                                                                                                                tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                       >  2:3540798/1‑150 (MQ=255)
                                                                                                                tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                       <  2:3433719/150‑1 (MQ=255)
                                                                                                                tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                       >  1:3498623/1‑150 (MQ=255)
                                                                                                                tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                       <  1:2768779/150‑1 (MQ=255)
                                                                                                                tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                       >  2:1386239/1‑150 (MQ=255)
                                                                                                                tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCATGATGT‑‑TAATGATACCCATTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                       >  2:3114088/1‑150 (MQ=255)
                                                                                                                tcattcatTGCAGGAATTAGCACACCAGGCCGGATGTGATCGGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                       <  1:1874451/150‑1 (MQ=255)
                                                                                                                  attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATa                                     >  2:2664813/1‑150 (MQ=255)
                                                                                                                   ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                    <  2:4050057/150‑1 (MQ=255)
                                                                                                                   ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                    >  1:1181812/1‑150 (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:3885029/1‑150 (MQ=255)
                                                                                                                     cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCa                                  >  2:1346615/1‑150 (MQ=255)
                                                                                                                      attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCaa                                 >  2:2065186/1‑150 (MQ=255)
                                                                                                                       ttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAAc                                <  1:3138736/150‑1 (MQ=255)
                                                                                                                             gCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACt                          >  2:87350/1‑150 (MQ=255)
                                                                                                                              cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTc                         <  2:4088333/150‑1 (MQ=255)
                                                                                                                               aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCa                        >  1:1652614/1‑150 (MQ=255)
                                                                                                                                ttAGCACACCTGGCGGGTTGTGATCCGGCCTGATTT‑‑TATTGTTACCCTTTCACACTTCTTCCGAAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTATTGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCaa                       <  1:3725171/150‑1 (MQ=255)
                                                                                                                                       atcAGGCCGGATCACATCCGGCCTGGTGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCAGAcc                >  2:3858720/3‑150 (MQ=255)
                                                                                                                                        tcAGGCCGGATCACATCCGGCCTGGTGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCg               >  2:882012/2‑150 (MQ=255)
                                                                                                                                           ggCCGGATCACATCCGGCCTGGTGTGCTAATGATACCCTGTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGc              >  1:503908/1‑150 (MQ=255)
                                                                                                                                                  tgtgATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATc     >  1:2441400/1‑150 (MQ=255)
                                                                                                                                                     gATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAg  <  1:3620414/150‑1 (MQ=255)
                                                                                                                                                     gATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAg  <  2:3948910/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
GCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAG  >  NC_000913/2028011‑2028309

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: