Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,250,992 G→C 7.7% intergenic (+38/‑32) 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,9920GC 7.7% 243.7 / 10.8 89intergenic (+38/‑32)yqjF/yqjGputative quinol oxidase subunit/putative S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (53/28);  new base C (4/3);  total (57/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.38e-01

AACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGGC  >  NC_000913/3250844‑3251137
                                                                                                                                                    |                                                                                                                                                 
aaCTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAg                                                                                                                                                   >  2:384944/1‑149 (MQ=255)
aaCTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAg                                                                                                                                                   <  2:594783/149‑1 (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                                                 >  2:384379/1‑149 (MQ=255)
                                                                                                  tAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAt                                                 >  1:1257358/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)
                                                                                                                                                    |                                                                                                                                                 
AACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGGC  >  NC_000913/3250844‑3251137

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