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

TTAAACAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGC  >  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)
                     gcCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑Ca                                                                                                                             <  1:191130/150‑1 (MQ=255)
                      ccTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CAtc                                                                                                                           >  1:624518/1‑151 (MQ=255)
                         cGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CAtct                                                                                                                          <  2:1102853/149‑1 (MQ=255)
                          gAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCtt                                                                                                                       <  2:433183/151‑1 (MQ=255)
                          gAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCtt                                                                                                                       >  2:268591/1‑151 (MQ=255)
                           aaaTGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAACTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCttt                                                                                                                      <  2:441074/151‑1 (MQ=255)
                            aaTGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCttt                                                                                                                      >  1:734520/1‑150 (MQ=255)
                              tGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAAc                                                                                                                   >  2:1098181/1‑151 (MQ=255)
                               gAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAAcc                                                                                                                  <  2:994939/151‑1 (MQ=255)
                                aGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACca                                                                                                                 <  2:467002/151‑1 (MQ=255)
                                      ggAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAata                                                                                                           <  2:802259/151‑1 (MQ=255)
                                       gAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAatat                                                                                                          >  2:237744/1‑151 (MQ=255)
                                         aaaTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATgtgt                                                                                                        <  1:550766/151‑1 (MQ=255)
                                         aaaTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGt                                                                                                        >  1:841632/1‑151 (MQ=255)
                                         aaaTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGt                                                                                                        >  1:649790/1‑151 (MQ=255)
                                         aaaTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGt                                                                                                        <  1:1082468/151‑1 (MQ=255)
                                          aaTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGt                                                                                                        >  1:563984/1‑150 (MQ=255)
                                           aTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTa                                                                                                       >  1:420994/1‑150 (MQ=255)
                                            tAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTa                                                                                                       >  1:574123/1‑149 (MQ=255)
                                             aaGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTa                                                                                                       >  1:475143/1‑148 (MQ=255)
                                             aaGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAg                                                                                                     <  1:913828/150‑1 (MQ=255)
                                               ggTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGACTTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATg                                                                                                         >  1:505590/1‑144 (MQ=255)
                                                gTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGACGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGAt                                                                                                   >  2:242223/1‑149 (MQ=255)
                                                 ttATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGAtct                                                                                                 >  2:499600/1‑150 (MQ=255)
                                                  tATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGCAAGAtctc                                                                                                >  1:477360/1‑150 (MQ=255)
                                                   aTCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCa                                                                                               >  2:1170073/1‑150 (MQ=255)
                                                      aGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTa                                                                                           <  1:438357/151‑1 (MQ=255)
                                                       gCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTAt                                                                                          >  2:476277/1‑151 (MQ=255)
                                                           tGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGc                                                                                       <  2:658224/150‑1 (MQ=255)
                                                            gTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTT‑ACCAACAATTTGTAAAAttttaa                                                                                              >  1:402607/1‑136 (MQ=255)
                                                             ttttCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCAt                                                                                     <  2:914245/150‑1 (MQ=255)
                                                              tttCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCAt                                                                                     >  2:687446/1‑149 (MQ=255)
                                                              tttCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATc                                                                                    <  2:744353/150‑1 (MQ=255)
                                                                  tcCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTgg                                                                               >  1:162602/1‑151 (MQ=255)
                                                                    ccTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGt                                                                                 <  1:628252/147‑1 (MQ=255)
                                                                    ccTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTg                                                                                <  2:421412/148‑1 (MQ=255)
                                                                          tACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AAtt                                                                       <  1:832846/150‑1 (MQ=255)
                                                                              tGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCaa                                                                    <  2:314871/149‑1 (MQ=255)
                                                                                aaGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATtat                                                                <  2:934762/151‑1 (MQ=255)
                                                                                aaGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATta                                                                 <  2:530876/150‑1 (MQ=255)
                                                                                  ggATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATtata                                                               >  1:1142842/1‑150 (MQ=255)
                                                                                   gATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAAc                                                             >  2:1123453/1‑151 (MQ=255)
                                                                                        gAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTtctc                                                        >  2:103354/1‑151 (MQ=255)
                                                                                        gAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTtct                                                         <  2:1181868/150‑1 (MQ=255)
                                                                                        gAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTtct                                                         <  2:76712/150‑1 (MQ=255)
                                                                                         aaGCTTAAACCCTATGTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTGATTAAGAACAATTTTATCTC‑CATCCCTTTAACCAACAATATGTAAGAGCTCCGCTATCGCATCCGGAGATT‑AATTCAATtatgattt                                                           <  2:1085319/147‑6 (MQ=255)
                                                                                         aaGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTtctct                                                       >  1:555254/1‑151 (MQ=255)
                                                                                         aaGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTtctct                                                       <  2:276336/151‑1 (MQ=255)
                                                                                           gCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTctctct                                                     <  1:852481/151‑1 (MQ=255)
                                                                                              aaaaCCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAAc                                                  >  2:711591/1‑151 (MQ=255)
                                                                                                     ttttttATAAAGCATTTGTCCGAATTCGGACATATTAAAAAAAGTTTATTTAAGTTCAATTT‑ATTTACCACTTTTTTAACCAACATTATGAAAGTTCTAACCAATCCAATCCG‑GGATTAAATTCAATTATACCTTCCTTCTAAGCCtgt                                             >  1:999788/1‑149 (MQ=255)
                                                                                                      tttttATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTgtgt                                           >  2:1149110/1‑150 (MQ=255)
                                                                                                      tttttATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTgtgt                                           >  1:380275/1‑150 (MQ=255)
                                                                                                            taAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTa                                     >  1:529386/1‑150 (MQ=255)
                                                                                                            taAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTa                                     >  1:211104/1‑150 (MQ=255)
                                                                                                              aaGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACg                                  <  1:1062577/151‑1 (MQ=255)
                                                                                                               aGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACgg                                 <  1:1181415/151‑1 (MQ=255)
                                                                                                               aGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACgg                                 <  1:807409/151‑1 (MQ=255)
                                                                                                                 cattaatCCGACTTCGGACACAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCCATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGt                                >  2:1092694/7‑150 (MQ=255)
                                                                                                                 cATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGt                                >  1:170519/1‑150 (MQ=255)
                                                                                                                    ttGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAac                             >  2:1073032/1‑150 (MQ=255)
                                                                                                                    ttGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAac                             >  1:539357/1‑150 (MQ=255)
                                                                                                                     tGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAGCCAACAATATGTAAGATCTCAACTATCGCATGCGTGGATT‑AATTCAATTATAAATTCTCTCTAACGCAGTGTATCGTAACGGTAaca                            <  2:15643/150‑1 (MQ=255)
                                                                                                                          gAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGt                        >  1:829966/1‑149 (MQ=255)
                                                                                                                          gAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTag                      <  1:1044692/151‑1 (MQ=255)
                                                                                                                                ggACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGa                <  2:273878/151‑1 (MQ=255)
                                                                                                                                  aCATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGc              >  1:332079/1‑151 (MQ=255)
                                                                                                                                         cAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtg        >  2:357015/1‑150 (MQ=255)
                                                                                                                                         cAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtg        >  2:421370/1‑150 (MQ=255)
                                                                                                                                         cAAAAAAGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtg        >  1:287849/1‑150 (MQ=255)
                                                                                                                                          aaaaaaGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACATTATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgat      >  1:843322/1‑151 (MQ=255)
                                                                                                                                          aaaaaaGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgat      >  1:265324/1‑151 (MQ=255)
                                                                                                                                           aaaaaGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgatg     >  2:959376/1‑151 (MQ=255)
                                                                                                                                            aaaaGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGATCACATtgatg     >  2:1109851/1‑150 (MQ=255)
                                                                                                                                              aaGCTTCATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGc  >  1:884980/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTAAACAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTA‑CATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATT‑AATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGC  >  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: ‑