Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 207,586 T→C 6.2% Y821H (TAC→CAC)  dnaE → DNA polymerase III alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913207,5860TC 6.2% 193.1 / 3.4 64Y821H (TAC→CAC) dnaEDNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base C (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

TTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAGGCGGCTACTTCCGCGAACTGTTTGATCTCTGCGCCCGTACC  >  NC_000913/207439‑207732
                                                                                                                                                   |                                                                                                                                                  
ttATGGCATAAAGAGCTCTATCCTGCCGAGTTTATGCCGGCTGTAATGACCGCCGATATTGACAACACCGAGAAGGTGGTGGGGCTGTTGGATGAGTGCTGGCGCATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCttt                                                                                                                                                    <  2:204086/148‑1 (MQ=255)
      cTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCAtt                                                                                                                                              <  2:70716/148‑1 (MQ=255)
        gAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTc                                                                                                                                            >  1:202609/1‑148 (MQ=255)
        gAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTc                                                                                                                                            <  4:251691/148‑1 (MQ=255)
         aaaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTc                                                                                                                                            <  5:178615/147‑1 (MQ=255)
            gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTCACCATTTCCACg                                                                                                                                        >  3:73590/1‑148 (MQ=255)
            gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTCACCATTTCCACg                                                                                                                                        >  3:68453/1‑148 (MQ=255)
                aCTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCaa                                                                                                                                    >  2:366623/1‑148 (MQ=255)
                 cTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAac                                                                                                                                   >  7:415685/1‑148 (MQ=255)
                  tATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTcc                                                                                                                                           >  4:340596/1‑139 (MQ=255)
                  tATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTcc                                                                                                                                           <  3:340596/139‑1 (MQ=255)
                     ccTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAacg                                                                                                                                  <  4:253916/145‑1 (MQ=255)
                     ccTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAacg                                                                                                                                  >  3:253916/1‑145 (MQ=255)
                         cGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAacgacg                                                                                                                               <  8:82833/144‑1 (MQ=255)
                         cGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAacgacg                                                                                                                               >  7:82833/1‑144 (MQ=255)
                               tttaggcggcGGTTATGACCGCCGATACGGACAACACCGAGCAGTTGGTGGGTCTGGTGGATGACTGCTGGCGGCTGGGGCTGAAAAACCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAAt                                                                                                                     <  2:136024/144‑1 (MQ=255)
                                tATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAAt                                                                                                                     <  8:266326/147‑1 (MQ=255)
                                 aTGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGTATGGGGCTGAAAATCCTGCCACCAGATATAAACTCAGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCg                                                                                                                   >  4:169722/1‑148 (MQ=255)
                                  tGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCg                                                                                                                   <  1:366623/147‑1 (MQ=255)
                                    gcggcgGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgtgt                                                                                                                <  8:415685/148‑1 (MQ=255)
                                    gcggcgGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgtgt                                                                                                                <  7:224688/148‑1 (MQ=255)
                                     cggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGGTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCa                                                                                                                                          <  2:98414/121‑1 (MQ=255)
                                     cggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGGTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCa                                                                                                                                          >  1:98414/1‑121 (MQ=255)
                                             aTGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTg                                                                                                       >  8:32236/1‑148 (MQ=255)
                                               gACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTcc                                                                                                                                           >  8:257061/1‑110 (MQ=255)
                                               gACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTcc                                                                                                                                           <  7:257061/110‑1 (MQ=255)
                                               gACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGgc                                                                                                     >  5:365248/1‑148 (MQ=255)
                                               gACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGgc                                                                                                     >  7:417484/1‑148 (MQ=255)
                                                      gATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGt                                                                                                                                       <  8:171920/107‑1 (MQ=255)
                                                      gATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGt                                                                                                                                       >  7:171920/1‑107 (MQ=255)
                                                         atGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAg                                                                                           >  2:133027/1‑148 (MQ=255)
                                                           ggACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgtgt                                                                                                                <  7:121227/125‑1 (MQ=255)
                                                           ggACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgtgt                                                                                                                >  8:121227/1‑125 (MQ=255)
                                                              caacaCCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTc                                                                                      <  6:351954/148‑1 (MQ=255)
                                                                 cacCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGt                                                                                   >  3:93177/1‑148 (MQ=255)
                                                                  acCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgtg                                                                                                                 <  4:304015/117‑1 (MQ=255)
                                                                  acCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgtg                                                                                                                 >  3:304015/1‑117 (MQ=255)
                                                                    cGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGaa                                                                                <  2:202609/148‑1 (MQ=255)
                                                                           gtggtgGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTCACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGa                                                                         <  4:68453/148‑1 (MQ=255)
                                                                           gtggtgGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTCACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGa                                                                         <  4:73590/148‑1 (MQ=255)
                                                                            tggtggGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGAt                                                                        <  6:365248/148‑1 (MQ=255)
                                                                             ggtggGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGAtt                                                                       >  4:138112/1‑148 (MQ=255)
                                                                               tggGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGAtt                                                                       >  8:149971/1‑146 (MQ=255)
                                                                                      ggtggATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCACCGACGACGGCGAAATCGTGTATGGTTTTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCcat                                                               >  3:333409/1‑147 (MQ=255)
                                                                                       gtggATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCcatc                                                              <  7:206425/147‑1 (MQ=255)
                                                                                        tggATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCcatcat                                                            >  1:163243/1‑148 (MQ=255)
                                                                                         ggATGAGTGCTGGCGGATGTGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCatcatc                                                           <  2:269772/148‑1 (MQ=255)
                                                                                           aTGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGAGCAAAGGGGGCGGTGAAGGTCCGATTGAGGCCATCATCGa                                                         <  8:383273/148‑1 (MQ=255)
                                                                                                tGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGccc                                                    >  6:156735/1‑148 (MQ=255)
                                                                                                    ggCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGtaa                                                >  2:162243/1‑148 (MQ=255)
                                                                                                     gCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGtaat                                               <  3:138112/148‑1 (MQ=255)
                                                                                                       ggATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGtaataa                                             <  5:156735/148‑1 (MQ=255)
                                                                                                       ggATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGtaataa                                             <  4:93177/148‑1 (MQ=255)
                                                                                                       ggATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGtaataa                                             <  1:162243/148‑1 (MQ=255)
                                                                                                          tGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAgg                                          >  5:290580/1‑148 (MQ=255)
                                                                                                            gggCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAggcg                                        <  3:156189/148‑1 (MQ=255)
                                                                                                               cTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgtgt                                                                                                                <  1:36391/73‑1 (MQ=255)
                                                                                                               cTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgtgt                                                                                                                >  2:36391/1‑73 (MQ=255)
                                                                                                                   aaaTCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGTCATCATCGAAGCCCGTAATAAAGGCGGCTACtt                                 >  6:132808/1‑148 (MQ=255)
                                                                                                                     aTCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGAGTGAGCCCATCATCGAAGCCCGTAATAAAGGCGGCTACTTcc                               <  1:88503/148‑1 (MQ=255)
                                                                                                                              ccaGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAGGCGGCTACTTCCGCGAACTGt                      >  8:404353/1‑148 (MQ=255)
                                                                                                                              ccaGATATAAACTCCGGTCTATAACATATCCACGTCAACGACGACGGCTAAAACGTGTAAGGTATTGTCGCGATCAAACGGGTAGCTGAAGGTCCGATTGAGGCCATCAACGAAGCCCGTAATAAAGGCGGCTACTAACGCGAACTGt                      >  1:349212/1‑148 (MQ=255)
                                                                                                                                aGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAGGCGGCTACTTCCGCGAACTGttt                    <  2:163243/148‑1 (MQ=255)
                                                                                                                                     taAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAGGCGGCTACTTCCGCGAACTGTTTGAtct               <  6:290580/148‑1 (MQ=255)
                                                                                                                                              gTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAGGCGGCTACTTCCGCGAACTGTTTGATCTCTGCGCCCg      >  2:223454/1‑148 (MQ=255)
                                                                                                                                                  ttACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAGGCGGCTACTTCCGCGAACTGTTTGATCTCTGCGCCCGTAcc  >  4:119057/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAGGCGGCTACTTCCGCGAACTGTTTGATCTCTGCGCCCGTACC  >  NC_000913/207439‑207732

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑