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.9% 342.0 / ‑4.1 87intergenic (+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 (57/29);  total (58/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

TAAACAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACG  >  NC_000913/4133549‑4133843
                                                                                                                                                   |                                                                                                                                                   
tAAACAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCttcatt                                                                                                                                                  >  2:402571/1‑151 (MQ=255)
 aaaCAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCttcatt                                                                                                                                                  >  2:1190177/1‑150 (MQ=255)
 aaaCAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCttcatt                                                                                                                                                  <  2:154108/150‑1 (MQ=255)
  aaCAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCttcattaa                                                                                                                                                <  2:1003991/151‑1 (MQ=255)
  aaCAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCttcattaa                                                                                                                                                <  2:688924/151‑1 (MQ=255)
   aCAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAg                                                                                                                                               <  2:653944/151‑1 (MQ=255)
       agCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATc                                                                                                                                            >  1:1176790/1‑150 (MQ=255)
        gCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCaa                                                                                                                                          <  1:435941/151‑1 (MQ=255)
        gCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCa                                                                                                                                           >  1:204001/1‑150 (MQ=255)
         cacaAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAAt                                                                                                                                         >  1:669383/1‑151 (MQ=255)
           caAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAAttt                                                                                                                                       >  2:68375/1‑151 (MQ=255)
             aaaTGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATCATCAAAAAAGCTTCATTAAGATCAATTTg                                                                                                                                      >  2:734222/1‑150 (MQ=255)
             aaaTGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTg                                                                                                                                      >  1:251678/1‑150 (MQ=255)
               aTGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATc                                                                                                                                   >  1:1212136/1‑151 (MQ=255)
                 gctgcCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTa                                                                                                                                 <  2:1182413/151‑1 (MQ=255)
                    gcCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACAt                                                                                                                              <  2:412542/151‑1 (MQ=255)
                     ccTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACAtc                                                                                                                             <  2:815623/151‑1 (MQ=255)
                       tCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACAtctc                                                                                                                           >  2:899082/1‑151 (MQ=255)
                        cGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACAtct                                                                                                                            <  1:66766/149‑1 (MQ=255)
                        cGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATctct                                                                                                                          <  1:743351/151‑1 (MQ=255)
                             tGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAAc                                                                                                                     >  1:54339/1‑151 (MQ=255)
                              gAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAAcc                                                                                                                    >  1:571189/1‑151 (MQ=255)
                                  gCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACcaaca                                                                                                                <  2:65229/151‑1 (MQ=255)
                                   cGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACcaaca                                                                                                                >  2:582007/1‑150 (MQ=255)
                                        aaaTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGt                                                                                                          >  1:315411/1‑151 (MQ=255)
                                         aaTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGt                                                                                                          >  1:106103/1‑150 (MQ=255)
                                          aTAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTa                                                                                                         >  2:908917/1‑150 (MQ=255)
                                           tAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGCTCTACATCTCTTTAACCAACAATATGTa                                                                                                         >  2:739501/1‑149 (MQ=255)
                                           tAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTa                                                                                                         >  1:1257199/1‑149 (MQ=255)
                                            aaGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTa                                                                                                         >  1:310797/1‑148 (MQ=255)
                                            aaGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAg                                                                                                       <  2:546899/150‑1 (MQ=255)
                                             aGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTa                                                                                                         >  1:1024749/1‑147 (MQ=255)
                                              ggTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTa                                                                                                         >  2:363062/1‑146 (MQ=255)
                                                     aGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACt                                                                                              >  1:898848/1‑150 (MQ=255)
                                                       ccTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATc                                                                                           <  2:104406/151‑1 (MQ=255)
                                                        cTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCg                                                                                          >  1:705006/1‑151 (MQ=255)
                                                              ttCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATcc                                                                                     <  1:568900/150‑1 (MQ=255)
                                                              ttCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATcc                                                                                     <  1:39975/150‑1 (MQ=255)
                                                              ttCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCg                                                                                    >  2:345601/1‑151 (MQ=255)
                                                               tctcCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCg                                                                                    <  1:423802/150‑1 (MQ=255)
                                                               tctcCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGt                                                                                   >  1:1102754/1‑151 (MQ=255)
                                                               tctcCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATGCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCg                                                                                    <  1:1081905/150‑1 (MQ=255)
                                                                ctcCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTg                                                                                  >  2:485245/1‑151 (MQ=255)
                                                                 tcCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTgg                                                                                 >  1:759831/1‑151 (MQ=255)
                                                                  cccTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGAGCTACATCTCTGGAACCAACAAGAGGGAAGATCTCAACTATCGCATCCGTgt                                                                                 <  1:276063/150‑2 (MQ=255)
                                                                   ccTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGAt                                                                               >  2:820078/1‑151 (MQ=255)
                                                                     tCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGAtt                                                                              >  2:1111945/1‑150 (MQ=255)
                                                                      cATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTa                                                                             <  1:48604/150‑1 (MQ=255)
                                                                      cATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTa                                                                             >  1:1102163/1‑150 (MQ=255)
                                                                         tACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTc                                                                         >  2:915976/1‑151 (MQ=255)
                                                                          aCTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTc                                                                         <  1:1273150/150‑1 (MQ=255)
                                                                           cTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCa                                                                        >  1:985478/1‑150 (MQ=255)
                                                                           cTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCa                                                                        >  1:596847/1‑150 (MQ=255)
                                                                             tGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCaa                                                                       <  2:781523/149‑1 (MQ=255)
                                                                                 ggATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATtata                                                                  >  2:219636/1‑150 (MQ=255)
                                                                                      tGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTtct                                                            >  1:125795/1‑151 (MQ=255)
                                                                                       gAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTtctc                                                           <  1:12649/151‑1 (MQ=255)
                                                                                        aaGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTtctct                                                          >  2:59481/1‑151 (MQ=255)
                                                                                                aCCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCt                                                  <  1:885161/151‑1 (MQ=255)
                                                                                                 cccTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCtg                                                 <  2:962781/151‑1 (MQ=255)
                                                                                                  ccTTTTTTATCAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGCTTAATTCAATTATAACTTCTCTCTAACGCtgt                                                >  1:1018738/1‑151 (MQ=255)
                                                                                                  ccTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGCGGATTAATTCAATTATAACTTCTCTCTAACGCtgt                                                >  1:240987/1‑151 (MQ=255)
                                                                                                   cTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCtgtg                                               >  2:225083/1‑151 (MQ=255)
                                                                                                     tttttATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTgtgt                                              >  1:1214839/1‑150 (MQ=255)
                                                                                                       tttATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATc                                           >  1:1080762/1‑151 (MQ=255)
                                                                                                        ttATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCg                                          >  1:71279/1‑151 (MQ=255)
                                                                                                         tataAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGt                                         >  2:1240050/1‑151 (MQ=255)
                                                                                                          ataAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGt                                         <  1:10883/150‑1 (MQ=255)
                                                                                                          ataAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTa                                        <  2:571205/151‑1 (MQ=255)
                                                                                                           taAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTa                                        >  2:178924/1‑150 (MQ=255)
                                                                                                           taAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTa                                        >  1:384656/1‑150 (MQ=255)
                                                                                                           taAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTa                                        >  2:25065/1‑150 (MQ=255)
                                                                                                            aaaGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTaa                                       <  2:1260120/150‑1 (MQ=255)
                                                                                                            aaaGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAAc                                      >  1:336979/1‑151 (MQ=255)
                                                                                                             aaGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACg                                     <  1:1269514/151‑1 (MQ=255)
                                                                                                               gCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGt                                   <  2:204003/151‑1 (MQ=255)
                                                                                                                cATTTGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTa                                  <  1:698058/151‑1 (MQ=255)
                                                                                                                    tGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAacac                              >  2:605836/1‑151 (MQ=255)
                                                                                                                    tGTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAaca                               >  1:620806/1‑150 (MQ=255)
                                                                                                                     gTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACt                             >  1:320449/1‑151 (MQ=255)
                                                                                                                     gTCCGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACt                             >  1:1140467/1‑151 (MQ=255)
                                                                                                                        cGAATTCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGt                           >  1:501436/1‑150 (MQ=255)
                                                                                                                          aaTTCGGACATAATCAAAAAAGCTGCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGt                           >  2:1158369/1‑148 (MQ=255)
                                                                                                                             tCGGACATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAgg                      <  2:1122835/150‑1 (MQ=255)
                                                                                                                                  cATAATCAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGca                <  1:384382/151‑1 (MQ=255)
                                                                                                                                        cAAAAAAGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtg           >  2:1018713/1‑150 (MQ=255)
                                                                                                                                         aaaaaaGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgat         >  2:104305/1‑151 (MQ=255)
                                                                                                                                         aaaaaaGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgat         >  1:518996/1‑151 (MQ=255)
                                                                                                                                          aaaaaGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgatg        >  1:530438/1‑151 (MQ=255)
                                                                                                                                          aaaaaGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAGATACCAACTTTCGCAATCGTGGATTAAGTAAATTTTAAACTCTTCCTAACGCTGGGGATCCGAACGGTAAAAATGTTAAGGGGAAcaatttaatgg       >  2:475188/1‑142 (MQ=255)
                                                                                                                                           aaaaGCTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgatg        >  1:282026/1‑150 (MQ=255)
                                                                                                                                                cTTCATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACg  >  1:380837/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TAAACAGAGCACAAAATGCTGCCTCGAAATGAGGGCGGGAAAATAAGGTTATCAGCCTTGTTTTCTCCCTCATTACTTGAAGGATATGAAGCTAAAACCCTTTTTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACG  >  NC_000913/4133549‑4133843

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