Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,250,989 G→C 7.9% 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 7.9% 216.5 / 10.3 87intergenic (+35/‑35)yqjF/yqjGputative quinol oxidase subunit/putative S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (51/27);  new base C (4/3);  total (55/32)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

GTCAACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGGC  >  NC_000913/3250841‑3251137
                                                                                                                                                    |                                                                                                                                                    
gTCAACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGTTAAAAAACGaggagg                                                                                                                                                      >  2:317541/1‑149 (MQ=255)
gTCAACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGaggagg                                                                                                                                                      >  1:1545119/1‑149 (MQ=255)
   aaCTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAg                                                                                                                                                   <  2:594783/149‑1 (MQ=255)
   aaCTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAg                                                                                                                                                   >  2:384944/1‑149 (MQ=255)
     cTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGct                                                                                                                                                 >  2:1107825/1‑149 (MQ=255)
           gATGTTCATGAAACACCTGACAATTTCTTGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGt                                                                                                                                           <  2:167796/149‑1 (MQ=255)
             tGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGttt                                                                                                                                         <  1:1300096/149‑1 (MQ=255)
             tGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGGATAGCATCGACCGCCTGCTGAATAAAAAGTGTTAAGCACGCGCTATACTTAACGATAAAAACCGAGGAGGAAGCTCCTCGttt                                                                                                                                         <  2:1428864/149‑1 (MQ=255)
              gTTCATGAAAAACCTGACAATTTCTGGCGGATTCCGGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGtttt                                                                                                                                        >  1:1502931/1‑149 (MQ=255)
               ttCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCCCGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTGCTCGttttt                                                                                                                                       >  2:374744/1‑149 (MQ=255)
                 cATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGc                                                                                                                                     >  2:1136639/1‑149 (MQ=255)
                  aTGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATATCATAGACCTCATGCTGAATAAAAAGTTGTAATCACTCTCTATACTTAACGATAAAAAACGAGGAGGAAGATACTCGttttttat                                                                                                                                    >  2:1377321/1‑146 (MQ=255)
                  aTGAAAAACCTGACAATTTCTGGCGGATTCCTACTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCt                                                                                                                                    >  2:523721/1‑149 (MQ=255)
                  aTGAAAAACCTGACAATATCTGGCGGATTCCTGCTGCTTGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAACAGTGGTAAGCACGCGCTATACATAACGATAAAAAACGCGGAGCTTCCTCCTCGttttttat                                                                                                                                    <  1:626861/149‑4 (MQ=255)
                    gAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCTCCTGCTGAATAAAAAGTGGTAAGAACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAt                                                                                                                                  >  1:878591/1‑149 (MQ=255)
                    gAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATAGTTAACGATAAAAAACGAGTAGGAAGCTCCTCGTTTTTGCTAt                                                                                                                                  >  2:119323/1‑149 (MQ=255)
                    gAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAt                                                                                                                                  >  1:603245/1‑149 (MQ=255)
                     aaaaaCCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAtt                                                                                                                                 >  2:445433/1‑149 (MQ=255)
                      aaaaCCTGCAAATTTATGGCGGATTCCTGCTGCTGGCAATTTCCGGGCCGGGCGCGTCTAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTAAACTATAAAAAACGATGTGGAAGCTCCTCGTTTTTGCTATTg                                                                                                                                <  2:1051634/149‑1 (MQ=255)
                      aaaaCCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGACGCTCCTCGTTTTTGCTATTg                                                                                                                                >  1:499203/1‑149 (MQ=255)
                       aacccTGACAATATCTGGCGGATTCCTGCTGCTGGCAATTACCGGGCCTGGCGCGGATAGCATCGACCGCCTGCTGACTAAAAACTGGTAAGCACGCGCTAAACTTACCGCTAACCCACGCGGAGCAACCGCCTCGTTTTTGCTATTgg                                                                                                                               <  1:1235791/146‑1 (MQ=255)
                            tGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCCCGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGgaggag                                                                                                                          >  2:1453739/1‑149 (MQ=255)
                             gACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGATGAATAAAAAGTGGTAAGCACGAGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGgaga                                                                                                                         >  1:917820/1‑149 (MQ=255)
                              acacTTTCTGGCGGATTCCTGCTGCTGGCAATTACGGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCCCGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGgagag                                                                                                                        <  2:1418723/149‑1 (MQ=255)
                                 aTTTCTGGCGGATTCCTGCTTCTGGCAATTACAGGTCCGGGCGCGTATAGCATCTAGCGACTGCTGAATAAATAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGATGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGaaa                                                                                                                     <  1:982705/149‑1 (MQ=255)
                                     cTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATg                                                                                                                 <  2:1098496/149‑1 (MQ=255)
                                      tGGCGGTTTCCTGCTGCTGGGAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATgg                                                                                                                <  2:22912/149‑1 (MQ=255)
                                      tGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATgg                                                                                                                >  1:792859/1‑149 (MQ=255)
                                      tGGCGGATTCCTGCTGCTGGCAAATACCGGGCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGGAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATgg                                                                                                                <  2:1451163/149‑1 (MQ=255)
                                        gCGGATTCCTGCTGCTGGCAATTACCGGGCCGGGCGCGTATAGTATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGt                                                                                                              <  1:684407/149‑1 (MQ=255)
                                         cGGATTCCTGCTGCTGGCAATTACCGGGCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTc                                                                                                             <  1:1322390/149‑1 (MQ=255)
                                             ttCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACt                                                                                                         <  1:147289/149‑1 (MQ=255)
                                                    tgctgGCAATTACCGGTCCGGGCGCGTATAGCATCGCCCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGAc                                                                                                  >  1:577508/1‑149 (MQ=255)
                                                    tgctgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGAc                                                                                                  >  1:1232062/1‑149 (MQ=255)
                                                        ggCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTTTTGGAGGAGAGATAAATGGGTCAACTGATTGACGGCg                                                                                              >  1:921230/1‑149 (MQ=255)
                                                           aaTTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCt                                                                                           >  1:460295/1‑149 (MQ=255)
                                                           aaTTACCGGTCCGGGCGCGTATAGCATCGACCGCCGGCGGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGGTTTTGGGATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCt                                                                                           >  2:1099455/1‑149 (MQ=255)
                                                                ccGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGGGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCAt                                                                                      <  2:609904/149‑1 (MQ=255)
                                                                 cGGTACGTGCGCGTATAGCATGTACCGCCTGCTGAATAAATAGTGGTAACCACGCGATATACATAACGATAAAAACCGCGGAGTAAGCTCCTCGGTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATg                                                                                     <  2:1163727/149‑1 (MQ=255)
                                                                  ggTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGa                                                                                    >  2:164748/1‑149 (MQ=255)
                                                                      cGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACAcc                                                                                <  2:241470/149‑1 (MQ=255)
                                                                      cGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACAcc                                                                                >  2:454575/1‑149 (MQ=255)
                                                                              tataGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGCGGAGGAAGCTCCTCGTTTTTGCTATTGGTGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGa                                                                        >  1:966856/1‑149 (MQ=255)
                                                                              tataGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGa                                                                        >  1:1259239/1‑149 (MQ=255)
                                                                                   cATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCa                                                                   <  1:809991/149‑1 (MQ=255)
                                                                                    aTCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAACCGAGGAGGAAGCGCCTCGTTTTTGCTATTGGCGGAGAGAAAAATGGGTCAACTCATTGACGGCGTCTGGCATGACACCTGGTACGATACCaa                                                                  <  2:319328/149‑1 (MQ=255)
                                                                                     tCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCaaa                                                                 <  1:676976/149‑1 (MQ=255)
                                                                                        aCCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGATTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCt                                                              >  2:1454278/1‑149 (MQ=255)
                                                                                           gCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTAcc                                                           <  2:1502931/149‑1 (MQ=255)
                                                                                                  gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGGGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGAGATGACACCTGGTACGATACCAAATCTACcggcggga                                                    >  2:1055715/1‑147 (MQ=255)
                                                                                                  gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTa                                                    >  1:1204621/1‑149 (MQ=255)
                                                                                                  gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGACAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCTGCGGTa                                                    >  2:503724/1‑149 (MQ=255)
                                                                                                   aaTAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTAGCGGCGGTaa                                                   >  1:39030/1‑149 (MQ=255)
                                                                                                   aaTAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTaa                                                   >  2:890510/1‑149 (MQ=255)
                                                                                                   aaTAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGCAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTaa                                                   <  1:159891/149‑1 (MQ=255)
                                                                                                     tAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAt                                                 >  1:1257358/1‑149 (MQ=255)
                                                                                                     tAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAt                                                 >  2:384379/1‑149 (MQ=255)
                                                                                                      aaaaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACAGATTGACGGCGTCTGGCATGACACCTGGTACGATAACAAATCTACCGGCGGTAAAtt                                                >  1:57528/1‑149 (MQ=255)
                                                                                                        aaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTc                                              >  2:538468/1‑149 (MQ=255)
                                                                                                           gTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAAc                                           >  2:446690/1‑149 (MQ=255)
                                                                                                               tAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTc                                       >  1:774003/1‑149 (MQ=255)
                                                                                                                  gCAAGCGCTATACTTAACGATAAAAAACGAGGAGGAATCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGTGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGc                                    <  1:1377321/149‑1 (MQ=255)
                                                                                                                    aCGCGCTATACTTAACGATAAAAAACGATGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTAAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACAGGCGGTAAATTTAAACGTTCAGCtt                                  >  2:1267478/1‑149 (MQ=255)
                                                                                                                    aCGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCtt                                  >  2:958900/1‑149 (MQ=255)
                                                                                                                      gcgcTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTcc                                >  2:558542/1‑149 (MQ=255)
                                                                                                                          tataCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTAAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCAt                            <  2:878591/149‑1 (MQ=255)
                                                                                                                             aCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAGTTCAACGTTCAGCTTCCGCATTTc                         >  1:774671/1‑149 (MQ=255)
                                                                                                                              cTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCg                        >  1:852199/1‑149 (MQ=255)
                                                                                                                              cTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTCGAGGAGAGAAAAATGGGTCAACTGATTGACTGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCTGCGGTAAATTTCAACGTTCAGCTTCCGCAATTCg                        >  2:318681/1‑149 (MQ=255)
                                                                                                                               ttAACGATAAAAACCGAGGAGGAAGCTCCTCGTATTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCTTGACACCTGGTACGATACCAAATCTACCGCCGGTAAATTTAAACGATCAGCGTCCGCAATTCGt                       >  1:1343995/1‑149 (MQ=255)
                                                                                                                               ttAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAttt                                               <  2:356941/125‑1 (MQ=255)
                                                                                                                               ttAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAttt                                               >  1:356941/1‑125 (MQ=255)
                                                                                                                                tAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCTTa                      >  2:188653/1‑149 (MQ=255)
                                                                                                                                   ctaTAAAAAACGAGGAGTATGCTCCTCGTTTTTGGTATTGGAGGAGAGAAAAATGGGTCAACTTATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTGCAGCTTCCGCATTTCGTAACt                   <  2:1505565/147‑1 (MQ=255)
                                                                                                                                    gATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCCGCTTCCGCATTTCGTATCTg                  >  2:628907/1‑149 (MQ=255)
                                                                                                                                     aTAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTgg                 <  1:523721/149‑1 (MQ=255)
                                                                                                                                      tAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGc                <  1:119323/149‑1 (MQ=255)
                                                                                                                                       caaaaaCGAGGAGCTTCCTCCTCGTTTTTGCTATTTGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTAACGGCGGTAAATTTCAACGTTCAGCTTACGCATTTCGTAACTGGCt               >  2:1064650/2‑149 (MQ=255)
                                                                                                                                       caaaaaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAAGTGATTGACGGCGTCAGGCATGACACCTGGTACGAGACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGcg               >  2:502691/2‑148 (MQ=255)
                                                                                                                                       aaaaaaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCt               <  1:1136639/149‑1 (MQ=255)
                                                                                                                                       aaaaaaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCt               >  1:443936/1‑149 (MQ=255)
                                                                                                                                         aaaaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCa             >  1:575097/1‑149 (MQ=255)
                                                                                                                                         aaaaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCTTAACTGGCTCa             >  2:1569018/1‑149 (MQ=255)
                                                                                                                                          aaaCGCTGAGGAAGCGCCTCGTTTTTGCAATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCAc            <  2:961754/149‑1 (MQ=255)
                                                                                                                                          aaaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCAc            >  2:1118932/1‑149 (MQ=255)
                                                                                                                                          aaaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCAc            >  1:1376066/1‑149 (MQ=255)
                                                                                                                                           aaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGGCTGGCATTACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTGCCGCATTTCGTAACTGGCtaccg           >  2:1403213/1‑145 (MQ=255)
                                                                                                                                             cGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGc         <  1:384944/149‑1 (MQ=255)
                                                                                                                                              gaggagGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCACCGTTCAGATTCCGCATTTCGTAACTGGCTCACTTcc        >  1:519522/1‑149 (MQ=255)
                                                                                                                                                 gaggaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAGAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCAGAt     >  1:805275/1‑149 (MQ=255)
                                                                                                                                                 gacgaAGCTCCTCGTTTTTGGTATTGGAGGAGAGAAAAATTGGTCAACAGATTGACTGCGTCTGGCATGACACCAGGTACGATACCAAATCTACCGGCGGGAAATTTCAACGTACAGCTTCCGCATTTAGTAAATGGCTCACTTCCGAt     >  2:923642/4‑149 (MQ=255)
                                                                                                                                                   ggaAGCTCCTCGTTTTTGATATTGGAGGAGAGAAAAATGGGGCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATgg   <  1:1107825/149‑1 (MQ=255)
                                                                                                                                                    gaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATggt  >  2:187360/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCAACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGGC  >  NC_000913/3250841‑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: