Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,9530AC43.5% 111.2 / 23.7 70V661G (GTG→GGG) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (32/7);  new base C (0/30);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 3.31e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC  >  NC_000913/298803‑299103
                                                                                                                                                           |                                                                                                                                                      
ccGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCa                                                                                                                                                        <  1:993554/151‑1 (MQ=255)
     ggCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTc                                                                                                                                                    <  2:512726/150‑1 (MQ=255)
       cGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCAt                                                                                                                                                  >  1:1035391/1‑150 (MQ=255)
       cGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATa                                                                                                                                                 <  2:70515/151‑1 (MQ=255)
        gTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATAc                                                                                                                                                >  1:1294600/1‑151 (MQ=255)
            taCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCg                                                                                                                                            >  1:61483/1‑151 (MQ=255)
             acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGc                                                                                                                                           >  1:1366328/1‑151 (MQ=255)
                ccGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCCCCGACCCC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCACCCAGAAAAATCCCCTCCTTTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCat                                                                                                                                        <  2:765384/151‑1 (MQ=255)
                      tGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTg                                                                                                                                  >  2:409893/1‑151 (MQ=255)
                       ggCGATCGCCGCGCTCACGCCGCACGGCCCCAGCTCACCGACCCC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCACCCAGAAAAACCCCCCCCTTTTTG‑GGG‑ATATCCGCATGACCCGGCACCTCATCCCCCGCCATATCGt                                                                                                                                   <  2:1199028/149‑1 (MQ=255)
                               ccGCCCTCACGCCGAACGGCCCCAGCCCCCCGACCCC‑TTT‑GGCCTTCTCCGGGGGGGATTTGGGGC‑CGGTATCCTCCAGGAAAACCCCCTCGTTTTTT‑GTG‑ATACCCCCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGaa                                                                                                                         <  1:1239363/151‑1 (MQ=255)
                                cgcgCTCCCGCCGCACAGGCCCAGCCCCCCGACCCT‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGAAAAACCCCCTCCTTTTTG‑GGG‑ATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACg                                                                                                                           <  2:1282879/148‑1 (MQ=255)
                                   gcTTACCCCGCCCCGGCCCAGC‑CACCCTCCCCTTTT‑GGCCTTAACGGGGGGGGATATGGGTT‑CGTTATCAACCAGAAAAATCCCCTCTTTTTTT‑GGG‑ATACCCGCATGAACCGGCACCTAATCCCCCGCCATATCGTGATTAACGAAGTa                                                                                                                      <  1:210583/150‑1 (MQ=255)
                                      ccctccgcccggccccAGCTCCCCGACCCC‑TTG‑GTCCTTCAGCGGGGAGTAT‑TGGGGTAGGGTATAATCCAGGAAATTCCCCTCCTGTTTG‑GGG‑ATTTCCCCGTAACCCGGCTCCTCACCCCCCGCCATATCGTGATTAAc                                                                                                                            <  1:800100/128‑1 (MQ=255)
                                      cACGCCGCACAGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                  >  1:1335656/1‑151 (MQ=255)
                                        cgccgcACGGGCCCAGCTCACCGACCCT‑TTT‑GGCCTTCAGGGGGGAGGATAGGGGTT‑CGGTATCATCCAGAAAAACCCCCCCCTGTTTG‑GGG‑ATACCCGAATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaa                                                                                                                <  2:336451/151‑1 (MQ=255)
                                        cgccgcACAGGCCCAGCTCCCCGACCCC‑TTT‑GGCCTTCATCGGGGAGGATAGGGGGT‑CGGTATCATCCAGGAAAATCACCCCCTGTTTT‑GGG‑ATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa                                                                                                                 <  1:458696/150‑1 (MQ=255)
                                        cgccgcACAGGCCCAGCCCCCCGACCCC‑TTT‑GGCCTTCACGGGGGAGGATATGGGTT‑CGGTATCACCCAGAAAAATCACCCCCTTTTTG‑GGA‑ATACCCGCATAAACCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa                                                                                                                 <  2:362929/150‑1 (MQ=255)
                                        ccccgcACAGGCCCAGCTCCCCGACACC‑TTT‑GCCTTTCATCGGGGGGGATAGGGGGT‑CGGTATCCCCCAGAAAAACCCCCCCCTTTTTG‑GGG‑ATACCCGCATCACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGccct                                                                                                                 <  2:226225/148‑2 (MQ=255)
                                          ccgcACAGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                              >  1:302495/1‑151 (MQ=255)
                                          ccgcACAGGCCCAGCTCACCAACCCT‑TTT‑GGCTTTCATCGGGGAGGATATGGGGT‑CGTTATCACCCAGAAAAATCCCCTCTTTTTTT‑GGG‑ATACCCGCATGACCCGCCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                              <  2:1212741/151‑1 (MQ=255)
                                           cgcACAGGCCCAGCCCCCCGACCCC‑TTT‑GGCCTTCATCGGGGAGGATAGGGGGT‑CGGTATCATCCAGAAAAACCCCCCCCTGTTTT‑GGG‑ATATCCGCATGACCCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg                                                                                                             <  1:904326/151‑1 (MQ=255)
                                             cacaGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGt                                                                                                           >  1:1408878/1‑151 (MQ=255)
                                               cggccccAGCCCCCCGACACC‑TTT‑GGCTTTCATGGGGGAGGATAGGGGGT‑CGGTATCATCCAGAAAAACCCCCTCCTGTTTT‑GGG‑ATACCCGCATGAACCGCCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCaaaag                                                                                                             <  1:402819/143‑3 (MQ=255)
                                                aGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt                                                                                                        >  2:705650/1‑151 (MQ=255)
                                                 ggCCCAGCTCACCGACCCC‑TTT‑GGCCTTCAGCGGGGAGGATAGGGGGT‑CGGTATCACCCAGAAAAACCCCCTCCTGTTTT‑GGG‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                          <  2:649977/148‑1 (MQ=255)
                                                 ggCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATc                                                                                                       >  2:996069/1‑151 (MQ=255)
                                                  gCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc                                                                                                      >  2:1208990/1‑151 (MQ=255)
                                                    ccAGCCCACCGCCGCC‑TTT‑GGCCTTCAGGGGGGGGAATGGGGGTT‑TGGTTTCTCCAAGAAAAATCCCTCCCTTTTTT‑GGG‑ATACCCCTATGAACCGGCCTCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc                                                                                                      <  1:15302/149‑1 (MQ=255)
                                                        ccccccGACCCC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGTTATCATCAAGAAAAACCCCCTCCTTTTTT‑GGA‑ATACCCGCATGAACCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCg                                                                                                  <  2:242284/145‑1 (MQ=255)
                                                          cccccgcccc‑TTTGGGCCTTCATGGGGGAGGATGTGGGGT‑GGTTATCACCCAGAAAAATCCCCCCTTTTTTT‑GGA‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa                                                                                               <  1:1358314/143‑1 (MQ=255)
                                                                  cc‑TTT‑GCCTTTCATGGGGGAGGATTGGGGGT‑CGTTATCACCCAGAAAAACCCCCTCCTTTTTG‑GGG‑ATACCCGCATGAACCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcct                                                                                         <  1:416780/148‑1 (MQ=255)
                                                                      tt‑GCCCTTAATGGGGGGGAATGGGGGGT‑CGTTATCACCAAGAAAAACCCCCTCTTTTTTG‑GGG‑ATTCCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                                    <  2:1423891/150‑1 (MQ=255)
                                                                        tgcccTTCATGGGGGAGAATATGGGGT‑GGGTATCACCCAGAAAAATCCCCTCCTTTTTT‑GGG‑ATACCCGCAGAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCa                                                                                   <  1:77317/147‑1 (MQ=255)
                                                                           ccTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc                                                                               >  2:1218913/1‑151 (MQ=255)
                                                                             ttCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg                                                                             >  1:1303427/1‑151 (MQ=255)
                                                                               caggggggAGGATAGGGGGT‑CGGTACCA‑CCAAGAAAAACCCCCCCTTTTTTGGGG‑ATACCCGCATGACCCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGcc                                                                              <  2:1391577/144‑1 (MQ=255)
                                                                                   ggggAGGATATGGGGT‑CGTTATCACCCAGAAAAACCCCCTCCTGTTTT‑GGG‑ATATCCGCATGACCCGGCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccca                                                                       <  1:1236092/151‑1 (MQ=255)
                                                                                    gggAGGATAGGGGGT‑CGTTATCACCCAGAAAAACCCCCTCTTTTTTT‑GGA‑ATACCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                      <  1:99869/151‑1 (MQ=255)
                                                                                       aggaTATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcc                                                                   >  2:845458/1‑151 (MQ=255)
                                                                                           tatGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                >  2:875071/1‑150 (MQ=255)
                                                                                              ggggT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                            >  1:362932/1‑151 (MQ=255)
                                                                                              ggggC‑CGGTACCATCAAGAAAACTC‑CCTCCTTTTTT‑GGGAATCCCCCCAACACCCGCCCCCTTTCCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                            <  1:46964/151‑1 (MQ=255)
                                                                                                 gT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCAtt                                                         >  1:1217646/1‑151 (MQ=255)
                                                                                                     ggTATCACCCAGGAAAATCCCCTCCTTTTTT‑GGA‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgc                                                       <  1:347614/150‑1 (MQ=255)
                                                                                                        atcatcCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc                                                    >  1:1349403/1‑150 (MQ=255)
                                                                                                          catcCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                                 >  1:666627/1‑151 (MQ=255)
                                                                                                             ccAGAAAAATCCCCTCCTTTTTG‑GGA‑ATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg                                                <  1:1192428/149‑1 (MQ=255)
                                                                                                               aGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                            >  1:716873/1‑151 (MQ=255)
                                                                                                                gaaaaaTCCCCTCTTTTTTT‑GGG‑ATCCCCGCAGAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg                                                <  1:416711/144‑1 (MQ=255)
                                                                                                                  aaacccccctcttTTTTT‑GGG‑ATACCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGa                                               <  2:1274067/138‑1 (MQ=255)
                                                                                                                  aaaaTCCCCTCTTGTTTT‑GGA‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGc                                         <  1:409950/151‑1 (MQ=255)
                                                                                                                   aaacccccccccgTTTT‑GGGAATACCCGCAGAACCCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg                                                <  2:1469600/132‑1 (MQ=255)
                                                                                                                        aCCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgc                                    >  2:545054/1‑150 (MQ=255)
                                                                                                                         cctcttTTTTG‑GGG‑ATACCCGCATGACCCGCCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgtgc                                    <  1:798527/149‑1 (MQ=255)
                                                                                                                            cctGTTTG‑GGG‑ATACCCCCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCa                                 <  1:830755/149‑1 (MQ=255)
                                                                                                                              tGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                             >  2:1152811/1‑151 (MQ=255)
                                                                                                                              tGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                             >  2:72581/1‑151 (MQ=255)
                                                                                                                                ttgt‑GGG‑ATACCCGAATGACCCGCCCCCTAACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtt                              <  2:26017/145‑1 (MQ=255)
                                                                                                                                 ttt‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCg                          >  2:270525/1‑151 (MQ=255)
                                                                                                                                  tt‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                         >  2:1082861/1‑151 (MQ=255)
                                                                                                                                  tt‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                         >  1:385817/1‑151 (MQ=255)
                                                                                                                                   t‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                         >  2:646558/1‑150 (MQ=255)
                                                                                                                                      ggaataccccaaTGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCg                          <  2:1366225/137‑1 (MQ=255)
                                                                                                                                         atatCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc                     >  1:330757/1‑150 (MQ=255)
                                                                                                                                          gccccgcgggacccGGCCCCAAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc                     <  2:1035291/137‑1 (MQ=255)
                                                                                                                                            tCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                  >  2:1354167/1‑150 (MQ=255)
                                                                                                                                            tCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                  >  1:810606/1‑150 (MQ=255)
                                                                                                                                             ccGCATGACCCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGAt                <  2:1335549/151‑1 (MQ=255)
                                                                                                                                               gCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc              >  1:866963/1‑151 (MQ=255)
                                                                                                                                                       cGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc       >  2:1434655/1‑150 (MQ=255)
                                                                                                                                                           aCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac  >  1:92925/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                                      
CCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACC‑TTT‑GGCCTTCATCGGGGAGGATATGGGGT‑CGGTATCATCCAGGAAAATCACCTCCTGTTTT‑GGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC  >  NC_000913/298803‑299103

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑