Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,057,140 A→G 20.5% intergenic (‑191/+38) ibsC ← / ← serA toxic membrane protein/D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,057,1400AG20.5% 395.4 / 105.1 207intergenic (‑191/+38)ibsC/serAtoxic membrane protein/D‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (84/80);  new base G (27/15);  total (111/96)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

Only 100 of 210 total aligned reads displayed.

CAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTT  >  NC_000913/3057041‑3057238
                                                                                                   |                                                                                                  
 aaCGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAg                                                                                                   >  1:1326734/1‑100 (MQ=255)
   cGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTc                                                                                                 >  2:3265941/1‑100 (MQ=255)
    gCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCa                                                                                                <  2:2389102/100‑1 (MQ=255)
     cATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAg                                                                                               >  2:406178/1‑100 (MQ=255)
     cATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAg                                                                                               <  2:2740851/100‑1 (MQ=255)
     cATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAg                                                                                               >  2:1638187/1‑100 (MQ=255)
      aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGt                                                                                              <  2:171540/100‑1 (MQ=255)
       ttGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                                                                                             <  2:360581/100‑1 (MQ=255)
         gATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAc                                                                                           >  2:1874517/1‑100 (MQ=255)
          aTCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAcc                                                                                          <  1:3468148/100‑1 (MQ=255)
            cTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTg                                                                                        <  2:3350177/100‑1 (MQ=255)
             tGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGc                                                                                       >  2:1245826/1‑100 (MQ=255)
                cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCg                                                                                    <  1:1863307/100‑1 (MQ=255)
                cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCg                                                                                    <  1:735322/100‑1 (MQ=255)
                cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                    >  1:1178588/1‑100 (MQ=255)
                 tttGAATTATTTTCTGGAGCAGAATCGCAAAGTAGAATGAGCAAAGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGACAGTGACATGCCCGt                                                                                   >  1:2169406/1‑100 (MQ=255)
                   tGATTTATTTTCGGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACCGTGTGGAGAATGGATAAAAAAACGGGCAGGTCTCTGACTTGCCCGttt                                                                                 <  1:3751374/100‑2 (MQ=255)
                    gATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttt                                                                                <  1:3422171/100‑3 (MQ=38)
                    gATTTATTTTCTGGAGCAGACTCGAAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttt                                                                                 <  1:1518694/99‑2 (MQ=38)
                    gATGTATTTTCTGGAGCAGACTCGAAAAGAAGAAGTAACAACGCGCCAACAGTGTGGAGAAGGGATTAAAAAACTGGCATGTCACTGACTTGCCCGtttt                                                                                <  1:2910108/100‑3 (MQ=255)
                     atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGACGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGa                                                                                <  1:797020/99‑1 (MQ=255)
                     atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttttt                                                                               <  2:629037/100‑4 (MQ=38)
                     atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttttt                                                                               <  1:1722015/100‑4 (MQ=38)
                      tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAGCGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt                                                                              >  1:2443078/1‑100 (MQ=255)
                       ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAttt                                                                             <  2:1396960/100‑1 (MQ=255)
                       ttattgTCTGGAGTAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGCGAAGGGATAAAAAGACGGGCAAGTCAGTGACCTGCCCGTTgctct                                                                             >  1:2727847/1‑96 (MQ=255)
                        tatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttttttta                                                                            <  1:508890/100‑7 (MQ=37)
                        tatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtttt                                                                            >  2:2483928/1‑100 (MQ=255)
                            ttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCaga                                                                        <  1:1619814/100‑1 (MQ=255)
                            ttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCaga                                                                        <  2:2028489/100‑1 (MQ=255)
                              cTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAgaga                                                                      >  1:165262/1‑100 (MQ=255)
                              cTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGAGATAAAAGAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAgaga                                                                      >  1:1311273/1‑100 (MQ=255)
                                ggAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                    >  2:1003648/1‑100 (MQ=255)
                                ggAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                    <  1:2290462/100‑1 (MQ=255)
                                   gCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                 >  1:2217944/1‑100 (MQ=255)
                                      gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAt                                                              >  2:1791279/1‑100 (MQ=255)
                                       tcTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAtt                                                             <  2:3195521/99‑1 (MQ=255)
                                        cTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTa                                                            >  1:3221867/1‑100 (MQ=255)
                                          cGCAACGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                          >  2:1073975/1‑100 (MQ=255)
                                          cGCAAAGTAGAATGCGCAACGCGGCAACTGTGTGGAGAAGGGATAAAAATACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAg                                                           <  1:416137/99‑1 (MQ=255)
                                          cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAga                                                          >  1:3195504/1‑99 (MQ=255)
                                          cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                          >  1:2487403/1‑100 (MQ=255)
                                          cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                          >  2:1085685/1‑100 (MQ=255)
                                          cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                          >  2:2691738/1‑100 (MQ=255)
                                              aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                      >  2:3253560/1‑100 (MQ=255)
                                              aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                      <  2:3753082/100‑1 (MQ=255)
                                              aaGTAGAATACGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                      <  1:1058146/100‑1 (MQ=255)
                                               aGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagc                                                     <  1:4124238/100‑1 (MQ=255)
                                               aGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagc                                                     >  2:2081665/1‑100 (MQ=255)
                                                gTAGAATGCACAACGCTGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCCAGTCAGTGACCTGCCCGTTGATTATCAGAGAAGGGGAATTAGTAc                                                        >  1:4043891/1‑96 (MQ=255)
                                                   gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGAc                                                 >  1:1000867/1‑100 (MQ=255)
                                                   gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGAc                                                 <  2:1946873/100‑1 (MQ=255)
                                                   gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGAc                                                 <  1:3259433/100‑1 (MQ=255)
                                                      tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACggg                                              <  2:683699/100‑1 (MQ=255)
                                                       gcgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgc                                             <  2:4397313/100‑1 (MQ=255)
                                                       gcgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgc                                             >  2:3561413/1‑100 (MQ=255)
                                                            aCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGaa                                        >  1:1939398/1‑100 (MQ=255)
                                                             cgcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGaa                                        <  1:4247090/99‑1 (MQ=255)
                                                              gcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATg                                      <  2:2336302/100‑1 (MQ=255)
                                                               cgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATgg                                     <  1:1702527/100‑1 (MQ=255)
                                                                  cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAc                                  <  2:154966/100‑1 (MQ=255)
                                                                  cAACGGCGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAc                                  >  1:2697912/1‑100 (MQ=255)
                                                                     cGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAccc                                <  2:2813534/99‑1 (MQ=255)
                                                                      ggTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCgg                              <  2:3662791/100‑1 (MQ=255)
                                                                         gtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAAt                           <  1:2548301/100‑1 (MQ=255)
                                                                          tgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATa                          <  1:602839/100‑1 (MQ=255)
                                                                           ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAg                         >  2:1825613/1‑100 (MQ=255)
                                                                           ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAg                         <  1:3980659/100‑1 (MQ=255)
                                                                             agaAGGGATAAAAAAACGGGCAAGTCGGTGACCTGCCCGGTGATTCTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAAt                           <  1:2806315/96‑1 (MQ=255)
                                                                             agaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCt                       >  1:2086297/1‑100 (MQ=255)
                                                                                aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTc                    <  2:291503/100‑1 (MQ=255)
                                                                                   tgaaaatcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAtt                 >  2:3118894/9‑100 (MQ=37)
                                                                                   tgaaaatcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAtt                 >  2:96606/9‑100 (MQ=37)
                                                                                   gATAAGAAAACGGGCAAGTCAGTGCCCCTCCCGTTGATTTTCAGGGAAGGGTAATTAGTACAGCAGAGGGGCGCCAATGGTATCGGGTAAAGCCTTCAt                  >  1:2655106/1‑99 (MQ=255)
                                                                                     tAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc               >  1:81391/1‑100 (MQ=255)
                                                                                      aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGcc              <  2:2500838/100‑1 (MQ=255)
                                                                                       aaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCt             >  1:1130950/1‑100 (MQ=255)
                                                                                       aaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCt             >  1:3295713/1‑100 (MQ=255)
                                                                                         aaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc           <  2:1221565/100‑1 (MQ=255)
                                                                                         aaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc           <  2:2850045/100‑1 (MQ=255)
                                                                                         aaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGCGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc           >  1:461905/1‑100 (MQ=255)
                                                                                          caaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa          >  1:3791075/2‑100 (MQ=39)
                                                                                           aaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAg         >  2:890297/1‑100 (MQ=39)
                                                                                           aaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAg         >  2:1459589/1‑100 (MQ=39)
                                                                                           aaCGGGCAGGTCAACGACTTGCCCGTTGATTTTCAGCGCAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACACGGAATAGCCTTCATTGCCTGCa          >  1:2381233/1‑99 (MQ=255)
                                                                                            actggCAATGCAGAGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTgacgc        >  1:4116064/4‑96 (MQ=255)
                                                                                            aCGGGCAGGTCACTGTCTTGCCCGTTGATTTTCAGAGACGGGGAATTAGTACAGCAGACGGGCGCGGATGGTACGCGGAATAGCTTTCATTGcc              >  1:1226022/1‑94 (MQ=255)
                                                                                            aCGGGCAGGTCACTGACTTGCCTGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAGCCGGAATAGCTTTCATTGCCTGGAg         >  1:3075076/1‑99 (MQ=255)
                                                                                            aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc        >  2:2226646/1‑100 (MQ=39)
                                                                                            aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc        >  2:3136832/1‑100 (MQ=255)
                                                                                            aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc        >  2:768155/1‑100 (MQ=39)
                                                                                            aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc        >  2:1704840/1‑100 (MQ=255)
                                                                                             cGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcg       >  2:1721899/1‑100 (MQ=255)
                                                                                              gggCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcgc      >  2:798670/1‑100 (MQ=255)
                                                                                              gggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcgc      <  2:342288/100‑1 (MQ=255)
                                                                                               ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGAAATAGCTTTCATTGCCTGCAgcgc      >  1:3283562/1‑99 (MQ=255)
                                                                                                gCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACTCGGAATAGCTTTCATCGCCTTCTACGCtt    >  1:1390546/1‑100 (MQ=255)
                                                                                                gCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtt    <  2:3872178/100‑1 (MQ=255)
                                                                                                 cgaGTCAGTGCCCTCCGCGTTGATGTTTAGAGAGAGGGAATTAGTACAGCAGACGGGCGGGAATGGTACCCGGATTAGCTGTCATTGCCTGCAGCGCttt   <  1:2048135/98‑1 (MQ=255)
                                                                                                 cAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttt   >  1:301413/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTT  >  NC_000913/3057041‑3057238

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