Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,250,990 A→T 5.1% intergenic (+36/‑34) 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,9900AT 5.1% 318.3 / 3.0 107intergenic (+36/‑34)yqjF/yqjGputative quinol oxidase subunit/putative S‑transferase
Reads supporting (aligned to +/- strand):  ref base A (57/44);  new base T (2/4);  total (59/48)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

TCAACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATG  >  NC_000913/3250842‑3251135
                                                                                                                                                    |                                                                                                                                                 
tCAACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGaa                                                                                                                                                  <  2:562198/150‑1 (MQ=255)
tCAACTCGCTGATGTTCATGAAAAACCTGAAAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGaa                                                                                                                                                  <  2:1816618/150‑1 (MQ=255)
 cAACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAg                                                                                                                                                 >  2:468074/1‑150 (MQ=255)
  aaCTCGCTGATGTTCATGAAAAACCTGACAATGTCTGGCGGATTCCTGCTGCTGGCAATTACCGTTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTCGTAAGCACGCGCTATACTTAACGATAAAAAAAGAGGAGGAAGc                                                                                                                                                >  2:145662/1‑150 (MQ=255)
      cGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGctcct                                                                                                                                            >  2:1180534/1‑150 (MQ=255)
            tGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGtttt                                                                                                                                      <  1:212449/150‑1 (MQ=255)
             gTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGttttt                                                                                                                                     >  1:1025851/1‑150 (MQ=255)
             gTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGttttt                                                                                                                                     <  1:1257415/150‑1 (MQ=255)
               tCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTTCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGc                                                                                                                                   <  1:1808431/150‑1 (MQ=255)
                cATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCt                                                                                                                                  <  2:1754109/150‑1 (MQ=255)
                cATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCt                                                                                                                                  <  1:842517/150‑1 (MQ=255)
                cATGAAAAACCTGACAATTTCTGGCGGAATCCTGCTGCTCGCAATTACCGATCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTACCGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCt                                                                                                                                  >  1:552058/1‑150 (MQ=255)
                  tGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGGCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGttttttatcg                                                                                                                                <  2:947790/150‑6 (MQ=255)
                  tGAAAAACCTGACAATTTCTGCCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAt                                                                                                                                >  2:1662579/1‑150 (MQ=255)
                    aaacaCCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTg                                                                                                                              <  1:851749/150‑1 (MQ=255)
                    aaaaaCCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTg                                                                                                                              <  1:1766839/150‑1 (MQ=255)
                     aaaaCCTGACAATATCAGGCGGAGTCCTGCTGCTGGCAATTACCGGGCCGGGCGAGTATCGCATCGACCGCCTGCTGAATAAAAAGTGGGATGCCCGCGCTATACTTAACGAAAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTgg                                                                                                                             <  1:1201587/150‑1 (MQ=255)
                        aCCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTggagg                                                                                                                          <  1:1346807/150‑1 (MQ=255)
                          cTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGgaggag                                                                                                                        >  1:1330796/1‑150 (MQ=255)
                          cTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGCCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGgaggag                                                                                                                        >  1:1714375/1‑150 (MQ=255)
                           tGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGgaga                                                                                                                       >  2:1666665/1‑150 (MQ=255)
                            gACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCGGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGgagag                                                                                                                      <  2:1229191/150‑1 (MQ=255)
                             aCAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGgagaga                                                                                                                     >  2:885105/1‑150 (MQ=255)
                               aaTTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGaaa                                                                                                                   >  2:705680/1‑150 (MQ=255)
                                  ttCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAAt                                                                                                                >  1:1395415/1‑150 (MQ=255)
                                   tCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATg                                                                                                               >  1:715115/1‑150 (MQ=255)
                                   tCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATg                                                                                                               <  2:595119/150‑1 (MQ=255)
                                     tGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATggg                                                                                                             <  2:855052/150‑1 (MQ=255)
                                     tGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATggg                                                                                                             <  1:428842/150‑1 (MQ=255)
                                       gCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTc                                                                                                           <  1:396445/150‑1 (MQ=255)
                                          gATTCCTGCTGCTGGCAATTACCGCTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAAc                                                                                                        <  1:164622/150‑1 (MQ=255)
                                             tCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGa                                                                                                     >  1:795645/1‑150 (MQ=255)
                                              ccagctgctgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGAt                                                                                                    <  2:1616865/147‑1 (MQ=255)
                                              ccTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGAt                                                                                                    >  2:1250932/1‑150 (MQ=255)
                                               ctgctgctgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTTGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGAtt                                                                                                   <  2:473007/150‑1 (MQ=255)
                                                tgctgctgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTg                                                                                                  >  1:396256/1‑150 (MQ=255)
                                                 gctgctgGCAATTACCGGTCCGGGCGCGTATATCATCGACCGCCTGCTGAATAAAAAGAGGTAAGCACGCGCTATACTTAACGCTAAAAAACGAGGAGGAAGCTGCTCGTTGTTGCTATTGGAGTAGAGAAAAATGGGTATACTGGTTGa                                                                                                 >  1:1195257/1‑150 (MQ=255)
                                                       ggCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGt                                                                                           >  1:383982/1‑150 (MQ=255)
                                                        gCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACCTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTATGCTATTGGAGGAGAGAAAAATGTGTCAACTGATTGACTGCTTc                                                                                          >  2:907920/1‑150 (MQ=255)
                                                        gCAATTACCGGTCAGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGAGTc                                                                                          >  1:415299/1‑150 (MQ=255)
                                                          aaTTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTg                                                                                        >  1:216728/1‑150 (MQ=255)
                                                          aaTTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTg                                                                                        >  2:631170/1‑150 (MQ=255)
                                                           aTTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCTCTATACTTAACGATAAAAAACGCGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTgg                                                                                       <  2:1289505/150‑1 (MQ=255)
                                                            ttACCGGTGCGGGCGATTATATCATCGACCGCTTTCTGAAATAAAATTGGTATGAACGCTCAATAATTAACGATAAAAACGGAGGAGGTAGCTCCTCGTTTTTGCTATTGGAGGAGCTAAAACTGGGTCAACTGATTGACGGCGTCTgtc                                                                                      <  2:894420/150‑3 (MQ=255)
                                                             tACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCa                                                                                     <  2:1742758/150‑1 (MQ=255)
                                                              aCCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCAt                                                                                    <  2:1395415/150‑1 (MQ=255)
                                                               ccGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATg                                                                                   <  1:1047153/150‑1 (MQ=255)
                                                               ccGGTCCGGGCGCGTATACCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTAAACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATg                                                                                   <  1:1655843/150‑1 (MQ=255)
                                                                    ccGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACAcc                                                                              >  1:1259566/1‑150 (MQ=255)
                                                                       ggCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTgg                                                                           >  2:1211984/1‑150 (MQ=255)
                                                                          gcttataTCATCGACCGACTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGATAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTAc                                                                        <  1:646950/147‑1 (MQ=255)
                                                                          gcgTATAGGATCGACCGCCTGCTGAATAAAACGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGGCAACTGATTGACGGCGTCTGGCATGACACCTGGTAc                                                                        <  2:45942/150‑1 (MQ=255)
                                                                                 gCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCa                                                                 >  2:602868/1‑150 (MQ=255)
                                                                                 gCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCa                                                                 >  1:1590583/1‑150 (MQ=255)
                                                                                  cATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCaa                                                                >  1:1771742/1‑150 (MQ=255)
                                                                                  cATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCaa                                                                >  2:1216580/1‑150 (MQ=255)
                                                                                  cATAGACAGCATGCTGAATAAAAAGTGGTAAGCACGCGATATACTTAACGATAAAAAACTAGTAGGAAGCTCCGCGATGTTGATATTGGAGGAGAGAAAAATGGGTCAACTGATTGACTGAGTCTGGCATGACACCTGGTACGATAACaa                                                                >  2:836771/1‑150 (MQ=255)
                                                                                   aTCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCaaa                                                               >  2:583996/1‑150 (MQ=255)
                                                                                   aTCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCaaa                                                               >  2:1042847/1‑150 (MQ=255)
                                                                                      gACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCt                                                            <  1:1396358/150‑1 (MQ=255)
                                                                                          gCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACcg                                                        <  2:400077/150‑1 (MQ=255)
                                                                                             tgctgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACcggcg                                                     <  1:1506578/150‑1 (MQ=255)
                                                                                              gctgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACcggcgg                                                    <  1:854132/150‑1 (MQ=255)
                                                                                                tgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTa                                                  >  1:738418/1‑150 (MQ=255)
                                                                                                 gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTaa                                                 >  2:514795/1‑150 (MQ=255)
                                                                                                    tAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAtt                                              >  1:1234317/1‑150 (MQ=255)
                                                                                                    tAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAtt                                              <  1:682711/150‑1 (MQ=255)
                                                                                                      aaaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTc                                            <  1:649933/150‑1 (MQ=255)
                                                                                                      aaaaGTGGTAAGCAAGCGCTATACTTAACGATAAAAAACGACGAGGAAGCTCCTCGTATTTGCTATTGGAGGAGAGAAAAAACGGTCAACTGAATGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTc                                            >  2:1337023/1‑150 (MQ=255)
                                                                                                       aaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCa                                           >  1:596600/1‑150 (MQ=255)
                                                                                                        aaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCaa                                          <  2:1550940/150‑1 (MQ=255)
                                                                                                         aGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAAc                                         <  1:885105/150‑1 (MQ=255)
                                                                                                          gTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACg                                        >  1:138624/1‑150 (MQ=255)
                                                                                                            ggTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGtt                                      >  2:1392380/1‑150 (MQ=255)
                                                                                                             gTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGtta                                     >  1:1596623/1‑149 (MQ=255)
                                                                                                              tAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCa                                    >  2:365192/1‑150 (MQ=255)
                                                                                                                 gCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCt                                 >  1:1796972/1‑150 (MQ=255)
                                                                                                                 gCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCt                                 >  1:1289134/1‑150 (MQ=255)
                                                                                                                 gCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCt                                 <  2:1275759/150‑1 (MQ=255)
                                                                                                                     gcgcTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCg                             <  1:1180534/150‑1 (MQ=255)
                                                                                                                      cgcTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGc                            >  1:1018236/1‑150 (MQ=255)
                                                                                                                      cgcTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGc                            <  1:1666665/150‑1 (MQ=255)
                                                                                                                       gcTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCa                           <  1:1662579/150‑1 (MQ=255)
                                                                                                                          ataCTTAACGATAAAAAACGAGGAGGAAGCTCCTAGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGAATGACGGCGTCTGGCTTGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCAttt                        <  1:312905/150‑1 (MQ=255)
                                                                                                                           taCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGGAAATTTCAACGTTCAGATTCCGCATTTc                       >  1:105767/1‑150 (MQ=255)
                                                                                                                             cTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGCCACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGt                     >  1:1552546/1‑150 (MQ=255)
                                                                                                                             cTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGt                     >  2:804015/1‑150 (MQ=255)
                                                                                                                             cTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGAAGCCGGCAATCAGGACACCAGGTACGATAACAAATCTACCGTCGTTAAATATCAAAGTTCAGCGTCCGCAGTTctt                     >  2:1203911/1‑148 (MQ=255)
                                                                                                                             attAAAGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGt                     <  2:974568/149‑1 (MQ=255)
                                                                                                                                aaCGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAAc                  <  1:982844/150‑1 (MQ=255)
                                                                                                                                      aaaaaaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTc            <  1:705680/150‑1 (MQ=255)
                                                                                                                                       aaaaaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCa           >  2:632619/1‑150 (MQ=255)
                                                                                                                                       aaaaaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCa           >  1:1743380/1‑150 (MQ=255)
                                                                                                                                       aaaaaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTTAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCa           >  1:907425/1‑150 (MQ=255)
                                                                                                                                       aaaaaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCa           >  1:628394/1‑150 (MQ=255)
                                                                                                                                        aaaaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTccc          >  1:309451/1‑148 (MQ=255)
                                                                                                                                          aaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTg        >  2:1018704/1‑150 (MQ=255)
                                                                                                                                                gaggaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATg  <  2:618848/150‑1 (MQ=255)
                                                                                                                                                gaggaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATg  <  1:1211984/150‑1 (MQ=255)
                                                                                                                                                gaggaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATg  >  1:1735193/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCAACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATG  >  NC_000913/3250842‑3251135

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