Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,159 T→A 15.2% 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.2% 281.4 / 41.7 113intergenic (‑142/+29)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (48/47);  new base A (4/13);  total (52/61)
Fisher's exact test for biased strand distribution p-value = 6.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

Only 100 of 116 total aligned reads displayed.

GGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGA  >  NC_000913/2028062‑2028258
                                                                                                 |                                                                                                   
ggCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGa                                                                                                   >  1:2952296/1‑100 (MQ=255)
 gCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                  <  2:3398611/100‑1 (MQ=255)
   gATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                <  1:3796385/100‑1 (MQ=255)
      aGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc                                                                                             >  1:692868/1‑100 (MQ=255)
       gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGcc                                                                                            <  2:2233353/100‑1 (MQ=255)
       gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                            <  1:3533882/100‑1 (MQ=255)
        cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                           <  1:429333/100‑1 (MQ=38)
        cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                           <  1:1917978/100‑1 (MQ=38)
        cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                           <  1:1472465/100‑1 (MQ=255)
        cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTCCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                           <  2:2039966/100‑1 (MQ=255)
         gTACCGCAAAACAGGGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                           <  2:3166777/99‑1 (MQ=255)
         gTACCGCAAAACAGCGATAATTTTATAAGGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTg                                                                                          <  1:185266/100‑1 (MQ=255)
          tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATTTGATCCGGCCtta                                                                                         <  1:2589568/100‑3 (MQ=255)
          tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                         <  1:2034157/100‑2 (MQ=255)
          tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                         <  1:3384074/100‑2 (MQ=255)
          tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                         <  1:2855549/100‑2 (MQ=255)
           aCCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                         <  2:1985746/99‑2 (MQ=39)
            ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                       >  2:3558486/1‑100 (MQ=255)
            ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                       <  1:1102924/100‑1 (MQ=255)
             cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgt                                                                                      <  2:3570850/100‑1 (MQ=37)
             cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGt                                                                                      >  2:1940351/1‑100 (MQ=255)
              gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGtt                                                                                     >  1:3619601/1‑100 (MQ=255)
              gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGtt                                                                                     <  1:1368727/100‑1 (MQ=255)
                aaaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa                                                                                   <  2:1473888/100‑1 (MQ=255)
                aaaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa                                                                                   <  2:2166915/100‑1 (MQ=255)
                aaaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa                                                                                   <  2:2121232/100‑1 (MQ=255)
                 aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                  <  2:2402554/100‑1 (MQ=255)
                   aCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtgctaat                                                                                <  1:3350014/100‑7 (MQ=25)
                   aCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGa                                                                                <  2:767222/100‑1 (MQ=255)
                    cAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGAt                                                                               <  2:2244774/100‑1 (MQ=255)
                      gCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAc                                                                             <  2:159697/100‑1 (MQ=255)
                       cGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAcc                                                                            >  1:2873580/1‑100 (MQ=255)
                         aTAATTTTATAATGGTGGCATGATAAGCCGCGACAATCATTCATTGCAGGCATCAGAACACAAGTCCGGATGTGATCCGGCCACATGTTAAtacca                                                                              >  2:2229827/1‑92 (MQ=255)
                          tAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCtt                                                                         >  1:216294/1‑100 (MQ=255)
                          tAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCtt                                                                         <  2:3192355/100‑1 (MQ=255)
                            aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                       >  2:1635288/1‑100 (MQ=255)
                            aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                       <  2:1995702/100‑1 (MQ=255)
                             ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTca                                                                      >  1:3056141/1‑100 (MQ=255)
                             ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                       <  2:1482954/99‑1 (MQ=255)
                              tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac                                                                     <  1:1001265/100‑1 (MQ=255)
                                 ataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct                                                                  <  2:2964276/100‑1 (MQ=255)
                                 ataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct                                                                  <  1:1116662/100‑1 (MQ=255)
                                    aTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActtct                                                               >  1:1005291/1‑100 (MQ=255)
                                     tggtggCATGATTAGACGCGACTATTATTCCTTCCAGGGAGTAGCACACCAGGCCTGATGTGATCCGGCCTGATGTTAATGATACCGTTGCACActtc                                                                <  2:1543467/98‑1 (MQ=255)
                                     tggtggCATGATAAGACTCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTCATCCGGCCTGATGTTAATGATACCCTTTCACActtctt                                                              >  2:1043147/1‑100 (MQ=255)
                                      ggtggCATGATAAGACGCGACTATAATTCATTGCCGGCATTCGCACACCACGCAGGATGTGATCCGTCCTGATGTTAATGATCCCCTTCAACACtacttc                                                             >  2:3369734/1‑100 (MQ=255)
                                       gtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTcc                                                            <  2:534600/100‑1 (MQ=255)
                                       gtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTcc                                                            <  1:2733552/100‑1 (MQ=255)
                                         ggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc                                                          <  1:2553031/100‑1 (MQ=255)
                                           cATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAg                                                        >  2:2478661/1‑100 (MQ=255)
                                           cATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAg                                                        <  1:3213831/100‑1 (MQ=255)
                                            aTGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGa                                                       >  2:750585/1‑100 (MQ=255)
                                              gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                     <  2:2952275/100‑1 (MQ=255)
                                              gATAAGACGCGACGATAATTCAGTGCAGGCATTAGCACACCACGCCGGATGCGATTCGCCCTGATGATACTGATACCCTTTCGCGCTACTTCCGCAGagc                                                     <  2:2482104/100‑3 (MQ=255)
                                               aTAAGACGCGACTATCATTCATTGCAGGAATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACg                                                    <  2:3044182/100‑1 (MQ=255)
                                                  aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct                                                 >  1:3763163/1‑100 (MQ=255)
                                                  aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCAGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct                                                 >  2:3058044/1‑100 (MQ=255)
                                                   gACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActc                                                >  1:2690312/1‑100 (MQ=255)
                                                      gcgACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccc                                             <  1:1122050/100‑1 (MQ=255)
                                                       cgACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccct                                            >  1:1156128/1‑100 (MQ=255)
                                                         aCTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccctcc                                          <  2:1693957/100‑1 (MQ=255)
                                                         aCTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccctcc                                          <  1:1218284/100‑1 (MQ=255)
                                                          cTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccctccc                                         >  2:1129975/1‑100 (MQ=255)
                                                             tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCg                                      <  2:3706877/100‑1 (MQ=255)
                                                              cattcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCgg                                     <  2:2927949/100‑1 (MQ=255)
                                                               attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                    >  1:3175281/1‑100 (MQ=255)
                                                                ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGaa                                   <  1:3680381/100‑1 (MQ=255)
                                                                ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGaa                                   <  2:3040022/100‑1 (MQ=255)
                                                                 tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc                                  <  2:1562780/100‑1 (MQ=255)
                                                                 tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc                                  >  1:3469916/1‑100 (MQ=255)
                                                                  cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa                                 >  1:2668514/1‑100 (MQ=255)
                                                                   attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACaa                                >  2:2405229/1‑100 (MQ=255)
                                                                    ttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCCCACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAg                               >  2:649382/1‑100 (MQ=255)
                                                                     tGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGc                              >  2:2188365/1‑100 (MQ=255)
                                                                      gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGc                              >  2:3707706/1‑99 (MQ=255)
                                                                      gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                             >  2:3707244/1‑100 (MQ=255)
                                                                      gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                             >  1:1525075/1‑100 (MQ=255)
                                                                      gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                             >  1:1021268/1‑100 (MQ=255)
                                                                      gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCCTACGCAGGACAActacac                                             >  2:3412371/1‑80 (MQ=255)
                                                                       cAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTc                            >  2:399741/1‑100 (MQ=255)
                                                                         ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                          <  1:1492956/100‑1 (MQ=255)
                                                                             ttAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTc                      <  2:1422230/100‑1 (MQ=255)
                                                                              attaacatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTc                      >  1:2039635/5‑99 (MQ=37)
                                                                               aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTTGCCGCCCCTGGTACAAGCTCGGTTTctt                    >  1:1127764/1‑98 (MQ=255)
                                                                               aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTc                      >  2:469974/1‑98 (MQ=255)
                                                                               aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                    >  1:2650769/1‑100 (MQ=255)
                                                                               aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                    <  2:1598621/100‑1 (MQ=255)
                                                                                gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                   >  1:3326840/1‑100 (MQ=255)
                                                                                 cacaCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                    <  2:3543829/98‑1 (MQ=255)
                                                                                 aacatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGa                  >  1:2949636/2‑100 (MQ=37)
                                                                                   cacCAGGCCGGATGTGATCCGGCCTGATGTGAACGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGAGTTCATGaaa                >  2:2508725/1‑100 (MQ=255)
                                                                                     ccAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg              >  1:3515047/1‑100 (MQ=255)
                                                                                        ggCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCtt           >  1:2023610/1‑100 (MQ=255)
                                                                                         gCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTc          >  2:3015144/1‑100 (MQ=255)
                                                                                         gCCGGATGTGATCAGGCCTGATGTTAATGATACCGTTTCACACTTCTTCCGCAGAAGGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTc          >  2:2631765/1‑100 (MQ=255)
                                                                                           cGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTccc        >  2:138974/1‑100 (MQ=255)
                                                                                           cGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTccc        <  1:3770543/100‑1 (MQ=255)
                                                                                              aTGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt     <  1:699566/100‑1 (MQ=255)
                                                                                                gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg   <  1:3585232/100‑1 (MQ=255)
                                                                                                 tgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa  >  2:3159865/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: