Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,160 G→C 15.4% 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,1600GC15.4% 316.3 / 50.7 118intergenic (‑143/+28)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (48/51);  new base C (7/11);  total (56/62)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-01

Only 100 of 127 total aligned reads displayed.

GCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAC  >  NC_000913/2028063‑2028259
                                                                                                 |                                                                                                   
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)
     aGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc                                                                                              >  2:2929009/1‑100 (MQ=255)
      gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                             >  2:3260918/1‑100 (MQ=255)
       cGTACCGTAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                            <  1:1045430/100‑1 (MQ=255)
       cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                            <  1:943584/100‑1 (MQ=255)
        gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg                                                                                           <  2:1755241/100‑1 (MQ=255)
        gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg                                                                                           >  1:332264/1‑100 (MQ=255)
         tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                          <  2:293149/100‑1 (MQ=255)
         tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                          <  2:1093398/100‑2 (MQ=255)
          aCCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                         <  1:2488893/100‑1 (MQ=255)
           ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                        >  1:2462833/1‑100 (MQ=255)
           ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                        >  2:1911576/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)
             gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtg                                                                                      <  2:2277861/100‑2 (MQ=37)
             gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtg                                                                                      <  1:167463/100‑2 (MQ=37)
             gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtg                                                                                      <  2:2277860/100‑2 (MQ=37)
             gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGtt                                                                                      <  2:1058864/100‑1 (MQ=255)
             gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGt                                                                                       >  1:1779269/1‑99 (MQ=255)
                    aGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtgctaatgc                                                                               <  1:1831067/100‑9 (MQ=21)
                    aGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATa                                                                               >  2:2127605/1‑100 (MQ=255)
                     gCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtgctaatgcc                                                                              <  2:2758845/100‑10 (MQ=21)
                     gCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAc                                                                              <  1:2649124/100‑1 (MQ=255)
                       gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATATTAATGATAccc                                                                            >  1:3064904/1‑100 (MQ=255)
                           aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                        <  1:1575369/100‑1 (MQ=255)
                           aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                        <  2:1132938/100‑1 (MQ=255)
                             tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac                                                                      <  2:1649609/100‑1 (MQ=255)
                              ttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcaca                                                                     >  2:2623573/1‑100 (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)
                                      gtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTcc                                                             <  2:3141429/100‑1 (MQ=255)
                                       tggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCg                                                            <  2:572782/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)
                                             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)
                                                 aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct                                                  >  2:2720750/1‑100 (MQ=255)
                                                 aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct                                                  <  1:1730999/100‑1 (MQ=255)
                                                 aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct                                                  <  1:3240172/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)
                                                              attcattGCAGGCATTAGCCCACCAGGCAGGAGGTGATCCGGCCTGATGTTAATGATACCCTTTCAGAATACTTCCGGAGAACGACtgccctcccc                                         >  1:1273590/1‑96 (MQ=255)
                                                              attcattGCAGGCATTAGCACACCGGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                     <  1:2699170/100‑1 (MQ=255)
                                                              attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCgg                                      >  1:1771293/1‑99 (MQ=255)
                                                                tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGTCCTGATGTTAATGATCTCCTTTCACACGTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                     <  1:1852333/98‑1 (MQ=255)
                                                                tcattGCAGGCATTAGCACACCAGGCCGCATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACCCCCCTCCCCCGGAAc                                   <  1:2640664/100‑1 (MQ=255)
                                                                 cattGCAGGCATTAGCACACCAGGCCGGATGTGATACGGCCTGATGTTAATGATACCCATTCACACTTATTCCGCAGAACGACTCCACTCCCCCGGAACa                                  >  1:3374403/1‑100 (MQ=255)
                                                                   ttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAg                                <  1:1742552/100‑1 (MQ=255)
                                                                     gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                              >  2:781244/1‑100 (MQ=255)
                                                                     gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                              >  2:193523/1‑100 (MQ=255)
                                                                     gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                              >  1:11151/1‑100 (MQ=255)
                                                                      cAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTc                             >  1:2667325/1‑100 (MQ=255)
                                                                        ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                           >  2:194233/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                    <  1:1454779/100‑1 (MQ=255)
                                                                               gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                    >  2:236684/1‑100 (MQ=255)
                                                                               gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                    <  1:603713/100‑1 (MQ=255)
                                                                                 acatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                  >  2:416286/1‑100 (MQ=37)
                                                                                  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)
                                                                                    tcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg               >  1:224144/2‑100 (MQ=38)
                                                                                    ccAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg               >  2:302109/1‑100 (MQ=255)
                                                                                     cAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGc              <  2:2055713/100‑1 (MQ=255)
                                                                                       ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCtt            >  2:396673/1‑100 (MQ=255)
                                                                                       ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGCAACAAGCTCGGTTTCATGAAAAGCt             >  2:1364630/1‑99 (MQ=255)
                                                                                         ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTcc          <  1:1105873/100‑1 (MQ=255)
                                                                                         ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTcc          <  1:1070181/100‑1 (MQ=255)
                                                                                          cGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTccc         <  1:1007257/100‑1 (MQ=255)
                                                                                             aTGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt      <  1:440843/100‑1 (MQ=255)
                                                                                                tgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa   >  2:2490089/1‑100 (MQ=255)
                                                                                                tgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa   >  2:2398610/1‑100 (MQ=255)
                                                                                                 gATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAc  >  1:1405373/1‑100 (MQ=255)
                                                                                                 gATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAc  >  2:2544707/1‑100 (MQ=255)
                                                                                                 gATCCAGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCACCTCCCCCGGAACAAGCTCGGTTTAATGAAAAGCTTCCCGCTTAAc  >  1:608869/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAC  >  NC_000913/2028063‑2028259

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: