Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,159 T→A 15.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,1590TA15.3% 306.6 / 45.0 118intergenic (‑142/+29)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (49/51);  new base A (7/11);  total (56/62)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01

Only 100 of 123 total aligned reads displayed.

GCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGA  >  NC_000913/2028063‑2028258
                                                                                                |                                                                                                   
gCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                  <  2:3350418/100‑1 (MQ=255)
 cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc                                                                                                 <  2:2394590/100‑1 (MQ=255)
 cGATAGCGTACCGCAAAACAGAAATAATTTTATAATGGTGGCATGATAAGACGCAACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                  >  1:1233364/1‑99 (MQ=255)
   aTAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCg                                                                                               <  2:2894759/100‑1 (MQ=255)
    tAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                              <  2:2166184/100‑1 (MQ=255)
    tAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                              <  1:1592626/100‑1 (MQ=255)
      gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                            >  2:3260918/1‑100 (MQ=255)
       cGTACCGTAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                           <  1:1045430/100‑1 (MQ=255)
       cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                           <  1:943584/100‑1 (MQ=255)
        gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg                                                                                          >  1:332264/1‑100 (MQ=255)
        gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTg                                                                                          <  2:406834/100‑1 (MQ=255)
         tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                         >  1:157405/1‑100 (MQ=255)
         tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                         <  2:293149/100‑1 (MQ=255)
          aCCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                        <  1:2488893/100‑1 (MQ=255)
           ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                       >  2:1911576/1‑100 (MQ=255)
           ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                       >  2:238469/1‑100 (MQ=255)
            cGCACAACAGCGATAATTTCATAAAGGAGGCACCATAAGACGCCACTATCATTCATTACAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCCGATg                                                                                       >  1:1498744/1‑99 (MQ=255)
            cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGt                                                                                      >  1:2231290/1‑100 (MQ=255)
            cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGt                                                                                      <  1:307424/100‑1 (MQ=255)
            cGCAAAACAGCGATAATTTGAGAATGGTGGCATGTTAAGACGAGACTAACATCCACCGCCGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                           <  1:170746/95‑1 (MQ=255)
             gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtg                                                                                     <  2:2277861/100‑2 (MQ=37)
             gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtg                                                                                     <  2:2277860/100‑2 (MQ=37)
             gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtg                                                                                     <  1:2900003/100‑2 (MQ=37)
             gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGtt                                                                                     <  2:1058864/100‑1 (MQ=255)
                 aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                 <  1:984164/100‑1 (MQ=255)
                 aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                 <  2:157405/100‑1 (MQ=255)
                 aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                 <  2:181783/100‑1 (MQ=255)
                    aGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtgctaatgc                                                                              <  1:1831067/100‑9 (MQ=21)
                     gCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtgctaatgcc                                                                             <  2:2758845/100‑10 (MQ=21)
                       gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATATTAATGATAccc                                                                           >  1:3064904/1‑100 (MQ=255)
                           aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                       <  2:1132938/100‑1 (MQ=255)
                           aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                       <  1:1575369/100‑1 (MQ=255)
                           aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTAATTGCAGGCATTAGCACACCAGGCCGGCTGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                       >  1:1074064/1‑100 (MQ=255)
                             tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac                                                                     <  2:1649609/100‑1 (MQ=255)
                             tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac                                                                     <  1:342987/100‑1 (MQ=255)
                             tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac                                                                     <  1:3247160/100‑1 (MQ=255)
                             tttAGAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac                                                                     <  2:3241251/100‑1 (MQ=255)
                               tataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCacac                                                                   <  2:1791325/100‑1 (MQ=255)
                                ataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct                                                                  <  2:1234079/100‑1 (MQ=255)
                                 taATGGTGGCATGATAAGTCGCGACTATCATTCATTGCAGGCATTAGCACACTAGGCCGGACGTTATCCGGCCTGATGTTAATGATACCCTTTCACAgtt                                                                 >  1:10472/1‑100 (MQ=255)
                                 taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                 >  2:1869137/1‑100 (MQ=255)
                                 cttccgatctaTGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                 <  2:590296/90‑1 (MQ=255)
                                    tggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActtctt                                                              <  1:1927745/100‑1 (MQ=255)
                                     ggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACttcttc                                                             <  2:533325/100‑1 (MQ=255)
                                        ggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc                                                          >  2:1346905/1‑100 (MQ=255)
                                          cATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAg                                                        >  1:207418/1‑100 (MQ=255)
                                          cATGATAAGACGAGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCa                                                         >  1:2309026/1‑99 (MQ=255)
                                           aTGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActtagatcggaag                                                       >  1:590293/1‑90 (MQ=255)
                                           aTGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGa                                                       <  2:3153065/100‑1 (MQ=255)
                                            tGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGaa                                                      >  2:1969803/1‑100 (MQ=255)
                                             gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                     >  1:1251836/1‑100 (MQ=255)
                                             gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                     <  2:2990944/100‑1 (MQ=255)
                                              aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACg                                                    <  2:443175/100‑1 (MQ=255)
                                              aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACg                                                    <  1:3283231/100‑1 (MQ=255)
                                              aTAAGACGCGACTATCATTCATCGCAGGCATCAGCACCCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACg                                                    <  1:2282310/100‑1 (MQ=255)
                                               tAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGa                                                   >  2:2150774/1‑100 (MQ=255)
                                                 aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct                                                 >  2:2720750/1‑100 (MQ=255)
                                                 aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct                                                 <  1:1730999/100‑1 (MQ=255)
                                                     gcgACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccc                                             <  1:2424922/100‑1 (MQ=255)
                                                      cgACTATCATTCATTGCAGGAATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccct                                            >  1:2067376/1‑100 (MQ=255)
                                                         cTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccctccc                                         >  1:2505430/1‑100 (MQ=255)
                                                          tATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACtcccctcccc                                        >  2:603716/1‑100 (MQ=255)
                                                              attcattGCAGGCATTAGCACACCGGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                    <  1:2699170/100‑1 (MQ=255)
                                                              attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCgg                                     >  1:1771293/1‑99 (MQ=255)
                                                              attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                    >  2:575170/1‑100 (MQ=255)
                                                                tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGTCCTGATGTTAATGATCTCCTTTCACACGTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                    <  1:1852333/98‑1 (MQ=255)
                                                                tcattGCAGGCATTAGCACACCAGGCCGCATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACCCCCCTCCCCCGGAAc                                  <  1:2640664/100‑1 (MQ=255)
                                                                  attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGGTAATGATACCCTTTTACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACaa                                >  1:2387173/1‑100 (MQ=255)
                                                                   ttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAg                               <  1:1742552/100‑1 (MQ=255)
                                                                     gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                             >  2:193523/1‑100 (MQ=255)
                                                                     gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                             >  1:11151/1‑100 (MQ=255)
                                                                     gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                             >  2:781244/1‑100 (MQ=255)
                                                                      cAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTc                            >  1:2667325/1‑100 (MQ=255)
                                                                            ttAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTc                      <  1:3288179/100‑1 (MQ=255)
                                                                             tAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                     >  1:1174978/1‑100 (MQ=255)
                                                                             attaacatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGCAACAAGCTCGGTTTCa                     >  1:2091151/5‑100 (MQ=25)
                                                                              aggacacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                     <  1:1855192/96‑1 (MQ=255)
                                                                              aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                    >  2:1527683/1‑100 (MQ=255)
                                                                               taacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                   >  2:143044/7‑100 (MQ=37)
                                                                               gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                    >  1:887000/1‑99 (MQ=255)
                                                                               gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                   <  2:2799745/100‑1 (MQ=255)
                                                                               gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                   <  1:603713/100‑1 (MQ=255)
                                                                               gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                   <  1:1454779/100‑1 (MQ=255)
                                                                                 acatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                 >  2:416286/1‑100 (MQ=37)
                                                                                 acatCAGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                 >  2:3012458/1‑100 (MQ=38)
                                                                                  catcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa                >  2:2045939/4‑100 (MQ=37)
                                                                                  catcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCCGAACCAGCTCGGTTTCATGac                 >  1:91977/4‑98 (MQ=25)
                                                                                  cacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa                >  2:2015189/1‑100 (MQ=255)
                                                                                   acCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa               >  2:3339532/1‑100 (MQ=255)
                                                                                    ccAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg              >  2:302109/1‑100 (MQ=255)
                                                                                     cAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGc             >  2:1243817/1‑100 (MQ=255)
                                                                                     cAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGc             <  2:2055713/100‑1 (MQ=255)
                                                                                         ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTcc         <  1:1070181/100‑1 (MQ=255)
                                                                                         ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCACCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTcc         >  1:2520495/1‑100 (MQ=255)
                                                                                          cGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTccc        <  1:1007257/100‑1 (MQ=255)
                                                                                             aTGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt     <  1:440843/100‑1 (MQ=255)
                                                                                               gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg   <  1:2202209/100‑1 (MQ=255)
                                                                                               gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg   <  2:2505449/100‑1 (MQ=255)
                                                                                                tgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa  >  2:2490089/1‑100 (MQ=255)
                                                                                                tgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa  >  2:2398610/1‑100 (MQ=255)
                                                                                                |                                                                                                   
GCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGA  >  NC_000913/2028063‑2028258

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