Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,926,3660AC23.3% 169.6 / 10.1 73G30G (GGT→GGGmioCflavoprotein MioC
Reads supporting (aligned to +/- strand):  ref base A (42/14);  new base C (0/17);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 1.50e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAA  >  NC_000913/3926227‑3926498
                                                                                                                                            |                                                                                                                                     
tCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTaaaa                                                                                                                                       >  2:380249/1‑140 (MQ=255)
    ttCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCg                                                                                                                                  <  1:1092046/140‑1 (MQ=255)
      ccTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAACCCCGcc                                                                                                                                <  2:433037/140‑1 (MQ=255)
            aTGCTTCATAGAAAGGAGAAAGGTTGCCCGAAATATCTCCGGCACCGGGGGTGGGGCTGATAACCAGCCAGACCCCTGAGGCAGGAAAATCTTCTAACAGCGGACCGGGCAGCGTTTCGGTGGTAAAACCCGCCTCTTcc                                                                                                                          <  2:885278/140‑1 (MQ=255)
             tGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCa                                                                                                                         <  2:564256/140‑1 (MQ=255)
               cTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGc                                                                                                                       >  1:409809/1‑140 (MQ=255)
                   aTAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAACCCCGCCTCTTCCAGCtttt                                                                                                                   <  2:557110/140‑1 (MQ=255)
                     agaaaggagaaagtTTTCCCGGATTACCCCGGGCGGGGTGGGTGGAGCGGTAAACCAGCGAGATCCCGGAGGCAGTTAAATCTTTAAAAGGGGCACGGGGACGCTTTTCGGTGTAAAAACCCGCCTCTTCCAGCTTTTCa                                                                                                                 <  2:486494/140‑1 (MQ=255)
                        aaggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAAAAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGcc                                                                                                              <  1:819493/140‑1 (MQ=255)
                         aggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCa                                                                                                             >  2:81558/1‑140 (MQ=255)
                              aaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGgtgt                                                                                                        >  1:1057098/1‑140 (MQ=255)
                               aaggTTGCCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGtt                                                                                                       <  2:75524/140‑1 (MQ=255)
                                aggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTc                                                                                                      <  2:994959/140‑1 (MQ=255)
                                 ggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCt                                                                                                     <  2:1057098/140‑1 (MQ=255)
                                     gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTa                                                                                                 >  2:602842/1‑140 (MQ=255)
                                     gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTa                                                                                                 >  2:1182655/1‑140 (MQ=255)
                                      tCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTAc                                                                                                >  2:728289/1‑140 (MQ=255)
                                      cccGGAATATCTCCGGCACCGTGGGTGGAGCTGATACCCAGCAAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGCCCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTAc                                                                                                <  1:1040161/139‑1 (MQ=255)
                                           aaTATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATAtt                                                                                           >  1:593315/1‑140 (MQ=255)
                                            atatCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAg                                                                                                            >  1:668853/1‑122 (MQ=255)
                                            atatCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAg                                                                                                            <  2:668853/122‑1 (MQ=255)
                                             tatCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCg                                                                                         >  2:767771/1‑140 (MQ=255)
                                             tatCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCg                                                                                         >  2:767010/1‑140 (MQ=255)
                                                ctcCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCa                                                                                      >  1:225422/1‑140 (MQ=255)
                                                 cccGGCCCGGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGCCCGGGCGGCGTTTCGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCAc                                                                                     <  2:409809/139‑1 (MQ=255)
                                                  ccGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCAcc                                                                                    >  2:831216/1‑140 (MQ=255)
                                                   cGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCAccg                                                                                   >  1:676342/1‑140 (MQ=255)
                                                     gCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCAccgcc                                                                                 >  1:419927/1‑140 (MQ=255)
                                                       aCCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGa                                                                               <  1:130592/140‑1 (MQ=255)
                                                           tgggtggAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAgggg                                                                           >  2:329097/1‑139 (MQ=255)
                                                           tgggtggAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAgggg                                                                           >  1:1214448/1‑139 (MQ=255)
                                                           tgggtggAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGt                                                                           >  1:979545/1‑140 (MQ=255)
                                                           tgggtggAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGt                                                                           >  2:701366/1‑140 (MQ=255)
                                                             ggtggAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGtgc                                                                         >  1:872272/1‑140 (MQ=255)
                                                             ggtggAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGtgc                                                                         >  1:393457/1‑140 (MQ=255)
                                                                 gAGCTGAAAACCAGCCAGACCCCTGGGGCAGGTAAATCTTCAAACAGCGGCCGGGGCAGCTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGcc                                                                     <  2:88067/140‑1 (MQ=255)
                                                                  ggCTGATACCCAGCCAGACCCCTGGGGAGGGTAAATCTTCTAACAGCGGCCGGTGCGGCGTTTCGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCg                                                                    <  2:177327/139‑1 (MQ=255)
                                                                   gCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGc                                                                   >  1:26152/1‑140 (MQ=255)
                                                                       aTACCCAGCCAGACCCCTGAGGCAGGAAAATCTTCTAACAGCGGCCGGGGCAGCTTTTCGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGAt                                                               <  2:953817/140‑1 (MQ=255)
                                                                       aTAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGAt                                                               >  1:771793/1‑140 (MQ=255)
                                                                       aTAACCAGCCAGATCCCTGAGGCAGGAAAATCTTCAAACAGCGGCCCGTGAAGGTTTTCGGGGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGAt                                                               <  2:706605/140‑1 (MQ=255)
                                                                        aacccagccagACCCTTGAGGCAGGAAAATTTTCTACCAGCGGCCGGTGCGGCTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATa                                                              <  2:225422/137‑1 (MQ=255)
                                                                        aacccagccaaATCCCGGGGGCAGTTAATTTTTCTAACGGCGGCCCGTGAAGCTTTCCGGGGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATa                                                              <  2:979545/137‑1 (MQ=255)
                                                                         acccacccagTCCCGTGAGGCAGTAAATTTTTAAAACAGCGGACGGTGCGGCGTTTGGGTGAAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATaa                                                             <  2:419927/138‑1 (MQ=255)
                                                                              gccagATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTg                                                        >  1:953817/1‑140 (MQ=255)
                                                                                  gaccccGGAGGCGGGTAAATTTTCTAACAGGGGCCGGGGCGGCTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGatat                                                    <  2:872272/137‑1 (MQ=255)
                                                                                   aTCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATc                                                   >  2:479212/1‑140 (MQ=255)
                                                                                   aTCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATc                                                   >  1:377304/1‑140 (MQ=255)
                                                                                             cAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAAc                                         >  1:1188856/1‑140 (MQ=255)
                                                                                             cAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAAc                                         >  1:316252/1‑140 (MQ=255)
                                                                                               ggTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCg                                       >  2:314425/1‑140 (MQ=255)
                                                                                               ggTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCg                                       >  1:125630/1‑140 (MQ=255)
                                                                                                gaaattcttcaACCAGCGGCCGGTGAAGGTTTTCGGTGATAAAGCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTTTTCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                      <  2:875873/135‑1 (MQ=255)
                                                                                                gTAAATCTTCTAACAGCGGCCCGTGCAGCGTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                      <  1:730269/140‑1 (MQ=255)
                                                                                                gTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                      >  1:353152/1‑140 (MQ=255)
                                                                                                gTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                      >  1:360697/1‑140 (MQ=255)
                                                                                                gTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                      >  2:724721/1‑140 (MQ=255)
                                                                                                   attctcttAACGGCGACCCGTGCGGCTTTCGGGTGTAAAACCCCGCCTCTTCCAGCTTTTCACCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCAc                                   <  2:417153/133‑1 (MQ=255)
                                                                                                        tctAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGCGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTtta                              >  1:875873/1‑140 (MQ=255)
                                                                                                          tAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                      >  2:1175285/1‑130 (MQ=255)
                                                                                                          tAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                      <  1:1175285/130‑1 (MQ=255)
                                                                                                          tAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTttatt                            >  2:1179317/1‑140 (MQ=255)
                                                                                                              aGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAaga                        >  2:509371/1‑140 (MQ=255)
                                                                                                                 ggACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTg                     >  2:719064/1‑140 (MQ=255)
                                                                                                                 ggACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTg                     >  1:520725/1‑140 (MQ=255)
                                                                                                                 ggACCGTGCAGCGTTTCGGTGGAAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTg                     <  1:831216/140‑1 (MQ=255)
                                                                                                                    ccggTGCAGCTTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGGGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGc                   <  1:565180/137‑1 (MQ=255)
                                                                                                                    cccGTGCAGCTTTTGGGTGTTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGc                   <  1:440329/139‑1 (MQ=255)
                                                                                                                       gTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTa                >  1:16895/1‑140 (MQ=255)
                                                                                                                            gCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTa           <  1:701366/140‑1 (MQ=255)
                                                                                                                            gCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTa           >  1:1128899/1‑140 (MQ=255)
                                                                                                                             tttttGGGTGTAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTAc          <  2:771793/138‑1 (MQ=255)
                                                                                                                             cGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTAc          >  1:210992/1‑140 (MQ=255)
                                                                                                                                  cGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCtgt     >  2:473803/1‑140 (MQ=255)
                                                                                                                                    gtggtAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGa   <  1:1179317/140‑1 (MQ=255)
                                                                                                                                    gtggtAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGa   <  1:509371/140‑1 (MQ=255)
                                                                                                                                    gtggaAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGa   <  2:26152/140‑1 (MQ=255)
                                                                                                                                     tggtAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGaa  <  1:767010/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                     
TCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAA  >  NC_000913/3926227‑3926498

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