Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,160 G→C 5.1% 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.1% 959.5 / 33.4 269intergenic (‑143/+28)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (129/122);  new base C (6/8);  total (136/133)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 272 total aligned reads displayed.

CTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCA  >  NC_000913/2028012‑2028308
                                                                                                                                                    |                                                                                                                                                     
ctgcACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGa                                                                                                                                                      <  1:2258862/150‑1 (MQ=255)
      ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc                                                                                                                                                >  2:2692115/1‑150 (MQ=255)
      ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc                                                                                                                                                <  1:2871118/150‑1 (MQ=255)
      ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc                                                                                                                                                <  2:1261288/150‑1 (MQ=255)
      ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc                                                                                                                                                <  1:4294832/150‑1 (MQ=255)
      ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc                                                                                                                                                <  2:332192/150‑1 (MQ=255)
      ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGc                                                                                                                                                <  1:2207278/150‑1 (MQ=255)
       cTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                                                                               >  1:3306691/1‑150 (MQ=255)
       cTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTA‑‑ACATCAGGCCGGATCACATCCGGCCTg                                                                                                                                             <  2:802782/150‑1 (MQ=255)
         gacgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTg                                                                                                                                             <  2:1366828/150‑1 (MQ=255)
         gacgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTg                                                                                                                                             <  1:3775120/150‑1 (MQ=255)
           cgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                                                                           <  1:1765269/150‑1 (MQ=255)
            gacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                                                                          >  1:3301434/1‑150 (MQ=255)
             acTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGt                                                                                                                                         <  1:4906211/150‑1 (MQ=255)
              cTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtg                                                                                                                                        <  1:2644914/150‑2 (MQ=255)
               tCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtgc                                                                                                                                       <  1:3978845/150‑3 (MQ=255)
               tCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTa                                                                                                                                       >  1:797564/1‑150 (MQ=255)
                 ccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                                                                     >  2:3168768/1‑150 (MQ=255)
                 ccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                                                                     <  1:4329589/150‑1 (MQ=255)
                   aaCCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGa                                                                                                                                   <  1:1652036/150‑1 (MQ=255)
                       tGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAcc                                                                                                                               <  2:1980732/150‑1 (MQ=255)
                           tttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCttt                                                                                                                           >  2:2886647/1‑150 (MQ=255)
                           tttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCttt                                                                                                                           >  1:1932392/1‑150 (MQ=255)
                             tttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTca                                                                                                                         <  2:4194986/150‑1 (MQ=255)
                               tATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcaca                                                                                                                       >  2:3884428/1‑150 (MQ=255)
                                 tGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct                                                                                                                     >  1:1077238/1‑150 (MQ=255)
                                  gATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                                                                    >  2:3338018/1‑150 (MQ=255)
                                  gATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                                                                    <  1:1063588/150‑1 (MQ=255)
                                    ttCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActtct                                                                                                                  >  1:255670/1‑150 (MQ=255)
                                      cTGGCATAACCTGGCGATAGCGTACCGCAACACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCTTTCTTTGCAGGCATTAGCAGACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACttcttc                                                                                                                <  2:4306338/150‑1 (MQ=255)
                                            tAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGa                                                                                                          <  2:4017904/150‑1 (MQ=255)
                                            tAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGa                                                                                                          >  2:1681427/1‑150 (MQ=255)
                                             aaCCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATGCCCTTTCACACTTCTTCCGCAGaa                                                                                                         <  2:70799/150‑1 (MQ=255)
                                                cTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGa                                                                                                      >  1:476890/1‑150 (MQ=255)
                                                 tGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAc                                                                                                     <  2:4399712/150‑1 (MQ=255)
                                                   gCGATAGCGTACCGCAAACCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActc                                                                                                   <  2:109554/150‑1 (MQ=255)
                                                     gATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActccc                                                                                                 >  1:939652/1‑150 (MQ=255)
                                                     gATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActccc                                                                                                 >  2:2643811/1‑150 (MQ=255)
                                                      aTAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccc                                                                                                <  2:370828/150‑1 (MQ=255)
                                                          cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccct‑ccc                                                                                           >  1:54707/1‑150 (MQ=255)
                                                          cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccct‑ccc                                                                                           >  1:2193467/1‑150 (MQ=255)
                                                               cGCAAACCAGCGATATTTTTATAATGGTGGCATGATAAGACGCGTCTATAATTCTTTGGAGGAATTAGCACACCTGGCCGGTTGTGATCGGGCCTGATGTTTATTATAACCTTTCACTCTTCTTCCGCAGAACGACTCCCCT‑CCCCCggt                                                                                      <  2:2293741/150‑2 (MQ=255)
                                                               cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGa                                                                                      >  1:1900811/1‑150 (MQ=255)
                                                                        gCGATTATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCg                                                                             <  2:2066509/150‑1 (MQ=255)
                                                                         cGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCgg                                                                            >  1:311775/1‑150 (MQ=255)
                                                                          gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGt                                                                           >  1:4848426/1‑150 (MQ=255)
                                                                           aTAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGtt                                                                          >  1:3710971/1‑150 (MQ=255)
                                                                           aTAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGtt                                                                          >  1:1003941/1‑150 (MQ=255)
                                                                            tAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGttt                                                                         <  2:909984/150‑1 (MQ=255)
                                                                               ttttATAATGCTGGCATTATAAGACGAGACTATCTTTCATTGCATGAATTAGCACACCAGGCCGGATTTTATCCGTCCTGTTGTTAATGATACCCTTTCACACTTATTCGTAAGAACGACTCCACT‑CAACCGTAAGAAGCTCGGTTTCAt                                                                      <  2:4378236/150‑1 (MQ=255)
                                                                                tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAACCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATg                                                                     >  2:5111533/1‑150 (MQ=255)
                                                                                 ttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACGCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGa                                                                    >  1:3832980/1‑150 (MQ=255)
                                                                                  tataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATTTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGaa                                                                   <  2:1563607/150‑1 (MQ=255)
                                                                                  tataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGaa                                                                   >  1:1597618/1‑150 (MQ=255)
                                                                                   ataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGaaa                                                                  <  1:2376152/150‑1 (MQ=255)
                                                                                   ataATGGTGCAATGAAAACACGCTAATAACATTCATTGCAGCCATTACCACACCAGGCCGCATGTCATCCGCCCTGATGTTAATGATACCCTTTCACACTTCTACCGCAGAACGACTCCCCT‑CCCCCCGAACAAGCTCGGTTTCATGaaa                                                                  <  1:3396362/150‑1 (MQ=255)
                                                                                        ggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCtt                                                             >  2:2718194/1‑150 (MQ=255)
                                                                                        ggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTGTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCtt                                                             <  2:4557062/150‑1 (MQ=255)
                                                                                              aTGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt                                                       >  2:3274015/1‑150 (MQ=255)
                                                                                              aTGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt                                                       <  2:1878339/150‑1 (MQ=255)
                                                                                                  tAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAc                                                   <  2:1295562/150‑1 (MQ=255)
                                                                                                   aaGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcc                                                  <  2:476890/150‑1 (MQ=255)
                                                                                                   aaGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcc                                                  <  2:98466/150‑1 (MQ=255)
                                                                                                   aaGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcc                                                  >  1:3772565/1‑150 (MQ=255)
                                                                                                   aaGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcc                                                  <  2:1616289/150‑1 (MQ=255)
                                                                                                   aaGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACATCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcc                                                  >  1:2153757/1‑150 (MQ=255)
                                                                                                   aaGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGACGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcc                                                  >  2:379493/1‑150 (MQ=255)
                                                                                                     gACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcccc                                                >  1:1889008/1‑150 (MQ=255)
                                                                                                          gACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAAc                                           >  1:2329516/1‑150 (MQ=255)
                                                                                                          gACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAAc                                           >  2:117480/1‑150 (MQ=255)
                                                                                                             tATCATTCATTGCAGGCATTAGCACACCAGGCCGGCTGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGTACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAACCCGTAACCGt                                        >  2:1415991/1‑150 (MQ=255)
                                                                                                              aTCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGt                                        <  2:539325/150‑1 (MQ=255)
                                                                                                               tcattcatTGCAGGCATTAACACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACCCTTCGTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                      <  1:2014243/150‑1 (MQ=255)
                                                                                                                cattcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCAt                                     >  2:4586424/1‑150 (MQ=255)
                                                                                                                 attcattGCAGGCATTAGCACACCAGGCGGGTTGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATa                                    <  1:1746206/150‑1 (MQ=255)
                                                                                                                 attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATa                                    <  2:4026695/150‑1 (MQ=255)
                                                                                                                  ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                   >  2:1031378/1‑150 (MQ=255)
                                                                                                                  ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                   <  2:41593/150‑1 (MQ=255)
                                                                                                                  ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                   >  1:2745669/1‑150 (MQ=255)
                                                                                                                   tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAcc                                  <  1:4594160/150‑1 (MQ=255)
                                                                                                                     attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCaa                                <  2:2138187/150‑1 (MQ=255)
                                                                                                                      ttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAAc                               <  1:526309/150‑1 (MQ=255)
                                                                                                                        gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGa                             >  1:1408137/1‑150 (MQ=255)
                                                                                                                          tggaaTTAGAACACCAGGCGGGTTGTTATCGGGCATGATTTTATTTATACCATTTAACACTTCTTCCGAAGAACGACTCCCCT‑CCCCGGTAACATGCTCGGTTTATTTAAAAGCTTCCCGCTTGACCCCGTAACCGTAATACCAACGaca                           <  1:1369400/146‑1 (MQ=255)
                                                                                                                          aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTTATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGaca                           <  1:272372/150‑1 (MQ=255)
                                                                                                                          aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGaca                           >  1:4797782/1‑150 (MQ=255)
                                                                                                                            gCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATTTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACt                         <  2:3493099/150‑1 (MQ=255)
                                                                                                                            gCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACt                         >  1:3895019/1‑150 (MQ=255)
                                                                                                                             cATTAGCACACCTGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTc                        <  2:636156/150‑1 (MQ=255)
                                                                                                                              aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCa                       >  1:1142404/1‑150 (MQ=255)
                                                                                                                               ttAGCACACCAGGCCGGATGTGATCCGGCCTTATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCaa                      <  2:2276669/150‑1 (MQ=255)
                                                                                                                                        cAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGc             >  2:4832611/1‑150 (MQ=255)
                                                                                                                                         aGGCCGGATGTGATCCGGCCTGATTTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCt            <  2:2468056/150‑1 (MQ=255)
                                                                                                                                          ggCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTg           >  1:1216080/1‑150 (MQ=255)
                                                                                                                                            ccGGATTTGATCCGGCCTGTTGTTAATGATACCCTTACACACTCCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCTTAACCGTCATACCAACGACACTCAACCATACCGCTGGa         <  1:4428259/150‑1 (MQ=255)
                                                                                                                                            ccGGATTTGATCCGGCCTGATTTTAATGATACACTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGa         <  2:525447/150‑1 (MQ=255)
                                                                                                                                            ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGa         <  1:1795919/150‑1 (MQ=255)
                                                                                                                                               gATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTa      >  1:4029889/1‑150 (MQ=255)
                                                                                                                                                aTGTGATCCGGCCTAATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTAt     >  2:3801785/1‑150 (MQ=255)
                                                                                                                                                 tgtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATc    <  2:2569107/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCT‑CCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCA  >  NC_000913/2028012‑2028308

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: ‑