Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,906,231 C→A 7.5% L182I (CTC→ATC)  mntP → Mn(2(+)) exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,906,2310CA 7.5% 266.5 / 8.4 80L182I (CTC→ATC) mntPMn(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base C (30/44);  new base A (3/3);  total (33/47)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

ACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAAT  >  NC_000913/1906082‑1906374
                                                                                                                                                     |                                                                                                                                               
aCATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATcc                                                                                                                                                 >  1:1238117/1‑150 (MQ=255)
 cATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCct                                                                                                                                                <  1:1214112/150‑1 (MQ=255)
  aTTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATGCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCctc                                                                                                                                               <  1:1068254/150‑1 (MQ=255)
     aTCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTCTTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTgg                                                                                                                                            <  1:1610813/150‑1 (MQ=255)
     aTCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTgg                                                                                                                                            <  1:616829/150‑1 (MQ=255)
     aTCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTgg                                                                                                                                            >  2:57042/1‑150 (MQ=255)
      tCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGa                                                                                                                                           >  2:862552/1‑150 (MQ=255)
         cgACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGc                                                                                                                                        <  2:1235746/150‑1 (MQ=255)
         cgACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAGGATCTGGACGc                                                                                                                                        <  2:414747/150‑1 (MQ=255)
         cgACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGGGCTGATCGGCATCGGCGTCCAGAGGATCTGGACGc                                                                                                                                        <  1:204926/150‑1 (MQ=255)
            ccGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTCTTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACt                                                                                                                                     <  1:1035749/150‑1 (MQ=255)
            ccGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTCTTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGGGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACt                                                                                                                                     <  1:1041656/150‑1 (MQ=255)
             cGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCCCtt                                                                                                                                    >  2:648661/1‑150 (MQ=255)
              gCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTc                                                                                                                                   >  1:534917/1‑150 (MQ=255)
              gCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTc                                                                                                                                   >  1:1011040/1‑150 (MQ=255)
               cATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGAGGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTcc                                                                                                                                  <  1:912968/150‑1 (MQ=255)
               cATTGGCCATTGGGTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGGTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTcc                                                                                                                                  <  1:918973/150‑1 (MQ=255)
                   ggCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACgg                                                                                                                              >  2:263922/1‑150 (MQ=255)
                   ggCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACgg                                                                                                                              >  1:443147/1‑150 (MQ=255)
                    gCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGt                                                                                                                             >  1:478483/1‑150 (MQ=255)
                     ccATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGtt                                                                                                                            <  1:338071/150‑1 (MQ=255)
                         tGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaata                                                                                                                        <  1:1060734/150‑1 (MQ=255)
                             tgtgCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGaa                                                                                                                    <  2:370575/150‑1 (MQ=255)
                               tgCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACg                                                                                                                  <  1:1253953/150‑1 (MQ=255)
                               tgCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACg                                                                                                                  >  1:409005/1‑150 (MQ=255)
                                 cAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCt                                                                                                                >  2:1295563/1‑150 (MQ=255)
                                 cAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCt                                                                                                                >  2:1350712/1‑150 (MQ=255)
                                         aTTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAgg                                                                                                        <  2:1439533/150‑1 (MQ=255)
                                           tATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTg                                                                                                      >  1:1148843/1‑150 (MQ=255)
                                           tATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCGGGTGCTGATCGGCAGCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTg                                                                                                      <  2:173403/150‑1 (MQ=255)
                                           tATGGCAACATTAGGGCTGCTGGGTGGGCCCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTg                                                                                                      <  2:539545/150‑1 (MQ=255)
                                               tCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAata                                                                                                  <  2:963430/150‑1 (MQ=255)
                                               tCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAata                                                                                                  <  1:949704/150‑1 (MQ=255)
                                                  aCATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTa                                                                                               <  2:962337/150‑1 (MQ=255)
                                                  aCATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTa                                                                                               <  1:57042/150‑1 (MQ=255)
                                                   cATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTaa                                                                                              <  1:940963/150‑1 (MQ=255)
                                                   cATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTaa                                                                                              <  1:445687/150‑1 (MQ=255)
                                                     ttAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAg                                                                                            <  1:617037/150‑1 (MQ=255)
                                                     ttAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAg                                                                                            <  2:770233/150‑1 (MQ=255)
                                                             gatgGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTg                                                                                    >  2:453667/1‑150 (MQ=255)
                                                                ggttggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCa                                                                                 >  1:1469241/1‑150 (MQ=255)
                                                                    ggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCg                                                                             <  1:255701/150‑1 (MQ=255)
                                                                     gtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGa                                                                            <  1:1386057/150‑1 (MQ=255)
                                                                      tCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGaa                                                                           <  2:534917/150‑1 (MQ=255)
                                                                       cGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGGGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAac                                                                          <  2:775551/150‑1 (MQ=255)
                                                                         cTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAacac                                                                        >  1:1401602/1‑150 (MQ=255)
                                                                         cTTTATCCGAACAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGCCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGAAAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAacac                                                                        >  2:599882/1‑150 (MQ=255)
                                                                          tttATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACAct                                                                       <  2:1238117/150‑1 (MQ=255)
                                                                            tATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACActtc                                                                     >  2:1636824/1‑150 (MQ=255)
                                                                               cGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCttt                                                                  <  2:1529241/150‑1 (MQ=255)
                                                                               cGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCttt                                                                  >  1:663972/1‑150 (MQ=255)
                                                                                        tattGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAg                                                         >  2:543083/1‑150 (MQ=255)
                                                                                         attGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAgt                                                        >  2:1651384/1‑150 (MQ=255)
                                                                                           tGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTtg                                                                >  2:1425603/1‑140 (MQ=255)
                                                                                           tGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTtg                                                                <  1:1425603/140‑1 (MQ=255)
                                                                                            gggAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGtt                                                     <  2:409005/150‑1 (MQ=255)
                                                                                                aaacaGCGGAAATTCTCGGCGGGCTGGCGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGcc                                                 <  2:1351873/150‑1 (MQ=255)
                                                                                                aaaaaGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGcc                                                 <  2:443147/150‑1 (MQ=255)
                                                                                                aaaaaGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGcc                                                 <  2:993305/150‑1 (MQ=255)
                                                                                                         aaaTTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTg                                        <  1:1350712/150‑1 (MQ=255)
                                                                                                          aaTTCTCGGCGGGCGGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGGCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTgg                                       <  2:1222466/150‑1 (MQ=255)
                                                                                                           aTTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGc                                      <  1:134938/150‑1 (MQ=255)
                                                                                                            ttCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCt                                     >  1:237586/1‑150 (MQ=255)
                                                                                                               tcGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCg                                  <  2:454706/150‑1 (MQ=255)
                                                                                                                cggcggGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGc                                 <  2:237586/150‑1 (MQ=255)
                                                                                                                 ggcggGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCa                                >  1:1043555/1‑150 (MQ=255)
                                                                                                                       cTGGTGCTGAGCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAg                          <  1:391629/150‑1 (MQ=255)
                                                                                                                         ggTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTT‑CCAGACTTCTGGCTTCGCACGCCAGGCa                       >  2:1201082/1‑150 (MQ=255)
                                                                                                                           tGCTGATCGGCATCGGCGTCCACATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCaa                      >  2:298135/1‑150 (MQ=255)
                                                                                                                            gCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCaaa                     >  2:1411311/1‑150 (MQ=255)
                                                                                                                             cTGATCGGCATCGGCGGCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGGGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAAc                    <  2:1463577/150‑1 (MQ=255)
                                                                                                                                    gcatcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTTTCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCa             <  2:949226/150‑1 (MQ=255)
                                                                                                                                     gaagtgcGTCCAGAGGATCTGGACGCCCTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAg                          >  1:394714/6‑136 (MQ=255)
                                                                                                                                     gaagtgcGTCCAGAGGATCTGGACGCCCTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAg                          <  2:394714/131‑1 (MQ=255)
                                                                                                                                       tcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCt          >  1:888465/1‑150 (MQ=255)
                                                                                                                                        gtgcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTg         >  1:1089308/3‑150 (MQ=255)
                                                                                                                                         tgcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAATGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc        >  2:1033182/2‑150 (MQ=255)
                                                                                                                                         ggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc        <  2:1011040/150‑1 (MQ=255)
                                                                                                                                          gcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCa       <  1:648661/150‑1 (MQ=255)
                                                                                                                                          gcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCa       >  1:78525/1‑150 (MQ=255)
                                                                                                                                             tCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCaata    <  1:453667/150‑1 (MQ=255)
                                                                                                                                               cAGATCCTCTGGACGCACTTGCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGCCTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCaacaat  >  2:375167/1‑150 (MQ=255)
                                                                                                                                               cAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCaataat  >  1:1044937/1‑150 (MQ=255)
                                                                                                                                                    ccTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc        >  2:695457/1‑139 (MQ=255)
                                                                                                                                                    ccTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc        <  1:695457/139‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
ACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAAT  >  NC_000913/1906082‑1906374

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: ‑