Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,133,696 A→C intergenic (+190/‑139) metF → / → katG 5,10‑methylenetetrahydrofolate reductase/catalase‑peroxidase HPI, heme b‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,133,6960AC98.8% 326.4 / NA 84intergenic (+190/‑139)metF/katG5,10‑methylenetetrahydrofolate reductase/catalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (48/34);  total (49/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTAAACAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGC  >  NC_000913/4133548‑4133840
                                                                                                                                                    |                                                                                                                                                  
ttAAACAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCttcat                                                                                                                                                  >  2:999625/1‑151 (MQ=255)
 tAAACAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCttcat                                                                                                                                                  <  2:1074943/150‑1 (MQ=255)
     cAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGa                                                                                                                                             <  1:213754/151‑1 (MQ=255)
      agagCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCttcatta                                                                                                                                                >  1:273889/1‑147 (MQ=255)
         gCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCaa                                                                                                                                         <  1:332631/151‑1 (MQ=255)
            caAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAAttt                                                                                                                                      >  1:66763/1‑151 (MQ=255)
             aaaaTGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTg                                                                                                                                     >  1:934943/1‑151 (MQ=255)
             aaaaTGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTg                                                                                                                                     >  2:515078/1‑151 (MQ=255)
             aaaaTGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTg                                                                                                                                     >  2:1175407/1‑151 (MQ=255)
               aaTGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGAt                                                                                                                                   <  1:55963/151‑1 (MQ=255)
                aTGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGCCATAATCAAAAAAGTTTCATTAAGATCAAtt                                                                                                                                       <  2:197955/146‑1 (MQ=255)
                  gctgcCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCt                                                                                                                                 >  1:790464/1‑150 (MQ=255)
                     gcCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACa                                                                                                                             <  1:191130/150‑1 (MQ=255)
                      ccTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACAtc                                                                                                                           >  1:624518/1‑151 (MQ=255)
                         cGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACAtct                                                                                                                          <  2:1102853/149‑1 (MQ=255)
                          gAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCtt                                                                                                                       <  2:433183/151‑1 (MQ=255)
                          gAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCtt                                                                                                                       >  2:268591/1‑151 (MQ=255)
                           aaaTGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAACTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCttt                                                                                                                      <  2:441074/151‑1 (MQ=255)
                            aaTGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCttt                                                                                                                      >  1:734520/1‑150 (MQ=255)
                              tGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAAc                                                                                                                   >  2:1098181/1‑151 (MQ=255)
                               gAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAAcc                                                                                                                  <  2:994939/151‑1 (MQ=255)
                                aGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACca                                                                                                                 <  2:467002/151‑1 (MQ=255)
                                      ggAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAata                                                                                                           <  2:802259/151‑1 (MQ=255)
                                       gAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAatat                                                                                                          >  2:237744/1‑151 (MQ=255)
                                         aaaTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATgtgt                                                                                                        <  1:550766/151‑1 (MQ=255)
                                         aaaTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGt                                                                                                        >  1:841632/1‑151 (MQ=255)
                                         aaaTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGt                                                                                                        >  1:649790/1‑151 (MQ=255)
                                         aaaTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGt                                                                                                        <  1:1082468/151‑1 (MQ=255)
                                          aaTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGt                                                                                                        >  1:563984/1‑150 (MQ=255)
                                           aTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTa                                                                                                       >  1:420994/1‑150 (MQ=255)
                                            tAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTa                                                                                                       >  1:574123/1‑149 (MQ=255)
                                             aaGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTa                                                                                                       >  1:475143/1‑148 (MQ=255)
                                             aaGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAg                                                                                                     <  1:913828/150‑1 (MQ=255)
                                               ggTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGACTTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATg                                                                                                         >  1:505590/1‑144 (MQ=255)
                                                gTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGACGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGAt                                                                                                   >  2:242223/1‑149 (MQ=255)
                                                 ttATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGAtct                                                                                                 >  2:499600/1‑150 (MQ=255)
                                                  tATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGCAAGAtctc                                                                                                >  1:477360/1‑150 (MQ=255)
                                                   aTCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCa                                                                                               >  2:1170073/1‑150 (MQ=255)
                                                      aGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTa                                                                                           <  1:438357/151‑1 (MQ=255)
                                                       gCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTAt                                                                                          >  2:476277/1‑151 (MQ=255)
                                                           tGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGc                                                                                       <  2:658224/150‑1 (MQ=255)
                                                            gTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTACCAACAATTTGTAAAAttttaa                                                                                              >  1:402607/1‑136 (MQ=255)
                                                             ttttCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCAt                                                                                     <  2:914245/150‑1 (MQ=255)
                                                              tttCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCAt                                                                                     >  2:687446/1‑149 (MQ=255)
                                                              tttCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATc                                                                                    <  2:744353/150‑1 (MQ=255)
                                                                  tcCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTgg                                                                               >  1:162602/1‑151 (MQ=255)
                                                                    ccTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGt                                                                                 <  1:628252/147‑1 (MQ=255)
                                                                    ccTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTg                                                                                <  2:421412/148‑1 (MQ=255)
                                                                          tACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAAtt                                                                       <  1:832846/150‑1 (MQ=255)
                                                                              tGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCaa                                                                    <  2:314871/149‑1 (MQ=255)
                                                                                aaGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATtat                                                                <  2:934762/151‑1 (MQ=255)
                                                                                aaGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATta                                                                 <  2:530876/150‑1 (MQ=255)
                                                                                  ggATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATtata                                                               >  1:1142842/1‑150 (MQ=255)
                                                                                   gATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAAc                                                             >  2:1123453/1‑151 (MQ=255)
                                                                                        gAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTtctc                                                        >  2:103354/1‑151 (MQ=255)
                                                                                        gAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTtct                                                         <  2:1181868/150‑1 (MQ=255)
                                                                                        gAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTtct                                                         <  2:76712/150‑1 (MQ=255)
                                                                                         aaGCTTAAACCCTATGTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTGATTAAGAACAATTTTATCTCCATCCCTTTAACCAACAATATGTAAGAGCTCCGCTATCGCATCCGGAGATTAATTCAATtatgattt                                                           <  2:1085319/147‑6 (MQ=255)
                                                                                         aaGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTtctct                                                       >  1:555254/1‑151 (MQ=255)
                                                                                         aaGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTtctct                                                       <  2:276336/151‑1 (MQ=255)
                                                                                           gCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTctctct                                                     <  1:852481/151‑1 (MQ=255)
                                                                                              aaaaCCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAAc                                                  >  2:711591/1‑151 (MQ=255)
                                                                                                     ttttttATAAAGCATTTGTCCGAATTCGGACATATTAAAAAAAGTTTATTTAAGTTCAATTTATTTACCACTTTTTTAACCAACATTATGAAAGTTCTAACCAATCCAATCCGGGATTAAATTCAATTATACCTTCCTTCTAAGCCtgt                                             >  1:999788/1‑149 (MQ=255)
                                                                                                      tttttATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTgtgt                                           >  2:1149110/1‑150 (MQ=255)
                                                                                                      tttttATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTgtgt                                           >  1:380275/1‑150 (MQ=255)
                                                                                                            taAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTa                                     >  1:529386/1‑150 (MQ=255)
                                                                                                            taAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTa                                     >  1:211104/1‑150 (MQ=255)
                                                                                                              aaGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACg                                  <  1:1062577/151‑1 (MQ=255)
                                                                                                               aGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACgg                                 <  1:1181415/151‑1 (MQ=255)
                                                                                                               aGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACgg                                 <  1:807409/151‑1 (MQ=255)
                                                                                                                 cattaatCCGACTTCGGACACAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCCATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGt                                >  2:1092694/7‑150 (MQ=255)
                                                                                                                 cATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGt                                >  1:170519/1‑150 (MQ=255)
                                                                                                                    ttGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAac                             >  2:1073032/1‑150 (MQ=255)
                                                                                                                    ttGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAac                             >  1:539357/1‑150 (MQ=255)
                                                                                                                     tGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAGCCAACAATATGTAAGATCTCAACTATCGCATGCGTGGATTAATTCAATTATAAATTCTCTCTAACGCAGTGTATCGTAACGGTAaca                            <  2:15643/150‑1 (MQ=255)
                                                                                                                          gAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGt                        >  1:829966/1‑149 (MQ=255)
                                                                                                                          gAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTag                      <  1:1044692/151‑1 (MQ=255)
                                                                                                                                ggACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGa                <  2:273878/151‑1 (MQ=255)
                                                                                                                                  aCATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGc              >  1:332079/1‑151 (MQ=255)
                                                                                                                                         cAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtg        >  2:357015/1‑150 (MQ=255)
                                                                                                                                         cAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtg        >  2:421370/1‑150 (MQ=255)
                                                                                                                                         cAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtg        >  1:287849/1‑150 (MQ=255)
                                                                                                                                          aaaaaaGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACATTATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgat      >  1:843322/1‑151 (MQ=255)
                                                                                                                                          aaaaaaGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgat      >  1:265324/1‑151 (MQ=255)
                                                                                                                                           aaaaaGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgatg     >  2:959376/1‑151 (MQ=255)
                                                                                                                                            aaaaGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGATCACATtgatg     >  2:1109851/1‑150 (MQ=255)
                                                                                                                                              aaGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGc  >  1:884980/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTAAACAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGC  >  NC_000913/4133548‑4133840

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