Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,250,991 A→T 5.2% intergenic (+37/‑33) yqjF → / → yqjG putative quinol oxidase subunit/putative S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,250,9910AT 5.2% 429.9 / 8.6 144intergenic (+37/‑33)yqjF/yqjGputative quinol oxidase subunit/putative S‑transferase
Reads supporting (aligned to +/- strand):  ref base A (77/59);  new base T (3/5);  total (80/64)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 144 total aligned reads displayed.

AACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGGC  >  NC_000913/3250844‑3251137
                                                                                                                                                   |                                                                                                                                                  
   tCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGctcc                                                                                                                                               <  1:1487145/150‑1 (MQ=255)
    cGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAGGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGctccc                                                                                                                                              >  1:426463/1‑149 (MQ=255)
    cGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGGAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGctcct                                                                                                                                              >  1:1158372/1‑150 (MQ=255)
         aTGTTCATGAAAAACCTGACAATTTCTGGCTGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGttt                                                                                                                                         <  2:802563/150‑1 (MQ=255)
         aTGTGCATGAAAAACCTGACAATTTCTGTCGGATTCCTGCTGCTGGCAATTACCGGGCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGttt                                                                                                                                         <  2:880007/150‑1 (MQ=255)
          tGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGtttt                                                                                                                                        <  2:1877667/150‑1 (MQ=255)
          tGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGtttt                                                                                                                                        <  2:397896/150‑1 (MQ=255)
          tGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTAAACGATAAAAAACGAGGAGGAAGCTCCTCGtttt                                                                                                                                        <  2:272828/150‑1 (MQ=255)
           gTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGgtttt                                                                                                                                       >  2:1761510/1‑150 (MQ=255)
             tCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGc                                                                                                                                     <  2:117953/150‑1 (MQ=255)
             tCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGc                                                                                                                                     <  2:1725646/150‑1 (MQ=255)
             tCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTCGc                                                                                                                                     >  1:1803868/1‑150 (MQ=255)
               aTGAAAAACCTGACAATTTCTGGCTGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGttttttatc                                                                                                                                   <  2:821947/150‑5 (MQ=255)
               aTGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTa                                                                                                                                   >  1:747099/1‑150 (MQ=255)
                tGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAt                                                                                                                                  >  2:1910204/1‑150 (MQ=255)
                 gAAACACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCTGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAtt                                                                                                                                 <  1:1130798/150‑1 (MQ=255)
                   aaaaCCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTTTTgg                                                                                                                               >  1:564740/1‑150 (MQ=255)
                   aaaaCCTGACAATATCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTgg                                                                                                                               <  1:519200/150‑1 (MQ=255)
                    aaaCCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTgga                                                                                                                              <  1:899080/150‑1 (MQ=255)
                      aCCTGACAATTTCTGGCGGAATCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCGGCTGAATAAAAAGTGGTAAGCACGCGCTATTCTTAACGATAAAAAACGAGGAGTAAGCTCCTCGTTTTTGCTATTggagg                                                                                                                            <  1:1573672/150‑1 (MQ=255)
                          gACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGTGCGCGTATAGCATCGACCGACTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCACCTCGTTTTTGCTATTGGAGgagac                                                                                                                        >  2:1272104/1‑149 (MQ=255)
                            cAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGgagagag                                                                                                                      >  1:934512/1‑149 (MQ=255)
                            cAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGaa                                                                                                                      >  1:1840948/1‑150 (MQ=255)
                               tttCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGTGGTGAGaaaaa                                                                                                                   >  2:1772727/1‑150 (MQ=255)
                                  cTGGCTGATTCCTGCTTCTTGCAATTACCGGTCCGGGCGCGTATATCTTCGACCGCCTTCTGAATAAATAGTGGTAAGCACGCGCTATTCTTAACTATAAATAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATgg                                                                                                                <  2:558103/150‑1 (MQ=255)
                                  cTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATgg                                                                                                                <  1:251221/150‑1 (MQ=255)
                                   tGGCGTATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCTCCTGATGAATAAAAAGTGGTAAGCACGCTCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGATGAGAGAAAAATggg                                                                                                               <  2:317312/150‑1 (MQ=255)
                                    ggCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGt                                                                                                              >  1:1706407/1‑150 (MQ=255)
                                      cGGATTCCTGCTGCTGGCAATTCCCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAATAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCa                                                                                                            <  2:1463149/150‑1 (MQ=255)
                                      cGGATTCCTGCTGCTGGCAATTACCGGTCCTGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCa                                                                                                            <  2:1353484/150‑1 (MQ=255)
                                        gATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAAc                                                                                                          <  1:1883390/150‑1 (MQ=255)
                                         aTTCCTGCTGCTGGCAATTACCGGTCCGTGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCTCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACt                                                                                                         <  1:1494842/150‑1 (MQ=255)
                                         aTTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACt                                                                                                         >  1:982169/1‑150 (MQ=255)
                                         aTTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACt                                                                                                         >  1:689052/1‑150 (MQ=255)
                                         aTTCCGGCTGCGGGCAATTACCGGTCCGGGCGCGTATAGCCTCGAACGCCTGCTGAATCAAAAGTGGGAAGCACGCGCTATACTTAAAGATAAAAAACGAGGCGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAAcc                                                                                                         <  1:899611/150‑2 (MQ=255)
                                            ccTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGAt                                                                                                      >  1:1959297/1‑150 (MQ=255)
                                                ctgctgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGAc                                                                                                  <  2:77225/150‑1 (MQ=255)
                                                 tgctgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACg                                                                                                 >  2:1025511/1‑150 (MQ=255)
                                                  gctgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACgg                                                                                                >  2:1209709/1‑150 (MQ=255)
                                                           tACCGGTCCGGGCGCGGATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGATGAGGAAGCTCCGCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCa                                                                                       <  1:836283/150‑1 (MQ=255)
                                                           tACCGGTACGGGCGCTTATATCATCGACCGGCTGCTTAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCa                                                                                       <  2:276862/150‑1 (MQ=255)
                                                           tACCGGGCCGGGCGCGGAAAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGATGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCa                                                                                       <  1:1019639/150‑1 (MQ=255)
                                                             ccGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGTAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCAtt                                                                                     >  2:1524332/1‑149 (MQ=255)
                                                             ccGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATg                                                                                     >  2:1108318/1‑150 (MQ=255)
                                                             ccGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATg                                                                                     >  2:1565801/1‑150 (MQ=255)
                                                                    gggCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTg                                                                              >  2:953580/1‑150 (MQ=255)
                                                                    gggCGCGTATAGCATCGACCGCCTGATGAATAAAATGTGTTAAGCACGCGCTATACTTTTAGTTAAATTACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAATATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTg                                                                              <  2:1936501/150‑1 (MQ=255)
                                                                     ggCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTgg                                                                             >  2:1137005/1‑150 (MQ=255)
                                                                         cgTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACg                                                                         >  1:793755/1‑150 (MQ=255)
                                                                         cgTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCAAGACACCTGGTACg                                                                         <  2:952615/150‑1 (MQ=255)
                                                                               gCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCa                                                                   >  2:519668/1‑150 (MQ=255)
                                                                               gCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGAAACCa                                                                   >  2:728757/1‑150 (MQ=255)
                                                                                 aTCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAAATGATTGACGGCGTCTGGCATGACACCTGGTACGATACCaaa                                                                 >  2:895266/1‑150 (MQ=255)
                                                                                 aTCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGGTTTTGCTATTGGAGGAGATAAAAATGGGACAACTGATAGACGGCGGATGGCATGACACCTGGTACGATACCaaa                                                                 >  1:64004/1‑150 (MQ=255)
                                                                                    gACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCt                                                              >  1:847609/1‑150 (MQ=255)
                                                                                    gACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCt                                                              >  1:56319/1‑150 (MQ=255)
                                                                                      ccGCCTGCTGAATTAAAAGTGGTAAGCACGCGCTATAATTAACGATAAAACACGAGGAGGAAGCTCGTCGTTTTTGCTAATGGAGGAGAGAAAAATGGTTCAACTGATTGACGGCGACTGGCTTGACACATGGTACGATACCAATTCTAc                                                            >  1:186721/1‑150 (MQ=255)
                                                                                            gctgAATAAAAAGTGGTCAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACcggcgg                                                      <  1:1540952/150‑1 (MQ=255)
                                                                                              tgAATAAAAAGTGGTAAGCACGCGCTATACTTTACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTa                                                    >  2:238466/1‑150 (MQ=255)
                                                                                              tgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACCCCTGGTACGATACCAAATCTACCGGCGGTa                                                    >  2:1856841/1‑150 (MQ=255)
                                                                                              tgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTa                                                    >  2:199943/1‑150 (MQ=255)
                                                                                              tgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTa                                                    <  1:650926/150‑1 (MQ=255)
                                                                                              tgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTa                                                    >  1:1556693/1‑150 (MQ=255)
                                                                                               gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAATTGGGTCAACTGATTGACGGCGTATGGCATGACACCTGGTACGATACGAAATCTACCGGCGGTaa                                                   >  2:236575/1‑150 (MQ=255)
                                                                                               gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTaa                                                   >  1:961647/1‑150 (MQ=255)
                                                                                                 aTAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAt                                                 >  1:596309/1‑150 (MQ=255)
                                                                                                  aaaaaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGGCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAtt                                                <  1:1515896/149‑1 (MQ=255)
                                                                                                   aaaaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAttt                                               >  2:977006/1‑150 (MQ=255)
                                                                                                    aaaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAttta                                              >  1:432852/1‑149 (MQ=255)
                                                                                                    aaaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAttta                                              >  1:1673158/1‑149 (MQ=255)
                                                                                                    aaaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTc                                              <  1:280908/150‑1 (MQ=255)
                                                                                                        gTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACg                                          <  1:253239/150‑1 (MQ=255)
                                                                                                          ggTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGtt                                        >  2:1479166/1‑150 (MQ=255)
                                                                                                           gTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTc                                       <  1:673880/150‑1 (MQ=255)
                                                                                                           gTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTc                                       >  2:1814938/1‑150 (MQ=255)
                                                                                                            tAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCa                                      >  1:589409/1‑150 (MQ=255)
                                                                                                              aGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCATCTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATATCAACGTTCAGc                                    >  1:338447/1‑150 (MQ=255)
                                                                                                                    cgcTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGc                              <  1:2053043/150‑1 (MQ=255)
                                                                                                                        ataCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCAttt                          <  2:711379/150‑1 (MQ=255)
                                                                                                                        ataCTTAACGATAAAAAACGAGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTc                         >  2:1145235/1‑150 (MQ=255)
                                                                                                                         taCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTc                         <  1:1025511/150‑1 (MQ=255)
                                                                                                                         taCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGGCAACTGATTGTCGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTc                         <  1:1258675/150‑1 (MQ=255)
                                                                                                                          aCTTATCTATAAAAATCGAGGATGATTCTACTCTTTTTTGCTATTGGAGGAGAGAAAAATGCGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCg                        <  1:123138/150‑1 (MQ=255)
                                                                                                                          aCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCg                        <  1:311678/150‑1 (MQ=255)
                                                                                                                          aCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCg                        <  1:740097/150‑1 (MQ=255)
                                                                                                                           cTTAACTATAAAAAACGAGGAGGAAGAGCCTCGTTTTTGCTATTGGATGAGAGAAAAATGGGTCAACTGATAGACGTCGTCTGGAATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGt                       <  1:366633/150‑1 (MQ=255)
                                                                                                                               aCGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACt                   >  1:217967/1‑150 (MQ=255)
                                                                                                                                cGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTg                  >  1:1950248/1‑150 (MQ=255)
                                                                                                                                  agcaaaaaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGc                >  1:843230/4‑150 (MQ=255)
                                                                                                                                     aaaaaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCACCGTTCAGCTTCCGCATTTCGTAACTGGCTCa             >  2:1055619/1‑150 (MQ=255)
                                                                                                                                     aaaaaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGCTTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAGCGTTCAGCTACCGCATTTCGTAACTGGCTcc             >  1:1569432/1‑149 (MQ=255)
                                                                                                                                     aaaaaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCa             >  2:1906469/1‑150 (MQ=255)
                                                                                                                                         aCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGc         >  1:1810171/1‑150 (MQ=255)
                                                                                                                                          ctagtaggaATTTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGTTCAACTGATTGACGGATGCTGTCATTACACCTTGTACGATACCATATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGct        <  1:116199/145‑2 (MQ=255)
                                                                                                                                          cGAGGAGTAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGTTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTTCCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGcc        <  1:627873/150‑1 (MQ=255)
                                                                                                                                          cGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGcc        <  2:274637/150‑1 (MQ=255)
                                                                                                                                            cggagcaAGCGCCTCGTTTTTGCTATTGGAGGCGCGAAAACTGGGGCACCTGATTGACGGCGTCTGGCATGACACCTGGAACGATCCCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGa      <  2:1266107/149‑1 (MQ=255)
                                                                                                                                             ggaggaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGAt     >  2:1675168/1‑150 (MQ=255)
                                                                                                                                               aggaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATgg   <  2:982169/150‑1 (MQ=255)
                                                                                                                                                ggaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGgc  >  1:2066051/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGGC  >  NC_000913/3250844‑3251137

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