Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,906,229 T→G 9.2% I181S (ATC→AGC) ‡ mntP → Mn(2(+)) exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,906,2290TG 9.2% 252.2 / 16.5 87I181S (ATC→AGC) ‡mntPMn(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base T (39/40);  new base G (5/3);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGC  >  NC_000913/1906080‑1906376
                                                                                                                                                     |                                                                                                                                                   
cAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAGAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAt                                                                                                                                                     >  1:1164859/1‑150 (MQ=255)
cAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAt                                                                                                                                                     >  2:630295/1‑150 (MQ=255)
 aaCATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATc                                                                                                                                                    <  1:829356/150‑1 (MQ=255)
   cATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGGTGATCGGCATCGGCGTCCAGATCct                                                                                                                                                  >  1:1303060/1‑150 (MQ=255)
    aTTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCctc                                                                                                                                                 <  1:389329/150‑1 (MQ=255)
    aTTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCctc                                                                                                                                                 >  2:516299/1‑150 (MQ=255)
        tCGCGACCTCCTTGGCCATTGGTTTTGCAACCTTGATTATGTCTACTTTAGTGATGATGGTTGGTCGCTTTATCGTCTCAATTATTGGGAAAAACGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGa                                                                                                                                             <  1:997248/150‑1 (MQ=255)
        tCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAGGATCTGGa                                                                                                                                             <  1:73665/150‑1 (MQ=255)
        ccgcgACCGCATTGGCCATTGGGGGTGCAACCTTGATTATGTCAACATTAGGGATGACGGTTGGGCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGa                                                                                                                                             <  1:36109/149‑1 (MQ=255)
           cgACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAGGATCTGGACGc                                                                                                                                          <  1:225780/150‑1 (MQ=255)
            gACCGCATTGGCCATTGGTTGTGAAACCTTGATTATGTCAACATTAGGTATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAGGATCTGGACGcc                                                                                                                                         <  2:475686/150‑2 (MQ=255)
            gACCGCATTGGCCATTGCTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCa                                                                                                                                         <  1:547218/150‑1 (MQ=255)
              ccGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACt                                                                                                                                       <  1:385918/150‑1 (MQ=255)
                       ccATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGTAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGtt                                                                                                                              <  2:1474799/150‑1 (MQ=255)
                        cATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGATCAATTATTGGCAAAAAAGCTGAAATTCACGTCGGGATGCTGCTGATCGGCATCGGCGTCCAGATACTCTGGACGCACTTCCACGGTta                                                                                                                             >  1:1490463/1‑150 (MQ=255)
                         aTTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaa                                                                                                                            <  1:1449216/150‑1 (MQ=255)
                          ttGGTTGTGCAACGTTGATTATGTCAACTTTATGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAATAAGCGGAAATTCTCGGATGTGTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGTACGCACTTCCACGGTtaat                                                                                                                           <  2:252412/150‑1 (MQ=255)
                           tGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaata                                                                                                                          >  1:764095/1‑150 (MQ=255)
                            ggTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaataa                                                                                                                         >  2:783447/1‑150 (MQ=255)
                            ggTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaataa                                                                                                                         <  2:972248/150‑1 (MQ=255)
                               tgtgCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGaa                                                                                                                      >  1:890056/1‑150 (MQ=255)
                               tgtgCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGaa                                                                                                                      >  2:1164215/1‑150 (MQ=255)
                                  gCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGc                                                                                                                   >  2:237154/1‑150 (MQ=255)
                                  gCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGc                                                                                                                   >  2:231674/1‑150 (MQ=255)
                                  gCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGc                                                                                                                   <  1:704378/150‑1 (MQ=255)
                                    aaCCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTg                                                                                                                 >  1:248183/1‑150 (MQ=255)
                                     aCCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGc                                                                                                                >  1:334299/1‑150 (MQ=255)
                                      ccTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGcc                                                                                                               <  2:614288/150‑1 (MQ=255)
                                      ccTTGATTATGTCAACATTAGGGAGGATGGTTGGTCGCTTTATCGGCTCAATTATTGGTAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGcc                                                                                                               <  1:544641/150‑1 (MQ=255)
                                       cTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCa                                                                                                              >  1:813470/1‑150 (MQ=255)
                                       cTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCa                                                                                                              >  2:629392/1‑150 (MQ=255)
                                             tATTTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTg                                                                                                        <  1:783447/150‑1 (MQ=255)
                                                 tCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATcc                                                                                                                                                   <  2:1207193/103‑1 (MQ=255)
                                                 tCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATcc                                                                                                                                                   >  1:1207193/1‑103 (MQ=255)
                                                  cAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAatat                                                                                                   >  1:483881/1‑150 (MQ=255)
                                                    aCATTAGGCATGAGGGTTGGGCGATTTATCGGCTCAATTATTGGGAATTAAGCGGAAATTTTCGGCGGCCTGGTGCTGATCGGCTTGTGCGTCCAGATCCTCTGGACTCACTTCCACGGTTAATAAGATCGCTGCCAGAGGTTAATATTa                                                                                                 <  2:1166757/150‑1 (MQ=255)
                                                     cATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTaa                                                                                                >  1:965309/1‑150 (MQ=255)
                                                     cATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTaa                                                                                                <  2:179880/150‑1 (MQ=255)
                                                     cATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTaa                                                                                                <  1:308254/150‑1 (MQ=255)
                                                               gatgGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTg                                                                                      <  2:1076308/150‑1 (MQ=255)
                                                                 tggttggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAATTTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGc                                                                                    <  1:1004542/150‑1 (MQ=255)
                                                                 tggttggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGc                                                                                    <  1:1506130/150‑1 (MQ=255)
                                                                 tggttggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGc                                                                                    <  1:1500899/150‑1 (MQ=255)
                                                                  ggttggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCa                                                                                   >  2:62381/1‑150 (MQ=255)
                                                                            tttATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACAct                                                                         <  1:222037/150‑1 (MQ=255)
                                                                              tATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACActtc                                                                       >  2:204322/1‑150 (MQ=255)
                                                                                tCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACActtctt                                                                     >  1:439526/1‑150 (MQ=255)
                                                                                    cTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTtgc                                                                 <  2:439526/150‑1 (MQ=255)
                                                                                    cTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTtgc                                                                 <  2:248183/150‑1 (MQ=255)
                                                                                      cAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTaa                                                                                                <  1:1461477/117‑1 (MQ=255)
                                                                                      cAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTaa                                                                                                >  2:1461477/1‑117 (MQ=255)
                                                                                       aaTTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTtgctgc                                                              <  1:629392/150‑1 (MQ=255)
                                                                                       aaTTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTtgctgc                                                              <  2:869637/150‑1 (MQ=255)
                                                                                             tGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAgtgt                                                        >  1:766801/1‑150 (MQ=255)
                                                                                              gggAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGAAGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGtt                                                       >  1:262602/1‑150 (MQ=255)
                                                                                              gggAAAAAAGCGGAAATTCTCGGCGGCCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTACACGCTTTATAAGAACGCTGCCAGAGCTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGAGtt                                                       >  2:351510/1‑150 (MQ=255)
                                                                                               ggAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGttt                                                      >  2:782817/1‑150 (MQ=255)
                                                                                                gAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTg                                                     >  2:268393/1‑150 (MQ=255)
                                                                                                  aaaaaGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGcc                                                   <  1:516299/150‑1 (MQ=255)
                                                                                                   aaaaGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGGCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCa                                                  <  2:811609/150‑1 (MQ=255)
                                                                                                    aaaGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAg                                                 >  2:1503682/1‑150 (MQ=255)
                                                                                                      aGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGAc                                               >  2:1090560/1‑150 (MQ=255)
                                                                                                      aGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGAc                                               <  1:1379916/150‑1 (MQ=255)
                                                                                                      aGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGAc                                               >  2:372458/1‑150 (MQ=255)
                                                                                                          gAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCt                                           >  1:607459/1‑150 (MQ=255)
                                                                                                          gAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCt                                           >  1:155446/1‑150 (MQ=255)
                                                                                                           aaaTTCTCGGCGGGCTGGTGCGGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTg                                          <  1:1256770/150‑1 (MQ=255)
                                                                                                              ttCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCt                                       >  2:1460713/1‑150 (MQ=255)
                                                                                                                  cggcggGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGc                                   >  1:870755/1‑150 (MQ=255)
                                                                                                                        gCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCa                             <  1:204322/150‑1 (MQ=255)
                                                                                                                            gTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCa                         >  2:1120531/1‑150 (MQ=255)
                                                                                                                               cTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAAc                      >  1:938218/1‑150 (MQ=255)
                                                                                                                                 gATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACgg                    >  1:674212/1‑150 (MQ=255)
                                                                                                                                    cggcatcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGt                 <  1:1190404/150‑1 (MQ=255)
                                                                                                                                     ggcatcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGTCAGTCGAACACTTCTTTTGCTGCCAGTGTTTTCCAGACTTCTTGATTCTCACGCTAGGCAAAGgtcgtc                <  2:334299/150‑1 (MQ=255)
                                                                                                                                      gcatcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCa               <  1:925248/150‑1 (MQ=255)
                                                                                                                                      gcatcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCa               <  1:372458/150‑1 (MQ=255)
                                                                                                                                         agtgcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCt            >  1:966454/4‑150 (MQ=255)
                                                                                                                                          cggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTg           <  2:1164859/150‑1 (MQ=255)
                                                                                                                                          cggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGGCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTg           <  2:965309/150‑1 (MQ=255)
                                                                                                                                           tgcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc          >  1:427902/2‑150 (MQ=255)
                                                                                                                                           tgcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc          >  2:1090294/2‑150 (MQ=255)
                                                                                                                                           ggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc          <  1:630295/150‑1 (MQ=255)
                                                                                                                                            gcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCa         >  2:65179/1‑150 (MQ=255)
                                                                                                                                            gcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCa         >  1:524877/1‑150 (MQ=255)
                                                                                                                                            gcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCa         >  1:355889/1‑150 (MQ=255)
                                                                                                                                                ccAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCaataa     >  1:1280053/1‑150 (MQ=255)
                                                                                                                                                   gATCCTCTGGACGCACTTCCACGGTTAATAAGACCGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGc  >  2:1031101/1‑150 (MQ=255)
                                                                                                                                                   gATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGc  >  2:1067651/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGC  >  NC_000913/1906080‑1906376

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