Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,177,7480CA26.8% 65.7 / 10.3 41E55* (GAG→TAG) nudFADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (28/2);  new base A (0/11);  total (28/13)
Fisher's exact test for biased strand distribution p-value = 2.47e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCA  >  NC_000913/3177528‑3177982
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
aaCACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCCAAGGGAAGCAAAACTGCGGCGTG‑AACGCGCTCCAAAATTT‑Cccgcc                                                                                                                                                                                                                              >  2:629870/1‑235 (MQ=255)
             tCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTc                                                                                                                                                                                                                >  1:87164/1‑236 (MQ=255)
                                         aaTCCCTGCGCGACGGTCACCATCTTCCACCCTTGCCCCCCCTCCAATCACCCCGGCAACCATCTCCAGTAGCCAAGGGTTTTCGCTGGTGTCGTACGCGGCACTCCGAATCTGCTCAATCAACACAACTTCACCCCGCACTGGGTCAAGGGGTAGCAAGACTGCGGCGTG‑CCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCg                                                                                                                                                                                                       <  2:728279/217‑1 (MQ=255)
                                                tCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATa                                                                                                                                                                                                >  1:728279/1‑217 (MQ=255)
                                                                   ccACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACCCATTTGCCCGTTGAATAGACGATGACGGAATCTATAAAGA‑CCAATGAAAAAAAga                                                                                                                                                         >  2:561447/1‑235 (MQ=255)
                                                                   ccACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCCCATGACTCATTTGCCCGTTGAATCGACGATGACGAAATCTATAAAGATCTAATGTAAAAAAg                                                                                                                                                          >  2:542020/1‑236 (MQ=255)
                                                                                   ttCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCCCAAACATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATTACGAAATGTCTAAGAATCTATTGCAAAAAATCCCCGATTAAGTGatt                                                                                                                                          >  2:175131/1‑236 (MQ=255)
                                                                                           tcCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                                                                                                             >  2:464712/1‑193 (MQ=255)
                                                                                           tcCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                                                                                                             <  1:464712/193‑1 (MQ=255)
                                                                                                       tctcCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACa                                                                                                                      >  1:720956/1‑236 (MQ=255)
                                                                                                                gCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGGAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGCCGATGACGAAATATATACAGATCTACTGCAAAACAGCCTCCTTAAAGACTGTTTCGTGCAACAATTTCTCCAttgactttg                                                                                                             >  2:619313/1‑228 (MQ=255)
                                                                                                                cccAAGGGGTTTTGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACCCATTTGCCCTTTGAATAAACGATGACGAAATCTATAAAGATCTAATGAAAAAAAACCTCGATAAAGGGTTTCTCCTGCAATAATTTCTACATCGTTATTg                                                                                                             >  2:256951/2‑236 (MQ=255)
                                                                                                                     gggtTTCCGCTGGAGTCAGCCGCGACAATCCACAGCTGCACAACCACCCCACCTTCATCACGCATGGGACCAAAGGGTAGCAAGACTGCGGCGGC‑ACCGCGCTCAAAAATTC‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCaa                                                                                                         <  2:87164/235‑1 (MQ=255)
                                                                                                                     ggggTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCaaa                                                                                                        >  2:104726/1‑236 (MQ=255)
                                                                                                                                      gaaCGCGGAAATCGAAACTTGATCAATCAGCACACCTCCTTCACGAACTGGGTCAAGGGGTAGAAGGACTGCGGCTTG‑CCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGt                                                                                       <  2:21741/234‑1 (MQ=255)
                                                                                                                                           cgacAATCGGAGTGTGCTAACTCAGCACACCTGCATCACGAACGGCGTAAAAGTGTAGCAAGACTGGGGCGTC‑ACCGCTCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGc                                                                                  <  2:520790/233‑1 (MQ=255)
                                                                                                                                           cgGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGc                                                                                  >  2:38979/1‑236 (MQ=255)
                                                                                                                                            ggCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGc                                                                                  >  1:21741/1‑235 (MQ=255)
                                                                                                                                            ggCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAACCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCCAATGTAACTGGCAGGTTGttctgct                                                                                 >  2:640653/1‑230 (MQ=255)
                                                                                                                                                  ccaaatgcgcacAATCAGAACATCTCCATCCCGCCTGGGGCCAAGGGGTACCAAGCCTCCGGCGTG‑CCCGCGCTAAAAATTTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCa                                                                           <  2:108421/225‑1 (MQ=255)
                                                                                                                                                        cTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGGTTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCAtt                                                                         >  2:637490/1‑232 (MQ=255)
                                                                                                                                                           cTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTaa                                                                    >  1:24126/1‑234 (MQ=255)
                                                                                                                                                             aaatccaccaaaCTTCATCCCGCACTGGTCCAAAGGGTAGCAAGACTGCGGCGTA‑ACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTc                                                                <  1:629870/226‑1 (MQ=255)
                                                                                                                                                                 cAGCACGATTTCATCCCGCACGGGGTCAAGGGGTAGAAAGACTGCGG‑GTGAACCGCGTTAAAAAATTC‑CCCGACGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTgg                                                            <  1:256951/236‑1 (MQ=255)
                                                                                                                                                                     acaACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑CCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGc                                                                                                            <  2:82666/184‑1 (MQ=255)
                                                                                                                                                                     acaACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGc                                                                                                            >  1:82666/1‑184 (MQ=255)
                                                                                                                                                                         cTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATa                                                    >  2:486245/1‑236 (MQ=255)
                                                                                                                                                                         cTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATATATACATCGTTTTTGCCAAATGTAACGGGCAGTTTGTCTGGCTTAAGCATTGTTaa                                                                    >  2:650327/1‑220 (MQ=255)
                                                                                                                                                                          ctcagcccccACTGGTTAAATGGGAAGCAAGCTGGCGTCTTG‑ACCCGCCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAg                                                   <  1:394885/227‑1 (MQ=255)
                                                                                                                                                                            caccacgcacgGGGTCAAAGGGTACCAAGACTGCGGCGTG‑CCCGCGCTAAAAATTTC‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGt                                                                                       <  2:307432/195‑1 (MQ=255)
                                                                                                                                                                                   cACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTAAATCGTTTTTGCCAAATGTAACGGGCAGGTTGGCTGGCtt                                                                                >  1:307432/1‑198 (MQ=255)
                                                                                                                                                                                           cAAAGGGTAACAATACTGCGGCTTG‑ACCGCGCTAAAAAATTC‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGt                                                                                       <  2:421206/183‑1 (MQ=255)
                                                                                                                                                                                                  tAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCtt                                                                                >  1:421206/1‑183 (MQ=255)
                                                                                                                                                                                                  tAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGa                           >  1:716143/1‑236 (MQ=255)
                                                                                                                                                                                                    gaaatgtggcggctTG‑ACCGCG‑TAAAAAATTTCCCCGCCGCACATCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGa                           <  2:720956/226‑1 (MQ=255)
                                                                                                                                                                                                     cAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGt                        >  1:292208/1‑236 (MQ=255)
                                                                                                                                                                                                          cTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGt                        >  1:315037/1‑231 (MQ=255)
                                                                                                                                                                                                          cTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCCGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGCATTAAATGTGAATTTTAGCGACGt                        >  2:580942/1‑231 (MQ=255)
                                                                                                                                                                                                           tGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACt                  >  1:254316/1‑236 (MQ=255)
                                                                                                                                                                                                             cggcgTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGc                >  2:499108/1‑236 (MQ=255)
                                                                                                                                                                                                                            tCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCa  >  1:487232/1‑236 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
AACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTG‑ACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCA  >  NC_000913/3177528‑3177982

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