Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,810,873 A→G 11.7% H35R (CAC→CGC) ‡ mprA → DNA‑binding transcriptional repressor MprA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,810,8730AG11.7% 200.3 / 16.1 77H35R (CAC→CGC) ‡mprADNA‑binding transcriptional repressor MprA
Reads supporting (aligned to +/- strand):  ref base A (37/31);  new base G (4/5);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCAT  >  NC_000913/2810727‑2811019
                                                                                                                                                  |                                                                                                                                                  
cTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacat                                                                                                                                                 >  2:1193685/1‑150 (MQ=255)
  gTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatgc                                                                                                                                               <  1:1608750/150‑1 (MQ=255)
        aCTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaagac                                                                                                                                         <  1:186776/150‑3 (MQ=255)
        aCTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaagac                                                                                                                                         <  2:1315224/150‑3 (MQ=255)
        aCTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagca                                                                                                                                         >  1:1807447/1‑150 (MQ=255)
        aCTATATAGGCTGAAATTAATGATGTCATACCCATATGGATAGTTCGTTTACTCCGATTTAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTGTTTGCTTGTGCATGCAaagac                                                                                                                                         <  2:1974766/150‑3 (MQ=255)
         cTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagcaa                                                                                                                                        <  1:982478/150‑1 (MQ=255)
           atatCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaagacgag                                                                                                                                      <  2:838362/150‑6 (MQ=255)
            tatCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAgct                                                                                                                                     >  1:249037/1‑150 (MQ=255)
               cgacTGAAAATAATGAGCTCATACCCAAATGGATAGTTCGTTTACGCCCAAAGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGCCTAGTCTTTGCATGCACATGCAAAGCAAgctgct                                                                                                                                  >  2:1864052/4‑150 (MQ=255)
                ggCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGctgctg                                                                                                                                 >  2:165425/1‑150 (MQ=255)
                ggCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGctgctg                                                                                                                                 <  1:871547/150‑1 (MQ=255)
                 gCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTgg                                                                                                                                <  2:47123/150‑1 (MQ=255)
                  cTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGga                                                                                                                               >  1:1395376/1‑150 (MQ=255)
                  cTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAACCAAGCTGCTGga                                                                                                                               >  2:840078/1‑150 (MQ=255)
                     aaaTTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGaa                                                                                                                            >  1:1957106/1‑150 (MQ=255)
                         tAATGAGGTCTTACCCAAATGGATAGTTCGTTTACGCCGATTGAAGAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGTGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGc                                                                                                                        <  2:645099/150‑1 (MQ=255)
                         tAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGc                                                                                                                        <  2:1513548/150‑1 (MQ=255)
                           aTGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCaa                                                                                                                      >  2:1014255/1‑150 (MQ=255)
                           aTGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCaa                                                                                                                      <  1:1821103/150‑1 (MQ=255)
                           aTGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCaa                                                                                                                      >  2:1021144/1‑150 (MQ=255)
                             gAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATa                                                                                                                    >  2:1804105/1‑150 (MQ=255)
                                    tACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCt                                                                                                             >  2:1429123/1‑150 (MQ=255)
                                     aCCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaagac                                                                                                                                         >  2:1905639/1‑119 (MQ=255)
                                     aCCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaagac                                                                                                                                         <  1:1905639/121‑3 (MQ=255)
                                      cccAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGa                                                                                                           >  2:1564445/1‑150 (MQ=255)
                                        cAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTtgcatgca                                                                                                                                                    >  2:1034717/1‑107 (MQ=255)
                                        cAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTtgcatgca                                                                                                                                                    <  1:1034717/107‑1 (MQ=255)
                                        cAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCTTGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAg                                                                                                         <  2:1957106/150‑1 (MQ=255)
                                            tGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGGTTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTc                                                                                                     >  1:760893/1‑150 (MQ=255)
                                            tGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTc                                                                                                     <  2:518396/150‑1 (MQ=255)
                                               aTAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAgg                                                                                                  >  1:949039/1‑150 (MQ=255)
                                                 aGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAgggg                                                                                                >  2:122971/1‑150 (MQ=255)
                                                        ttACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACg                                                                                         <  2:1383354/150‑1 (MQ=255)
                                                         tACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACga                                                                                        <  1:1180266/150‑1 (MQ=255)
                                                          aCGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACgag                                                                                       <  1:1884671/150‑1 (MQ=255)
                                                             cccATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACg                                                                                    <  2:1807447/150‑1 (MQ=255)
                                                             cccATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACg                                                                                    >  1:109390/1‑150 (MQ=255)
                                                             cccATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACg                                                                                    <  1:1021144/150‑1 (MQ=255)
                                                             cccATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACg                                                                                    <  1:1014255/150‑1 (MQ=255)
                                                             cccATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCGGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACg                                                                                    >  1:1451926/1‑150 (MQ=255)
                                                              ccATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCTGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgt                                                                                   >  2:975997/1‑150 (MQ=255)
                                                               cATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgtt                                                                                  >  1:1520529/1‑150 (MQ=255)
                                                                 ttGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgt                                                                                <  2:249037/150‑1 (MQ=255)
                                                                   gAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGttt                                                                              >  2:1002244/1‑150 (MQ=255)
                                                                   gAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGttt                                                                              <  1:1595948/150‑1 (MQ=255)
                                                                    aaCAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTa                                                                             >  1:954518/1‑150 (MQ=255)
                                                                    aaCAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTa                                                                             <  2:1699400/150‑1 (MQ=255)
                                                                       aaaTGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATgg                                                                          >  1:1649934/1‑150 (MQ=255)
                                                                        aaTGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGc                                                                         <  1:1278496/150‑1 (MQ=255)
                                                                           gCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGtt                                                                      <  1:1193685/150‑1 (MQ=255)
                                                                            cTAAAATTTCGCGCCAGCCGCCACGAAGGTTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTg                                                                     <  2:760893/150‑1 (MQ=255)
                                                                              aaaaTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                                                                   >  1:456892/1‑150 (MQ=255)
                                                                               aaaTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAtt                                                                  >  1:1531884/1‑150 (MQ=255)
                                                                               aaaTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAtt                                                                  <  1:1564445/150‑1 (MQ=255)
                                                                                 aTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTAc                                                                >  1:435498/1‑150 (MQ=255)
                                                                                       cgcCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa                                                          <  2:648075/150‑1 (MQ=255)
                                                                                        gcCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAg                                                         >  2:451866/1‑150 (MQ=255)
                                                                                        gcCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAg                                                         <  1:1804105/150‑1 (MQ=255)
                                                                                            gccgccACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtctc                                                     <  2:954518/150‑1 (MQ=255)
                                                                                               gccACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg                                                  >  1:649547/1‑150 (MQ=255)
                                                                                                 cACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaa                                                >  2:1189906/1‑150 (MQ=255)
                                                                                                  aCGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaaa                                               <  1:75816/150‑1 (MQ=255)
                                                                                                   cGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaaaa                                              >  1:589376/1‑150 (MQ=255)
                                                                                                    gAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAAc                                             >  1:60396/1‑150 (MQ=255)
                                                                                                      aGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACca                                           >  1:1239304/1‑150 (MQ=255)
                                                                                                      aGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACca                                           <  2:1821722/150‑1 (MQ=255)
                                                                                                         ttttCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAg                                        <  2:1395376/150‑1 (MQ=255)
                                                                                                           ttCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTa                                      <  2:1649934/150‑1 (MQ=255)
                                                                                                           ttCCTTATAAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTAAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTa                                      >  2:1831835/1‑150 (MQ=255)
                                                                                                                   cAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCt                              >  1:1572958/1‑150 (MQ=255)
                                                                                                                   cAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCt                              >  1:1566001/1‑150 (MQ=255)
                                                                                                                            cTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTa                     <  2:850889/150‑1 (MQ=255)
                                                                                                                              tCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAg                   <  1:419066/150‑1 (MQ=255)
                                                                                                                                 gACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAACACCACAGTAATCAGCCTTCTGAATTAAGTtg                >  2:1764473/1‑150 (MQ=255)
                                                                                                                                 gACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTtg                <  2:949039/150‑1 (MQ=255)
                                                                                                                                  aCTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTtgt               >  2:1551280/1‑150 (MQ=255)
                                                                                                                                     ttgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGct            >  1:182840/4‑150 (MQ=255)
                                                                                                                                     cGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGct            >  1:1396621/1‑150 (MQ=255)
                                                                                                                                      tgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGctc           >  1:1199046/3‑150 (MQ=255)
                                                                                                                                      tgcTTTCCATGTGCAAGAAAAGCAAGCAGCTGGAAAAACTCAATAAAATCCTGAACAATCACCAGAATAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGctc           >  2:1152292/3‑150 (MQ=255)
                                                                                                                                           tgcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGa      >  2:557495/1‑150 (MQ=255)
                                                                                                                                            gcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCGTCTGAATTAAGTTGTGCTCTTCGat     >  2:1601797/1‑150 (MQ=255)
                                                                                                                                             catgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGatc    >  2:1547129/1‑150 (MQ=255)
                                                                                                                                             catgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGatc    >  1:1007226/1‑150 (MQ=255)
                                                                                                                                               tgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGatcat  >  1:980988/1‑150 (MQ=255)
                                                                                                                                               tgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGatcat  >  2:1083607/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCAT  >  NC_000913/2810727‑2811019

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