Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,810,925 G→A G52G (GGG→GGA mprA → DNA‑binding transcriptional regulator MprA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,810,9250GA95.6% 146.6 / ‑1.2 45G52G (GGG→GGAmprADNA‑binding transcriptional regulator MprA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (17/26);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

ACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCAT  >  NC_000913/2810679‑2811168
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
aCCCCCAAGAATCATTTTTCTAAAACAATACATTTCCTTTATTTGTCCCTGTCGTTACTAGATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCAAAAATTCCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCATGGAATTa                                                                                                                                                                                                                                                 <  2:53084/251‑1 (MQ=255)
              tttttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGttt                                                                                                                                                                                                                                   >  2:288971/1‑251 (MQ=255)
                     aaaaCAATACATTTTCTTTATTTGTCACTGTCGTTACTATATCGGCTTAAATTAATGAGGTCATCCCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTCCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTTAGGGAATTAACGAGACGTTGTTTATgg                                                                                                                                                                                                                               <  1:212234/248‑1 (MQ=255)
                       aaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTg                                                                                                                                                                                                                          >  2:294780/1‑251 (MQ=255)
                          aaTACATTTACTTTATTTGTCACGTTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAAAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGAtt                                                                                                                                                                                                                       <  2:460637/251‑1 (MQ=255)
                                        aTTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtctc                                                                                                                                                                                                          <  2:329076/250‑1 (MQ=255)
                                         tttGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtctc                                                                                                                                                                                                          >  1:309309/1‑249 (MQ=255)
                                           tGTCACGGTCGGTACTAAAGCGGCTGAAATTAAGGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATGTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg                                                                                                                                                                                                       <  2:554613/250‑1 (MQ=255)
                                                            tatCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTCTGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGc                                                                                                                                                                                     >  1:213023/1‑251 (MQ=255)
                                                              ccGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACCCCCATTAAACAAATGCTAAAATTTCGCGCCATCCGCCACGAAGATTTTCCGTATCAGGAGATCCTTCTGACTCGGCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGGTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATGCAGCCt                                                                                                                                                                                   <  2:551001/250‑1 (MQ=255)
                                                                gtcTGAAATTAATAAGGTCATACCCAAAGGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGATCCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTc                                                                                                                                                                                 <  2:291924/249‑1 (MQ=255)
                                                                      aaTTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAAtt                                                                                                                                                                           <  1:385948/251‑1 (MQ=255)
                                                                       aTTAATGAGGTCATACCCACATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCGGGAAAACCACAGTATTAAGCCTTCTGAAtt                                                                                                                                                                           >  2:538613/1‑250 (MQ=255)
                                                                           aTGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGtt                                                                                                                                                                      <  1:504700/251‑1 (MQ=255)
                                                                             gAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTtg                                                                                                                                                                     <  1:34412/250‑1 (MQ=255)
                                                                                    tACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTg                                                                                                                                                             >  1:478995/1‑251 (MQ=255)
                                                                                      ccaaaaatgaaagtctgttgaCGCCCATTGAACAAATGATAAAATTTCTCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACACGCAAAGCAAGCTGCGGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAACCCACAGTATTCAGCCTTCGGAATTAAGTTGTGCTCTTGGa                                                                                                                                                           <  2:395323/231‑1 (MQ=255)
                                                                                                                    aaCAAATGCGAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGcc                                                                                                                             <  1:567210/251‑1 (MQ=255)
                                                                                                                       aaaTGCTAAAATTCCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGACTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGAt                                                                                                                          <  2:324707/251‑1 (MQ=255)
                                                                                                                          tGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGcctctga                                                                                                                        >  2:326426/1‑245 (MQ=255)
                                                                                                                             tAAAATTTCGCGCCAGCCGCCTCGAAGATGTTCCTTATCAGGAGATCCTTTTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGcc                                                                                                                                                                                    <  1:323952/187‑1 (MQ=255)
                                                                                                                             tAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGCCTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGCACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGcc                                                                                                                                                                                    >  2:323952/1‑187 (MQ=255)
                                                                                                                             aaaaaTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTg                                                                                                                    <  1:317828/250‑1 (MQ=255)
                                                                                                                                         ccAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGTAAAAACGCGGt                                                                                                        >  1:235628/1‑251 (MQ=255)
                                                                                                                                          agggcgccACGAAGATATGCCGTATCAGGAGATCCTTCTGACTCGTCGTTGCATGCACAGGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGATGAAGGCGCAGGGAATTAACGAGACGTTGTTTAGGGCGTTGATTACGCTGGAGTCTCAGGAAAGCCACAGTATTCAGCCTTCTGAATTAAGTTGTGGTCTTGGAGCAACCCGTTCCAACGCGAGGCGTAGTGCCGATGAACTGGAAAAACGCGGtt                                                                                                       <  2:57719/247‑1 (MQ=255)
                                                                                                                                                 cACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCg                                                                                                 >  1:312299/1‑250 (MQ=255)
                                                                                                                                                             ccTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGaaa                                                                                     <  1:266598/250‑1 (MQ=255)
                                                                                                                                                                 aacaggatgtCCTTCTGACTCGTCTTTGCATGTACATTCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTTGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGAt                                                                                <  2:493083/242‑1 (MQ=255)
                                                                                                                                                                 aTCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGCCGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGa                                                                                 >  2:461495/1‑250 (MQ=255)
                                                                                                                                                                 aTCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGAt                                                                                >  1:544973/1‑251 (MQ=255)
                                                                                                                                                                                 gACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGctgcc                                                                >  1:215868/1‑251 (MQ=255)
                                                                                                                                                                                          tggcctgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCGGAAGGCTCAGGGAATTAACGAGACGTTGTGTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCGGCATCt                                                         <  2:253888/244‑1 (MQ=255)
                                                                                                                                                                                                        aaaGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCAACCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGaaaaa                                         <  2:364952/251‑1 (MQ=255)
                                                                                                                                                                                                           gcaagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGt                                      <  1:461495/251‑1 (MQ=255)
                                                                                                                                                                                                             aagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTTAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAACCCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTGGGATCATCCCGTACCAACGGGACGCGGATTGCCGATGAACTGGAAAAACGCGGTTGGACCGAACTCCGTGAAAGCGATAACGATCGCCGCTGCCTGCCTCTGCAATTAACGGAAAAAGGTCa                                    <  1:325720/251‑1 (MQ=255)
                                                                                                                                                                                                              agcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCAc                                   <  1:406461/251‑1 (MQ=255)
                                                                                                                                                                                                                  gctgGAGAACCGCAATAAAATGCGGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAg                                <  2:150664/250‑1 (MQ=255)
                                                                                                                                                                                                                         aaCCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAACCGTTCCTTATACCGCCGATTAAGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCCTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGATGGATCGAACGTCGGGCAACCGATAACGATCCCCGCTGACTGACTCTGAAATTAAAGGAACAAGGACATGAAtttttccg                        >  2:243319/1‑248 (MQ=255)
                                                                                                                                                                                                                                aTAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGtt                 <  1:294780/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTa              <  2:20155/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTaccacc         >  2:345505/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATAATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTACCGGAAAAAGGTCACGAGTTTTTGCGCAAGGTTTTaccacc         >  1:343041/1‑251 (MQ=255)
                                                                                                                                                                                                                                         tGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCg        >  1:133092/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGAAGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTaccacagc       >  1:47517/1‑246 (MQ=255)
                                                                                                                                                                                                                                                 cAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCAt  <  1:517544/249‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
ACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCAT  >  NC_000913/2810679‑2811168

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: