Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,177,7480CA29.3% 38.6 / 13.8 41E55* (GAG→TAG) nudFADP‑sugar pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (25/4);  new base A (1/11);  total (26/15)
Fisher's exact test for biased strand distribution p-value = 4.55e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTA  >  NC_000913/3177510‑3177989
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
gggCTTGCCAGGAAACTTAACACCGGTTCGGTCCGTTTAACTATCACTCTCGCCTCTTCAATGGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGCCTGCGGCGTGACC‑GCGCTCAAAAATTTCccg                                                                                                                                                                                                                                       >  1:120048/1‑251 (MQ=255)
 ggCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCccgc                                                                                                                                                                                                                                      >  1:373365/1‑251 (MQ=255)
                            ggccccgttgaaCTACCAGCCCCGCCTCTTCAACCGCTCCGCGCCGGCCAACATCTTCCACCCTTCCCCCCTCTTCAACCACCCCGGCACCCACCCCCAGTAGCAAAGGGTTTTCGTGGGTGTCGTCCGCGGCAACCGGAATCTGCTCAATCAGCACAACTTCATCACGCACGGGGTCAAAGGGAAGCAAGCCTGGGGCTTGCCC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGtt                                                                                                                                                                                                            <  1:507921/239‑1 (MQ=255)
                                      gcctaccgcccccGCCTCTTCAACCGCTCCGCGACGGGCAACATCTCCCACCCTTCCCCCCTCTCCAATCACCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCGGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa                                                                                                                                                                                                  <  1:98932/240‑1 (MQ=255)
                                                 cctcccccacaaccccttcccgccGGGCACCTTCTTCCACCCTTCCCCCCTCTCAAATCACCCCGGCACCCATCCCCAGTAGCCAGGGGGTTTCGCTGGTGTCGGACGCGGCAACCCGAATCTGCCCAATAAGCCCAACTTCATCCCGCACTGGGTCAAAGGTAAGCAAGACGGCGGCGTGCCC‑GCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCta                                                                                                                                                                                      <  2:24460/228‑1 (MQ=255)
                                                   gCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCACTCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGACTAGACGATGACGAAATCtatc                                                                                                                                                                                    >  2:269151/1‑250 (MQ=255)
                                                         tCAATCGCTTCGCGCCGGGCAACATCTCCCACACTTCCCCCCTCTTCAATCACCCCGGCATCCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTCCGCGGCAACCCGAATCTGCTCAATCAGCACAACTTCACCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCC‑GCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGAt                                                                                                                                                                               <  1:411438/250‑1 (MQ=255)
                                                           aaTCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCt                                                                                                                                                                             >  1:117357/1‑250 (MQ=255)
                                                              cgctccccgACGGGAAACATCTCCCGCACTTCCCCCCTCTTCAACCAACCCGGCAACAATCGCAAGTAGCCAAGGGGTTTCGCTGGTTTCGGACGCGGCAATCCGAATCTGCTAAACCAGCACAACTTCATCACGCACTGGTTCAAAGGGTAGAAAGCCTGCGGCGTGGCC‑GGTCCAAAAAATTTACCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATg                                                                                                                                                                         <  1:269151/244‑1 (MQ=255)
                                                                cTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCACTCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGa                                                                                                                                                                        >  2:363838/1‑250 (MQ=255)
                                                                   cgcgACGGGCAACATCTTCCACACTTTCACCCTCTTCCATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCCATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCACGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGAATCATTTGCCCGTTGAATAGACGATGACGCAATCTATAAAGATCTAATGaaaa                                                                                                                                                                     >  1:281951/1‑250 (MQ=255)
                                                                      gACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaa                                                                                                                                                                   >  2:559013/1‑249 (MQ=255)
                                                                         ggccAACATCTTCCCCACTTCCCCCCTTTTCAATCACCCCGGCACCCATCTCCAGTAGCCAGGGGGTTTCGCTGGTGCCGTACGCGGCAACCCAAATCTGCTCAAGCACCACACCTCCAGCACGCATGGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcg                                                                                                                                                              <  1:363838/248‑1 (MQ=255)
                                                                           gCAACATATTCCACACTTTCCCCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGCAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                                                                                                                            >  1:564456/1‑251 (MQ=255)
                                                                            cAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                                                                                                                            >  1:231161/1‑250 (MQ=255)
                                                                             aaCATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcgct                                                                                                                                                          >  2:169105/1‑249 (MQ=255)
                                                                                  cTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCACGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGCCGAAATCTATAAAGATCTAATGACAAAAAGCCGCGATAAAgt                                                                                                                                                     >  1:631858/1‑251 (MQ=255)
                                                                                         acTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCAAAGGGGATTCGCTGGTGTCGTCCGCGGCAATCCGAATCTGCTCAATCAGCACACCTTCATCCCGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACA‑GCGCTCAAAAATTTCCCGCAGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGCAATATATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtctc                                                                                                                                              >  1:670829/1‑251 (MQ=255)
                                                                                                    ctacaaccaccccgccacccATCCCCAGTAGCCAAGGGGTTTCGCTGGGGTCGTACGCGGCAACCCGAACCTGCTCAATCAGCACAACTCCACCACGCACTGGGTAAAGGGGTAGCAAGCCGGCGGCGTGACC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAATTGTTTCTCGTGCaataa                                                                                                                                     <  1:559013/231‑1 (MQ=255)
                                                                                                      tcaaccaccccgccacccATCTCCAGTAGCCAGGGGGTTTCGGTGGTGTCGTACGCGCCAACCCAATTCGGCCCAATCAGCACACCTTCACCACGCTGGGGGCCAAAGGATAGCAGAACGGCGGCGTGACC‑GCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAAttt                                                                                                                                  <  2:51073/234‑1 (MQ=255)
                                                                                                         atcataCCGGCAACCCTATCCAATAACCAAGGGGTATAGTTGGTGTCGTCCGCGGCAATCCGAATATGCTCAATAAGCACAACTTCTTCACGCACTGGGTCAAAGGGAAGCAAGACTGCGGAGTGACC‑GCGCTCAAAAATTTACCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAGCTCTATAAAGATCTAATGAAAAAAAGCCGCGATAACGTGTTTCTCGTGCAATACTTTCTAc                                                                                                                              >  2:626018/1‑251 (MQ=255)
                                                                                                              cccGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGAGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGt                                                                                                                         >  1:405169/1‑251 (MQ=255)
                                                                                                                 ggCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGAACGCGGAAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCACGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGtttt                                                                                                                      >  1:365439/1‑251 (MQ=255)
                                                                                                                    aaCCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGGGTTTATCGTGCAATAATTTCTACATCGTTTTTGc                                                                                                                   >  2:383087/1‑251 (MQ=255)
                                                                                                                                aaGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGCCTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACgg                                                                                                       <  1:241059/250‑1 (MQ=255)
                                                                                                                                     aaGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGAAGCAAGACTGCGGCGTGAAC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAACCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAgg                                                                                                  >  1:598147/1‑251 (MQ=255)
                                                                                                                                                   gtgtCGGACGCGGCAACCCTAATCTGCACAAACCACAAAACTACATCACGCATGGGGAAAAAGGGTAAAAAACCGGGGGCGTGACC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAg                                                                                    <  2:373365/251‑1 (MQ=255)
                                                                                                                                                         tACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGt                                                                              >  1:524492/1‑251 (MQ=255)
                                                                                                                                                         tACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGt                                                                              >  2:526428/1‑251 (MQ=255)
                                                                                                                                                             cgGCAATCCGAATCTGCTCAATCCGCACAACTTCATCCCGCACTGGGACAACGGGTAGCAAGACTGCGGCGTGACC‑GAGCTCAAAACTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAACCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTCCATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAAt                                                                          >  2:462980/1‑251 (MQ=255)
                                                                                                                                                                cAATCCGAAGCTGCTCAACCACCAAAACTCCGCCACGCGGGGGGTCAAAGGGTAGCAAGACGGCGCCG‑GACCCGCGCTAAAAAGTTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTc                                                                       <  1:671646/251‑1 (MQ=255)
                                                                                                                                                                      gAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGc                                                                  >  2:567629/1‑250 (MQ=255)
                                                                                                                                                                                    aGCACAACTTCATCACGCACTGGGTCAACGGTTAGCACGACTGCGTCGCGACC‑GCGCTCAAAAATTTCCCGCCGCACCTAATGAATCATTTGCCCGTTGAATAGAAGATGACGAAATCTATCAAGGTCTAATGAAAAAAAACCGCGATAAAGTGTTTCTCGCACAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTACTAGTGGAtt                                                    >  2:674034/1‑250 (MQ=255)
                                                                                                                                                                                              catcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTCGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAAtt                                         >  1:43545/1‑251 (MQ=255)
                                                                                                                                                                                               atcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTAAAAAATTTCCCGCCGCACCCCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATACTTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAAttt                                        >  1:35814/1‑251 (MQ=255)
                                                                                                                                                                                                 caCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTc                                       >  2:218268/1‑250 (MQ=255)
                                                                                                                                                                                                      aCTGGGTAAAAGGGAAGCAAGCCTGCGGCGTGCCC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGAc                                 <  1:68779/251‑1 (MQ=255)
                                                                                                                                                                                                                aGGGTAGCAAGACTGCGGCTTGACC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGc                       <  1:295896/251‑1 (MQ=255)
                                                                                                                                                                                                                     aGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCCACTCATGACTCATTTGCCCGTTGAATAGACGATGACGAACTCTATAAAGATCTAATGAAAAAAAGCCGCGATTAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTTGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCCGCGACGTTTGACTGCCGttt                  >  2:603086/1‑251 (MQ=255)
                                                                                                                                                                                                                           cccgcggcgTCACC‑GCGCTCAAAAATCTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTGTTTGCCAAATGTAACGGGCAGGTTGTCGGGCTTAATCATTTTGAATGTCCTTGCACGAATAGGGAATTTAAGGTGAATTTCAGCGACGTTTGTCTTCCGTTTGAGCAg            <  2:354954/248‑1 (MQ=255)
                                                                                                                                                                                                                                     caCC‑GCTCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGAGCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGGTCTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTa  <  1:58675/250‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTA  >  NC_000913/3177510‑3177989

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