Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,160 G→C 17.9% 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,1600GC17.9% 275.8 / 55.2 112intergenic (‑143/+28)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (46/46);  new base C (7/13);  total (53/59)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

Only 100 of 120 total aligned reads displayed.

GGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGA  >  NC_000913/2028062‑2028258
                                                                                                  |                                                                                                  
ggCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGa                                                                                                   <  2:3963847/100‑1 (MQ=255)
  ngATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc                                                                                                 >  2:1020120/2‑100 (MQ=255)
  cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc                                                                                                 >  2:692961/1‑100 (MQ=255)
  cGATAGCGTACCGCAAAACAGCGATAATTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                >  2:3884619/1‑100 (MQ=255)
       gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                            <  2:2099253/100‑1 (MQ=255)
       gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                            >  2:638767/1‑100 (MQ=255)
       gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                            <  1:2236168/100‑1 (MQ=255)
        cGTACCGCAAAACCGCGATAATTTTATACTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                           <  2:649023/100‑1 (MQ=255)
        cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                           <  2:1644696/100‑1 (MQ=38)
        cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                           <  1:4398074/100‑1 (MQ=38)
        cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                           <  1:2190504/100‑1 (MQ=255)
        cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                           <  2:691410/100‑1 (MQ=255)
        cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                           <  2:48496/100‑1 (MQ=255)
        cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                           <  1:2676646/100‑1 (MQ=255)
         gTACCGCTAAAAAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg                                                                                          >  1:2301938/1‑100 (MQ=255)
         gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg                                                                                          <  1:2883879/100‑1 (MQ=255)
         gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTg                                                                                          <  2:2033701/100‑1 (MQ=255)
         gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTg                                                                                          <  2:1856450/100‑1 (MQ=255)
          tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg                                                                                          >  2:3602267/1‑99 (MQ=255)
          tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                         <  2:2820577/100‑1 (MQ=255)
          tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                         <  2:3067435/100‑2 (MQ=255)
           aCCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTggt                                                                                        <  2:3834196/100‑3 (MQ=38)
           aCCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                        >  2:1226618/1‑100 (MQ=255)
            ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                       >  1:3535739/1‑100 (MQ=255)
             cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgt                                                                                      <  1:2926269/100‑1 (MQ=38)
              gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGtt                                                                                     >  2:1408098/1‑100 (MQ=255)
              gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGtt                                                                                     >  2:1469136/1‑100 (MQ=255)
               cAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtgc                                                                                    <  2:3409330/100‑3 (MQ=37)
                aaagCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtgc                                                                                    <  1:2988943/99‑3 (MQ=25)
                 aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtgc                                                                                    <  1:2989159/98‑3 (MQ=37)
                 aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                  <  2:1141156/100‑1 (MQ=255)
                  aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                 >  2:1933332/1‑100 (MQ=255)
                  aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                 <  1:2202674/100‑1 (MQ=255)
                     aGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATa                                                                              >  2:2466846/1‑100 (MQ=255)
                     aGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATa                                                                              >  1:2247231/1‑100 (MQ=255)
                      gCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAc                                                                             >  1:3876494/1‑100 (MQ=255)
                       cGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAcc                                                                            >  1:2633473/1‑100 (MQ=255)
                       cGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAcc                                                                            <  2:344895/100‑1 (MQ=255)
                       cGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAcc                                                                            <  1:2936116/100‑1 (MQ=255)
                        gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAccc                                                                           <  2:2169879/100‑1 (MQ=255)
                          tAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCtt                                                                         <  2:343057/100‑1 (MQ=255)
                             ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTca                                                                      <  1:3435953/100‑1 (MQ=255)
                             ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTca                                                                      >  1:1423564/1‑100 (MQ=255)
                             ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTACAGGCATTAGCACACAAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTca                                                                      >  1:3720755/1‑100 (MQ=255)
                                tataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCacac                                                                   <  2:221489/100‑1 (MQ=255)
                                  taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                 >  1:261660/1‑100 (MQ=255)
                                  taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                 <  2:474376/100‑1 (MQ=255)
                                      ggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGGCTGATGTTAATGATACCCTTTCACACttcttc                                                             <  2:3816103/100‑1 (MQ=255)
                                       gtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTcc                                                            >  2:2558212/1‑100 (MQ=255)
                                          gCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc                                                          >  1:3391536/1‑99 (MQ=255)
                                                    aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                               <  1:769265/100‑1 (MQ=255)
                                                     cgcgACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActccc                                              <  1:2743405/100‑1 (MQ=255)
                                                       cggcTATCATTCATTGCAGGCATTAGCTCACCAGGCCGGTTGTGCTCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccn                                             <  1:1963506/96‑2 (MQ=255)
                                                           tATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccctccc                                         <  1:1700645/99‑1 (MQ=255)
                                                             tcattcatGCAGCCATTAGCACACCAGGCCGGATGTGCTCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTccccc                                       <  1:1765022/98‑1 (MQ=255)
                                                              cattcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCgg                                     >  2:2863279/1‑100 (MQ=255)
                                                              cattcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCgg                                     >  2:1425555/1‑100 (MQ=255)
                                                               attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                    >  1:2799394/1‑100 (MQ=255)
                                                               attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                    >  1:1705135/1‑100 (MQ=255)
                                                                ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGaa                                   <  2:865235/100‑1 (MQ=255)
                                                                  cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa                                 >  1:2641234/1‑100 (MQ=255)
                                                                   attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACaa                                <  1:844368/100‑1 (MQ=255)
                                                                      gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                             >  1:3128241/1‑100 (MQ=255)
                                                                       cAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTc                            >  2:2234322/1‑100 (MQ=255)
                                                                        aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg                           <  2:808838/100‑1 (MQ=255)
                                                                        aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg                           >  2:4095552/1‑100 (MQ=255)
                                                                        aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg                           <  1:3522252/100‑1 (MQ=255)
                                                                         ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACTCTTCTTCCGCAGAACGTCTCCCCTCCCCCGGAACAAGCTCgg                          >  1:3213935/1‑100 (MQ=255)
                                                                         ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                          >  1:1159589/1‑100 (MQ=255)
                                                                            aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGttt                       <  2:499914/100‑1 (MQ=255)
                                                                            aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGttt                       <  2:496919/100‑1 (MQ=255)
                                                                             ttAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTc                      <  1:3886375/100‑1 (MQ=255)
                                                                              tggCACACCAGTCCGAATGTGATCCTGCCTGATGTTAATGATACCCTTTCACACGTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGt                         <  1:1474450/94‑1 (MQ=255)
                                                                              tAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                     >  2:2131630/1‑100 (MQ=255)
                                                                              attaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                     >  1:1388873/9‑100 (MQ=25)
                                                                               ttaacatCAGGCCGGATCACATCCGGCCTGGTGTTAATGCTACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                    >  2:1789811/4‑100 (MQ=25)
                                                                               aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                    >  2:1481311/1‑100 (MQ=255)
                                                                               aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                    >  2:188609/1‑100 (MQ=255)
                                                                               aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                    >  1:171274/1‑100 (MQ=255)
                                                                               aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                    <  2:3609688/100‑1 (MQ=255)
                                                                                 cacaCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGa                  >  2:3194845/1‑100 (MQ=255)
                                                                                  acatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTCTCATGaa                 >  1:4027514/1‑100 (MQ=37)
                                                                                  acacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                 <  2:3201923/100‑1 (MQ=255)
                                                                                   catcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGGAGGACGACTGCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa                >  1:3087463/4‑100 (MQ=25)
                                                                                   catcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa                >  1:642046/4‑100 (MQ=38)
                                                                                   catcAGGCCGGATCACATCCGGCCTGATATTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                     >  1:1365402/4‑95 (MQ=38)
                                                                                    cccAGGCCGCCAGTGCTCCGGCCTGAGTTGAATGATGCCCCTTCACACTTCGTCCGCAGAACGACTCCCCCCCCCCGGAACAAGCTCGGTTTCAt                    <  1:3631013/94‑1 (MQ=255)
                                                                                    acCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACCTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa               <  1:3946861/100‑1 (MQ=255)
                                                                                     ccAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg              >  1:4104327/1‑100 (MQ=255)
                                                                                        ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCtt           >  1:3380942/1‑100 (MQ=255)
                                                                                          ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTcc         >  2:2678958/1‑100 (MQ=255)
                                                                                            ggATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCg       >  1:2801964/1‑100 (MQ=255)
                                                                                              aTGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt     <  2:1362776/100‑1 (MQ=255)
                                                                                               tgtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCtt    <  1:4198192/100‑1 (MQ=255)
                                                                                                gtgATCGGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCGCGGTTTCATGAAAAGCTTCCCGCTTg   <  1:3407807/100‑1 (MQ=255)
                                                                                                gtgATCCGGCCTGATGTTAATGATACCCTTTGACACTGCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCGCGGTTTCATGAAAAGCTTCCCGCTTg   <  1:2968679/100‑1 (MQ=255)
                                                                                                gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg   >  2:3502914/1‑100 (MQ=255)
                                                                                                gtgATCCGGCCTGATGTTAATGACACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg   >  1:2917003/1‑100 (MQ=255)
                                                                                                 tgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa  >  1:3299696/1‑100 (MQ=255)
                                                                                                 tgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa  >  1:3953593/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
GGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGA  >  NC_000913/2028062‑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: