Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,159 T→A 5.3% intergenic (‑142/+29) 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,1590TA 5.3% 684.7 / 14.3 190intergenic (‑142/+29)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (83/95);  new base A (7/4);  total (91/99)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

GCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATC  >  NC_000913/2028011‑2028306
                                                                                                                                                    |                                                                                                                                                     
    cACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                                                                  <  1:2465563/150‑1 (MQ=255)
    cACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                                                                  <  2:1299801/150‑1 (MQ=255)
    cACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                                                                  >  1:2182269/1‑150 (MQ=255)
     aCCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCg                                                                                                                                                 <  1:42246/150‑1 (MQ=255)
     aCCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCg                                                                                                                                                 <  1:1778656/150‑1 (MQ=255)
      cccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                                                                                >  1:2070021/1‑150 (MQ=255)
       ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACATCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGc                                                                                                                                               <  2:1192776/150‑1 (MQ=255)
        cTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                                                                              <  1:3765610/150‑1 (MQ=255)
            cgacTCCGAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCTCAACACAGCGATAATTTTATACTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCAGTAGCACACCAGGCCGGCTGTGATCCGGCCTGAt                                                                                                                                          <  2:2299109/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)
                    aaCCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACACCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGa                                                                                                                                <  1:3266145/150‑1 (MQ=255)
                           ttttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCtt                                                                                                                         >  1:125867/1‑150 (MQ=255)
                            tttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCttt                                                                                                                        <  2:1106919/150‑1 (MQ=255)
                                tATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTcaca                                                                                                                    <  1:866724/150‑1 (MQ=255)
                                 aTGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCacac                                                                                                                   <  2:1343267/150‑1 (MQ=255)
                                 aTGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCacac                                                                                                                   <  1:532058/150‑1 (MQ=255)
                                   gATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACActt                                                                                                                 <  2:2168510/150‑1 (MQ=255)
                                   gATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACActt                                                                                                                 >  2:2044488/1‑150 (MQ=255)
                                     ttCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACActtct                                                                                                               >  1:2854365/1‑150 (MQ=255)
                                      tCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACActtctt                                                                                                              <  2:3848871/150‑1 (MQ=255)
                                      tCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACActtctt                                                                                                              <  2:1513552/150‑1 (MQ=255)
                                      tCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACActtctt                                                                                                              <  2:302097/150‑1 (MQ=255)
                                       cTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACttcttc                                                                                                             <  1:3306674/150‑1 (MQ=255)
                                        tGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTcc                                                                                                            <  1:3793486/150‑1 (MQ=255)
                                         ggCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCg                                                                                                           <  2:1920039/150‑1 (MQ=255)
                                          gCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGc                                                                                                          >  1:1125690/1‑150 (MQ=255)
                                           catatacTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATATTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACTCCAGGCCGGATGTGATCCGTACTGATGT‑‑TATTGTTAGCCTTTCACACTTCTTCCGCa                                                                                                         <  2:1439077/144‑1 (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)
                                                     cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                                                                               <  2:1540097/150‑1 (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)
                                                              aCCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAAGT‑‑TAATGATACCCTTTCACACTACTTCCGCAGAACGACTCCCCTCCCCCg                                                                                      >  2:2507458/1‑150 (MQ=255)
                                                                cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                                                                    >  1:2975109/1‑150 (MQ=255)
                                                                cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                                                                    >  1:4084494/1‑150 (MQ=255)
                                                                cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                                                                    >  1:2807854/1‑150 (MQ=255)
                                                                cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                                                                    >  1:2697434/1‑150 (MQ=255)
                                                                cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGAACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                                                                    >  2:3169694/1‑150 (MQ=255)
                                                                 gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGaa                                                                                   >  2:3742401/1‑150 (MQ=255)
                                                                  cAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc                                                                                  >  1:3492187/1‑150 (MQ=255)
                                                                  cAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc                                                                                  >  1:3948910/1‑150 (MQ=255)
                                                                    aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACaa                                                                                <  2:674097/150‑1 (MQ=255)
                                                                     aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑AAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAg                                                                               >  2:2560124/1‑150 (MQ=255)
                                                                           gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGt                                                                         >  2:3616793/1‑150 (MQ=255)
                                                                           gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGt                                                                         >  1:156209/1‑150 (MQ=255)
                                                                            aTAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGtt                                                                        >  2:1382814/1‑150 (MQ=255)
                                                                               aTTTTATAATGGTGGCATTATAAGACGCGACTATCATTCATTGCAGGAATTAGCACACCAGGCCGGATGTGATCCGGCCTTATTT‑‑TAATTATACCCTTTCACACTTCTTCCGCACAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                                     <  1:1294512/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)
                                                                                 tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                                                                   >  1:853587/1‑150 (MQ=255)
                                                                                   tataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGAAACAAGCTCGGTTTCATGaa                                                                 >  1:2875680/1‑150 (MQ=255)
                                                                                      aaTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg                                                              >  1:2904367/1‑150 (MQ=255)
                                                                                        tggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt                                                            <  1:3701682/150‑1 (MQ=255)
                                                                                        tggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt                                                            <  1:761551/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)
                                                                                               attaTAAGACGCGACTATCATTCATTGCAGGAATTAGCACACCAGGCCGTATGTGATCCGGCCTTATGT‑‑TAATGATACCCTTTCAAACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt                                                     <  1:2923617/147‑1 (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:3723700/150‑1 (MQ=255)
                                                                                                    aaGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcc                                                <  2:125867/150‑1 (MQ=255)
                                                                                                       aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCg                                             >  1:612355/1‑150 (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)
                                                                                                                tcattcatTGCAGGAATTAGCACACCAGGCCGGATGTGATCGGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                    <  1:1874451/150‑1 (MQ=255)
                                                                                                                  attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATa                                  <  1:3541467/150‑1 (MQ=255)
                                                                                                                  attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATa                                  >  2:2664813/1‑150 (MQ=255)
                                                                                                                  attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATa                                  >  2:2421755/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)
                                                                                                                     cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCa                               >  2:3885029/1‑150 (MQ=255)
                                                                                                                     cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCa                               <  1:3017765/150‑1 (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:3236363/1‑150 (MQ=255)
                                                                                                                               aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCa                     >  1:1652614/1‑150 (MQ=255)
                                                                                                                                ttAGCACACCTGGCGGGTTGTGATCCGGCCTGATTT‑‑TATTGTTACCCTTTCACACTTCTTCCGAAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTATTGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCaa                    <  1:3725171/150‑1 (MQ=255)
                                                                                                                                   gCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCa                 >  1:2181237/1‑150 (MQ=255)
                                                                                                                                    cacaCCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCAt                >  1:3607100/1‑150 (MQ=255)
                                                                                                                                      cacCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATAc              <  1:3544581/150‑1 (MQ=255)
                                                                                                                                       atcAGGCCGGATCACATCCGGCCTGGTGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCAGAcc             >  2:3858720/3‑150 (MQ=255)
                                                                                                                                        tcAGGCCGGATCACATCCGGCCTGGTGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCg            >  2:882012/2‑150 (MQ=255)
                                                                                                                                        tcAGGCCGGATCACATCCGGCCTGGTGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCg            >  2:2999018/2‑150 (MQ=255)
                                                                                                                                          aGGCCGGATCACATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCt          >  1:996199/1‑150 (MQ=255)
                                                                                                                                           ggCCGGATCACATCCGGCCTGGTGTGCTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGc           >  1:503875/1‑150 (MQ=255)
                                                                                                                                           ggCCGGATCACATCCGGCCTGGTGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTg         >  1:575717/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGT‑‑TAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATC  >  NC_000913/2028011‑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: