Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,343,1940CT78.4% 62.4 / 15.8 38D428D (GAC→GAT) melAalpha‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base T (13/16);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAGGCGGGCGTTGGTAACAGCGACCCGATACCCTATGAGCATTTCAATGACTACAAAACTCAGTTATGGATTTGGAGCGTTCGGGAAGGATTTTGCGATCGGCATTGTGT  >  NC_000913/4342946‑4343444
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          
aaCGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACg                                                                                                                                                                                                                                                           >  1:129828/1‑250 (MQ=255)
          aCGATGGTTTGATTGATACCCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCAGTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGTACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGct                                                                                                                                                                                                                                                <  1:126584/251‑1 (MQ=255)
            gATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGctct                                                                                                                                                                                                                                              >  1:402869/1‑251 (MQ=255)
                  ttgatggaTAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCTCGATGATGGACCCTCATACTGCCGCCGTGCTGGGCATTGATGAAAGATATGCTCttgttg                                                                                                                                                                                                                                         <  2:281523/250‑1 (MQ=255)
                                     aaGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCAc                                                                                                                                                                                                                     >  1:390505/1‑251 (MQ=255)
                                           gTTGCGAGGAAGTATGCTGTCTTGTTGATGATAATGGCATTCAGCGGACCAACGTGGGTACGCAACCTTCGCATCAGGCCGCCCTGATGCAGACCAACATCAACGTACAGCCGCGGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGAGATGGTGGCCCCCCATACTGCCGCCGTGCTGGGCGTTTAGGAAATATATGATCTTGTTAACGACCTGCTTGCCTCCCTCGGCGAc                                                                                                                                                                                                               >  2:104732/1‑251 (MQ=255)
                                                    aaGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGCTGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGcc                                                                                                                                                                                                       >  2:31700/1‑250 (MQ=255)
                                                      gTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTCCGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCAAGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAgg                                                                                                                                                                                                    >  2:25698/1‑251 (MQ=255)
                                                             gTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTg                                                                                                                                                                                             <  1:310227/251‑1 (MQ=255)
                                                                 ggTTGATGCTACTGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATCGCTGACATAAATGCTCTTGTTGACGACCTGATTGCCGCCCCCGGCGACTGGCTGCCAGGCTGGTTGCAcc                                                                                                                                                                                         >  2:259065/1‑251 (MQ=255)
                                                                              ggCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGAc                                                                                                                                                                            >  1:142839/1‑251 (MQ=255)
                                                                               gCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACt                                                                                                                                                                           >  1:280386/1‑251 (MQ=255)
                                                                                                     tACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGAGACTGGCTGCCAAGCTGGTTGCACCGTTAAAACGCGACTAAACGCCACTGCGCCGGGGGAt                                                                                                                                                     >  2:223928/1‑251 (MQ=255)
                                                                                                            ccTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTcc                                                                                                                                              <  1:31700/251‑1 (MQ=255)
                                                                                                                         ccgcaCTGATGCAAACCAACATCAACGTACAGACTCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCTCATACTGCCGCCGTGCTGGGCATGGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGGTGCGCCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTCTTCCGGCGCACACctct                                                                                                                                 <  2:111654/251‑1 (MQ=255)
                                                                                                                          cgccCTGATGCAAACCAACATCAACTTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCTCACACCTCTg                                                                                                                                <  2:129828/251‑1 (MQ=255)
                                                                                                                            ccgatcggCAAACCAACATCAACGTACAGACGCTGCTGTCCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGAGTAAACGCTACGGCGCCGGGGGATTTATTCCGGCGCACACCTCTGAc                                                                                                                              <  1:43532/244‑1 (MQ=255)
                                                                                                                             ccTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACg                                                                                                                             >  1:236004/1‑251 (MQ=255)
                                                                                                                                   gCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCACAACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCAGCCGTGCTGGGCATTGATGCAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCAGGCGCACACCTCTgccagatagg                                                                                                                       >  2:43532/1‑242 (MQ=255)
                                                                                                                                             aTCAACGTACAGACGCTGCTGACCGAAGCCATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGGCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAg                                                                                                             <  1:145972/251‑1 (MQ=255)
                                                                                                                                                 ccGTACAGGCGCTGCGTACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGAGGATGGACCCGCATACGGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAggcgg                                                                                                         <  2:390505/250‑1 (MQ=255)
                                                                                                                                                  cGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAggcggg                                                                                                        <  1:160154/251‑1 (MQ=255)
                                                                                                                                                                        gCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCGGACGATACCAAGAACAGAAGGCGGGCGTTGGTAACAGCGACCCGATa                                                                                  <  1:210260/251‑1 (MQ=255)
                                                                                                                                                                               ttactgataaTCGCGACCGTGTTTCCCACGCCGCGATGATGGACTCGCATACCGCCGCCGTGCTGGTCATTGATGAAATATATGCTCTTGTGGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGTATTTATTCCGGCGCACTCCTCTGACGATACCTATACCAGAAGGCGGGCGTTGGTAACATCGTCCCGATACCCTAt                                                                            <  2:199936/242‑1 (MQ=255)
                                                                                                                                                                                  cGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAGGCGGGCGTTGGTAACAGCGACCCGATACCCTATGAGc                                                                        <  1:223928/251‑1 (MQ=255)
                                                                                                                                                                                    gAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATaa                                                                                                                   >  1:227399/1‑206 (MQ=255)
                                                                                                                                                                                    gAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACGGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATaa                                                                                                                   <  2:227399/206‑1 (MQ=255)
                                                                                                                                                                                                 gtgtTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGCCTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAGGCGGGCGTTGGTAACAGCGACCCGATACCCTATGAGCATTTCAATGACTAc                                                          >  1:248481/1‑250 (MQ=255)
                                                                                                                                                                                                           cgccgcGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAGGCGGGCGTTGGTAACAGCGACCCGATACCCTATGAGCATTTCAATGACTACAAAACTCAGtt                                               <  1:377257/251‑1 (MQ=255)
                                                                                                                                                                                                              cgcgATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACt                                                                                                                                                                 <  2:72244/134‑1 (MQ=255)
                                                                                                                                                                                                              cgcgATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACt                                                                                                                                                                 >  1:72244/1‑134 (MQ=255)
                                                                                                                                                                                                                    gatgGACCCGCATACTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCATGGCGGCTGTCTGTCAGTGTTGTTGCGCGGTTAAAACGCGATGAAACGCGTCGGCGCCTGGTGAGTTATTCCGGCGCACACCTCTGACGAGACCAATAACAGAAGGGGTGCGTTGGTAACAGCGGCCCGGGACCGTATGAGCATATCAATGACTACAAAACTCAGTTATGGAttt                                       <  1:404848/250‑1 (MQ=255)
                                                                                                                                                                                                                       ggcccagCATACTGCCGGCGTTCTGGGCATTGACGAAATATAGGCTCGTGTTGACGACCTGATTGCCGCCCACGGCGTCTGGCTGCCAGGCTGGTTGCAACGGTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAGGCGGGCCTTGGTAACAGCGACCCTATACCCTATGAGCATTTCAATGCCTACAAAACTCAGTTATGGATTTGGaa                                   <  2:122936/245‑2 (MQ=255)
                                                                                                                                                                                                                                 aCTGCCGCCGTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAGGCGGGCGTTGGTAACAGCGACCCGATACCCTATGAGCATTTCAATGACTACAAAACTCAGTTATGGATTTGGAGCGTTCGGGaa                         <  1:259065/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ccgccgTGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAGGCGGGCGTTGGTAACAGCGACCCGATACCCTATGAGCAGTTCAATGACTACAAAACTCAGTTATGGATTTGGAGCGTTCGGGAAGGAt                     <  2:280386/251‑1 (MQ=255)
                                                                                                                                                                                                                                           tGCTGGGCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAGGCGGGCGTTGGTAACAGCGACCCGATACCCTATGAGCATTTCAATGACTACAAAACTCAGTTATGGATTTGGAGCGTTCGGGAAGGATTTTGCg               >  1:248661/1‑251 (MQ=255)
                                                                                                                                                                                                                                             cggggCATTGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAGGCGGGCGTTGGTAACAGCGACCCGATACCCTATGAGCATTTCAATGACTACAAAACTCAGTTATGGATTTGGAGCGTTCGGGAAGGATTTTGCGAt             <  2:143665/249‑1 (MQ=255)
                                                                                                                                                                                                                                                ggCATGGATGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAGGCGGGCGTTGGTAACAGCGACCCGATACCCTATGAGCATTTCAATGACTACAAAACTCAGTTATGGATTTGGAGCGTTCGGGAAGGATTTTGCGATCgg          <  2:402869/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        cGAAATATATGCTCTTGTTTACGATCTGATTGCCGCCCTCGGCGACTGGCTGCCAGGCTGGTTGCATCGTCAAAACGCGACTAAACGCGTCTGCGCCTGGTGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAGGCTGGCGTTGGTAACAGCGACCCGATTCCCTATGAGCATTTCAATGACTACAAATCTCATTTATGGATTTGGAGCGTTCGGGAAGGATTTTGCGTTCGGCAgtgtgt  <  2:92612/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          
AACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGCATCTGGCCGCCCTGATGCAAACCAACATCAACGTACAGACGCTGCTGACCGAAGCTATTCTTACGGAAAATCGCGACCGTGTTTACCACGCCGCGATGATGGACCCGCATACTGCCGCCGTGCTGGGCATTGACGAAATATATGCTCTTGTTGACGACCTGATTGCCGCCCACGGCGACTGGCTGCCAGGCTGGTTGCACCGTTAAAACGCGACTAAACGCTACTGCGCCGGGGGATTTATTCCGGCGCACACCTCTGACGATACCAATAACAGAAGGCGGGCGTTGGTAACAGCGACCCGATACCCTATGAGCATTTCAATGACTACAAAACTCAGTTATGGATTTGGAGCGTTCGGGAAGGATTTTGCGATCGGCATTGTGT  >  NC_000913/4342946‑4343444

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