Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,250,989 G→C 5.6% intergenic (+35/‑35) 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,9890GC 5.6% 420.8 / 13.7 137intergenic (+35/‑35)yqjF/yqjGputative quinol oxidase subunit/putative S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (74/53);  new base C (4/4);  total (78/59)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

Only 100 of 146 total aligned reads displayed.

CGTCAACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGGC  >  NC_000913/3250840‑3251137
                                                                                                                                                     |                                                                                                                                                    
cGTCAACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCCATACTTAACGATAAAAAACGaggagg                                                                                                                                                      >  1:1240171/1‑150 (MQ=255)
 gTCAACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGCTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGTTAAGCACGCGCTATACTTAACGATAAAAAACGAggagga                                                                                                                                                     >  1:907771/1‑150 (MQ=255)
    aaCTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGc                                                                                                                                                  <  1:604028/150‑1 (MQ=255)
       tCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGctcc                                                                                                                                               <  1:1487145/150‑1 (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)
                 tCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGc                                                                                                                                     <  2:1725646/150‑1 (MQ=255)
                 tCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTCGc                                                                                                                                     >  1:1803868/1‑150 (MQ=255)
                  cATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCt                                                                                                                                    >  2:1920460/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)
                    tGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAt                                                                                                                                  >  1:1751708/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)
                              gACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGgagag                                                                                                                        >  1:2089195/1‑150 (MQ=255)
                                cAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGaa                                                                                                                      >  1:1840948/1‑150 (MQ=255)
                                   tttCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGTGGTGAGaaaaa                                                                                                                   >  2:1772727/1‑150 (MQ=255)
                                      cTGGCTGATTCCTGCTTCTTGCAATTACCGGTCCGGGCGCGTATATCTTCGACCGCCTTCTGAATAAATAGTGGTAAGCACGCGCTATTCTTAACTATAAATAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATgg                                                                                                                <  2:558103/150‑1 (MQ=255)
                                       tGGCGTATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCTCCTGATGAATAAAAAGTGGTAAGCACGCTCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGATGAGAGAAAAATggg                                                                                                               <  2:317312/150‑1 (MQ=255)
                                        ggCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGt                                                                                                              >  2:1793033/1‑150 (MQ=255)
                                         gCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTc                                                                                                             >  2:756733/1‑150 (MQ=255)
                                          cGGATTCCTGCTGCTGGCAATTACCGGTCCTGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCa                                                                                                            <  2:1353484/150‑1 (MQ=255)
                                            gATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAAc                                                                                                          <  1:1883390/150‑1 (MQ=255)
                                             aTTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACt                                                                                                         >  1:689052/1‑150 (MQ=255)
                                             aTTCCGGCTGCGGGCAATTACCGGTCCGGGCGCGTATAGCCTCGAACGCCTGCTGAATCAAAAGTGGGAAGCACGCGCTATACTTAAAGATAAAAAACGAGGCGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAAcc                                                                                                         <  1:899611/150‑2 (MQ=255)
                                                ccTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGAt                                                                                                      >  1:1959297/1‑150 (MQ=255)
                                                     tgctgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACg                                                                                                 >  2:1025511/1‑150 (MQ=255)
                                                      gctgGCAATTACCGGTCCGGGCGCGTTTAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACgg                                                                                                <  1:33764/150‑1 (MQ=255)
                                                      gctgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACgg                                                                                                <  1:15266/150‑1 (MQ=255)
                                                      gctgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACgg                                                                                                >  2:1209709/1‑150 (MQ=255)
                                                           cAATTACCGGTCCGGGCGCGTATAGCATCGCCCGCCTGCTGAATAAAAAGTGGTACGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGGCAACTGATTGACGGCGTCt                                                                                           >  1:944538/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)
                                                                 ccGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCCCGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATg                                                                                     >  2:1104634/1‑150 (MQ=255)
                                                                 ccGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATg                                                                                     >  2:1108318/1‑150 (MQ=255)
                                                                 ccGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATg                                                                                     >  2:1565801/1‑150 (MQ=255)
                                                                        gggCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTg                                                                              >  2:953580/1‑150 (MQ=255)
                                                                             cgTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACg                                                                         >  1:793755/1‑150 (MQ=255)
                                                                             cgTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCAAGACACCTGGTACg                                                                         <  2:952615/150‑1 (MQ=255)
                                                                                 taGCATCGACCGCCTGCTGAATAAAAAGTGGTACGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGGTATTGCTATTGGAGGAGAGAAAAATGGGTAAACTGATTGACGGCGTCTGGCATGACACCTCGTACGATAc                                                                     >  1:257744/1‑150 (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)
                                                                                       cGACAGCCTGCTGAATAAACAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGGCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATc                                                               <  1:398480/150‑1 (MQ=255)
                                                                                        gACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCt                                                              >  1:847609/1‑150 (MQ=255)
                                                                                        gACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCt                                                              >  1:56319/1‑150 (MQ=255)
                                                                                                  tgAATAAAAAGTGGTAAGCACGCGCTATACTTTACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTa                                                    >  2:238466/1‑150 (MQ=255)
                                                                                                  tgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACCCCTGGTACGATACCAAATCTACCGGCGGTa                                                    >  2:1856841/1‑150 (MQ=255)
                                                                                                  tgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTa                                                    >  1:1556693/1‑150 (MQ=255)
                                                                                                  tgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTa                                                    >  2:199943/1‑150 (MQ=255)
                                                                                                  tgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGAGACCAAATCTCCCGGCGGTa                                                    >  1:649250/1‑150 (MQ=255)
                                                                                                   gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTaa                                                   >  1:961647/1‑150 (MQ=255)
                                                                                                    aaTAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTaaa                                                  <  2:1870612/150‑1 (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)
                                                                                                        aaaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAttta                                              >  2:1658294/1‑149 (MQ=255)
                                                                                                            gTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACg                                          <  1:253239/150‑1 (MQ=255)
                                                                                                              ggTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGtt                                        >  1:995021/1‑150 (MQ=255)
                                                                                                              ggTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGtt                                        >  2:1479166/1‑150 (MQ=255)
                                                                                                               gTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTc                                       >  2:246820/1‑150 (MQ=255)
                                                                                                               gTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTc                                       <  1:199943/150‑1 (MQ=255)
                                                                                                               gTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTc                                       <  1:673880/150‑1 (MQ=255)
                                                                                                                tAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCa                                      >  1:589409/1‑150 (MQ=255)
                                                                                                                      cgcgCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTcc                                <  1:1910204/150‑1 (MQ=255)
                                                                                                                       gcgcTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGATAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCg                               <  2:2089195/150‑1 (MQ=255)
                                                                                                                            ataCTTAACGATAAAAAACGAGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTc                         >  2:1145235/1‑150 (MQ=255)
                                                                                                                             taCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTc                         <  1:1025511/150‑1 (MQ=255)
                                                                                                                             taCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGACGAGAGAAAAATGGGTCAACTGATTGACTGCGTCTGGCATGACACATGGTACGATACCAAATCTACCGGCGGAAAATTTAAACGTTCAGCTTCCGCaatac                         >  2:1284339/1‑146 (MQ=255)
                                                                                                                              aCTTATCTATAAAAATCGAGGATGATTCTACTCTTTTTTGCTATTGGAGGAGAGAAAAATGCGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCg                        <  1:123138/150‑1 (MQ=255)
                                                                                                                              aCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCg                        >  1:1447182/1‑150 (MQ=255)
                                                                                                                              aCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCg                        <  1:311678/150‑1 (MQ=255)
                                                                                                                              aCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCg                        <  2:1158372/150‑1 (MQ=255)
                                                                                                                              aCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCg                        <  1:740097/150‑1 (MQ=255)
                                                                                                                              aCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTGCAGCTTCCGCATTTCg                        <  1:1857927/150‑1 (MQ=255)
                                                                                                                               cTTAACTATAAAAAACGAGGAGGAAGAGCCTCGTTTTTGCTATTGGATGAGAGAAAAATGGGTCAACTGATAGACGTCGTCTGGAATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGt                       <  1:366633/150‑1 (MQ=255)
                                                                                                                               cTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGt                       <  2:1467481/150‑1 (MQ=255)
                                                                                                                                   aCGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACt                   >  1:217967/1‑150 (MQ=255)
                                                                                                                                      agcaaaaaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGc                >  1:843230/4‑150 (MQ=255)
                                                                                                                                      aTAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGc                <  1:1565801/150‑1 (MQ=255)
                                                                                                                                         aaaaaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGCTTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAGCGTTCAGCTACCGCATTTCGTAACTGGCTcc             >  1:1569432/1‑149 (MQ=255)
                                                                                                                                          aaaaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCAc            >  1:1582719/1‑150 (MQ=255)
                                                                                                                                             aCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGc         <  2:1240145/150‑1 (MQ=255)
                                                                                                                                              ctagtaggaATTTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGTTCAACTGATTGACGGATGCTGTCATTACACCTTGTACGATACCATATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGct        <  1:116199/145‑2 (MQ=255)
                                                                                                                                              cGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGcc        >  2:681083/1‑150 (MQ=255)
                                                                                                                                              cGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGcc        <  2:274637/150‑1 (MQ=255)
                                                                                                                                                cggagcaAGCGCCTCGTTTTTGCTATTGGAGGCGCGAAAACTGGGGCACCTGATTGACGGCGTCTGGCATGACACCTGGAACGATCCCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGa      <  2:1266107/149‑1 (MQ=255)
                                                                                                                                                   aggaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATgg   <  2:982169/150‑1 (MQ=255)
                                                                                                                                                    ggaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGgc  >  1:2066051/1‑150 (MQ=255)
                                                                                                                                                    ggaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGgc  <  2:1751708/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGTCAACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGGC  >  NC_000913/3250840‑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: