Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,760,3840AC30.3% 120.3 / 19.8 76V289G (GTG→GGG) selAselenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base A (45/8);  new base C (2/21);  total (47/29)
Fisher's exact test for biased strand distribution p-value = 2.76e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTC  >  NC_000913/3760147‑3760605
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                             
gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCAc                                                                                                                                                                                                                              <  2:1098495/239‑1 (MQ=255)
     gATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCCAACAACTTGTCGCCGGGGAAACTCAcccg                                                                                                                                                                                                                           >  1:1505722/1‑235 (MQ=255)
                ccTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGCTCATCTCTTTTTTTCCAACAATAAATCCTGCCTTCGGGCCCCCCAAAAACTTTTCGCCCGAGAAAATCACCAGaacg                                                                                                                                                                                                                       >  2:336733/1‑227 (MQ=255)
                    tgcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTTCGGCCCGCCTAAAAAATTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCaa                                                                                                                                                                                                          >  1:470117/1‑239 (MQ=255)
                    tgcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGGGAAACTCACCAGACTGAcgccccccccca                                                                                                                                                                                                          >  1:1231002/1‑235 (MQ=255)
                    tgcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCaa                                                                                                                                                                                                          >  1:763955/1‑239 (MQ=255)
                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGCTCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATTAACt                                                                                                                                                                                                    >  1:248528/1‑238 (MQ=255)
                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTTCCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACCCACCAGACTGACGCCCGCCGCAATCAAcc                                                                                                                                                                                                    >  2:621415/1‑237 (MQ=255)
                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACCGACCCCCGCCGCAATCAACt                                                                                                                                                                                                    >  1:1775088/1‑238 (MQ=255)
                                                         taaCTCGGGGCTTCAGGGTGTAAAGAGAGAGGAAACGTGGATGCCCGGCGCGGGAGTTTCGTTTCACCCGCGCATATTCCCGCCGAGAGCGGTTGGCTTGCCGGGGGGGCAATCATCTTTTTTTTCCAAACAATAATCCCGGCCGGGGCCGCGCCAAACACTTTGTCGCGGGAAAACCCCCCCAGACTGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  2:1774089/237‑1 (MQ=255)
                                                         ccACTCCGACGTTCGGGGTGTAAAAAAAGCCGAAGCGTGGCCTCCGGCGCGCGGGGGGTCTTTTTCCCCGCGCGTAGTGCACCCTTAAGCGGGTGGCTGGGCGGGCGGGCAACCTTCTCTTTTTTACAAACAATATTGCCGCCCGGCGGCCCGCCAAACATCTTGTCGGCGGAAAACCCCCCCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  2:248528/238‑1 (MQ=255)
                                                                 agCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAAAACTTCCTGCCTGCGGCCCGCCCAACAACCTGTCGCCCGAGAAAACCACCAGACCGCGGCCCGCCGCAATCAACCCCCCCCGGCCGCGCTCCTTCGGGAAACCGTGCTGGCtca                                                                                                                                                             >  2:1361053/1‑237 (MQ=255)
                                                                   cTTCAGGGTGAAAAAAAAGACGCAGCGTGGCTTCCGGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCGGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTGGTCGCGGGAAAACCTCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCa                                                                                                                                                                            <  2:661324/222‑1 (MQ=255)
                                                                             taaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTTCCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAAATCACCAGACCGACCCCCGCCCCAATCAACTCCCGCGGGCTTGGCCCTTTTGGGAAACCCTCCTGGCTAAAATCCACCAgtc                                                                                                                                                 >  2:157296/1‑237 (MQ=255)
                                                                             taaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTg                                                                                                                                                                  >  1:661324/1‑222 (MQ=255)
                                                                                  aaaGCCGCAACGTGGCTTCCAGCGCCGCGGGGGTCATTTTACCCGCGCGTAATGCACGCTTCAGGGGGTGGCTTTGCGGGCGGGCGATCATCTCTTTTTTCCCAAAAATAATCCCTCCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                            <  1:1613221/239‑1 (MQ=255)
                                                                                   aaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATACCTTTTTTAACAACAATAATTCCTGCCTGCGGCCCCCCCAACCACCTTTTGCAGGAGAAACTCCCCACACTGCCACCCGCCGCAATCAACTCCt                                                                                                                                                                                                 >  1:216216/1‑185 (MQ=255)
                                                                                     gACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCCATCAACTCCTGCGGCATTGGCCCTTTCGGCCAACCGTACTGGCTAAGATCGACCAGCGAGcc                                                                                                                                             >  1:1452460/1‑235 (MQ=255)
                                                                                         cAACGTGGCTTCCAGCGCCGCGAGGGGCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCTTCTCTTTTTTTCCAACAATAATTCCTGCCTGCGGCCCGCCCAAAAACCTGTCGCCCGAGAAACTCACCaca                                                                                                                                                                                                                          >  2:613302/1‑152 (MQ=255)
                                                                                                 cTTCCCGCGCCGAGGGGGTCATTTGCTCCGAGGGTTAAGCACCATGAGGCGTGTGGCTTTGGAGGCGGGCATCCCTCTTTTTTTAACAAAAAATATTCCCCGCCGGCGGCCCGCATAACAACTTGTCGCGGGAAAAACCCACCAGACGACCCCCCGCCGCAATCAACTCCTGCGGCATGGGCTCTTTCGGCAAACCGTACTGGCAAAGATCGCCCAGCGAGCCACCGCCTAAATCAGt                                                                                                                              <  2:402997/238‑1 (MQ=255)
                                                                                                    ccAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTaaaa                                                                                                                                  >  1:1608043/1‑230 (MQ=255)
                                                                                                          ccgcGAGGGTCATTTTCTCCGCGTGTAAGCCCCCCTTCGGCGGGTGGCTTTGCGGGCGGCCGATCATCTTTTTTTTACCAAAAATAATCCCTCCCGGCGGCCCGCCTAACAACTTGCCGCCGGAAAAACTCACCAGCTTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                    <  1:1241039/239‑1 (MQ=255)
                                                                                                                               cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCCAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCCGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTAc                                                                                                                        >  1:191007/1‑214 (MQ=255)
                                                                                                                                  gTAATGCGCGTTTCAGGGGGTGGCTTTGGAGGGGGACAATCATCTTTTTTTTCCAAAAAAAACTCCCTGCCAGGGGCCCGCCTAACATCTTGTCGCCGGAAAAACCCCCCAGCCGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  2:1494922/166‑1 (MQ=255)
                                                                                                                                      tGCCCGCTACAGGGGGTGGTTTTGCAGGCGGGCGATCATCTCTTTTTTACCAAAAAAAATTCCTGCCGGGGGCCCGCCAAACAACTTGGCGCGGGAGAAACCCCCCAGACTGACGCCCGCCGCAATCAAc                                                                                                                                                                                                     <  1:1334220/130‑1 (MQ=255)
                                                                                                                                          cGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCCACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCcgcg                                                                                                                                                                                              >  1:1236022/1‑133 (MQ=255)
                                                                                                                                          cGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCCCCCGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCCCTTTCGGCAAACCGTACTGGCCAAAATCGACCAGCGAGCCACTCCCTAAATCAGTCACTACGGGAACATCCCGCTCCTTTCCGAGCCcca                                                                                           >  2:580337/1‑232 (MQ=255)
                                                                                                                                            cTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACCATAATTCCTGCCTGCGGCCCGCGTAACAACTTGTCCCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGCACTGGCTaaaa                                                                                                                                                           >  1:1788448/1‑164 (MQ=255)
                                                                                                                                            cTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCCAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACCCCTGCGGCATTGGCCCTTTCGGCCAACCGTACCGGCTAAGa                                                                                                                                                           >  1:1494922/1‑166 (MQ=255)
                                                                                                                                             ttCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGGCCGCCTAACAACTTTTTGCCGGGGAAACTCACCCGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAAATCGACCAGCGAGCCCCCGCCTAAATCAGTCACCACGGGAAAATCCCGCTCTTTGCCGCGCGccacccagtccc                                                                                  >  2:507052/1‑231 (MQ=255)
                                                                                                                                             ttCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTcc                                                                                   >  1:590481/1‑237 (MQ=255)
                                                                                                                                              tCAGCGGGGGGCTTTGGCGGGGGGGCAATATTCTTTTTTTACCCAACAATATCCCCTCCCGGCGCCCCGCCAAACATCTTGCCGCGGGAAAACCTCCCCAGACTGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAAc                                                                                                                  <  2:1739540/206‑1 (MQ=255)
                                                                                                                                                aGCGGGTGGCTTTGCAGGCGGGGGATCATTTCTTTTTTACCAACAAAAAATCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGGGAAAATCACCAGACTGACGCCCGCCCCCAACAAATCCTGCGGCATTGGCTCTTTTGGCAAACCGTACTGGCTAAGATCGACCAACGAGCCCCCGCCTCAATCAACCACACTGGGAAAATCCCGATCTTTGTCCAGCGCCCCCACTTTCGcct                                                                               >  2:422147/1‑236 (MQ=255)
                                                                                                                                                aGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTTCCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCTCCGGAGAAACTCACCAGACTGACCCCCGCCGCCATCAACTCCTGCgtcc                                                                                                                                                                                           >  2:1334220/1‑127 (MQ=255)
                                                                                                                                                aGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAAAAATAATTCCTGCCTTCGGCCCCCCCAACAACTTGTTGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGTCATTGGCTCTTTCGGCAAACCGGACTGGCTAAAATCGACCAGCGAGTCCCTCGCCAAATCCGTCACTACCGGAAAATCCCGCTCTTTTCCGAGCGCCACCCGTTgcgccc                                                                               >  2:1376591/1‑236 (MQ=255)
                                                                                                                                                    ggTGCCTTGGCAGGCGGGCGTTCATCTTTTTTTTCCAAACAAAAATCCCTGCCTGCGGCCCGCCTAACAACTTGCCGCCGGAGAAACTCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                            <  1:624759/173‑1 (MQ=255)
                                                                                                                                                      tGGTTTTGCAGGGGGGCGATCATCTTTTTTTCACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAAAAACCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTa                                                                         <  1:507052/238‑1 (MQ=255)
                                                                                                                                                       ggCTTTGCAGGTGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGctct                                                                                                        >  1:1739540/1‑206 (MQ=255)
                                                                                                                                                       ggCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGGACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTAAGGGAAAATCCAGCTCTTTGCCCAGCGCCACCAGTTCCGc                                                                                 >  1:83186/1‑229 (MQ=255)
                                                                                                                                                               cAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCCGTTCCCCCTCATCTATCGCTTTgg                                                                >  1:1584227/1‑239 (MQ=255)
                                                                                                                                                               cAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCCTACTGTCTAAGATCGACCAGCGAGCCCCTGGCTAAAt                                                                                                                                  >  2:624759/1‑173 (MQ=255)
                                                                                                                                                                  gcgggccACCATCTCTTTTACAAAAACAATTACCCCGGCGTGCCGCCCGCCAAACATCTTGCCCGCGGAAAACCCCCCCAGCCAGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                              <  2:1452460/238‑1 (MQ=255)
                                                                                                                                                                    gggccATCATCTCTTTTTGACCAAAAATAATCCCGGCCGGCGGCCCGCCTAACAATTTGCCCCCGGAGAAACCCACCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                            <  1:1361053/238‑1 (MQ=255)
                                                                                                                                                                      gcgATCATCTTTTTTTCCCCAAAAAAAATCCCCGCCGGCGCCCCCCCAACCATCTTGTCCCAGGAAAAACCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                            <  2:1738377/156‑1 (MQ=255)
                                                                                                                                                                        gATCATCTTTTTTTCCCCAAAAAAACTTCCTGCCGGCGCCCCGCCTAACACCTTGTCGCCGGAGAAACTCCCCAGACTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                       <  1:157296/239‑1 (MQ=255)
                                                                                                                                                                         atcatcTTTTTTACACAAACAAATACCCCCGCATGCGCCCCGCCTAAAACTTTGCCCCGGGAAAACCCCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                       <  2:1775088/238‑1 (MQ=255)
                                                                                                                                                                                ctTTTTTACCATCAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAAt                                                                                                                                  >  1:1738377/1‑156 (MQ=255)
                                                                                                                                                                                ctTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                    >  1:603229/1‑170 (MQ=255)
                                                                                                                                                                                 ttttttCCCCACAAAAATCCCGCGCCCGCCGCCCGCCAAACAATTTGCCGCGGGAGAACCCCCCCACCCGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTg                                               <  2:763955/239‑1 (MQ=255)
                                                                                                                                                                                 tttttcACCAACAAAAATCCCGGCCGGGGCCCCGCCTAACAACTTGTCGCCGGAAAACCTCCCAAGCCGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGt                                              <  1:580337/239‑1 (MQ=255)
                                                                                                                                                                                         aaaaaaaaatCCCCGCCGGCGGCCCGCCTAACACCTTGCCCCCGGAAAACCCCACCGGACTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTg                                      <  1:422147/232‑1 (MQ=255)
                                                                                                                                                                                           acaaTAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     >  2:538955/1‑110 (MQ=255)
                                                                                                                                                                                           acaaTAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  1:538955/110‑1 (MQ=255)
                                                                                                                                                                                           acaaTAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                    >  1:90376/1‑239 (MQ=255)
                                                                                                                                                                                             aataatCCCTGCCTGCGCCCCCCCTAACAACTTGTCGCCGGAGAACCCCCCCAGCCTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  2:432735/108‑1 (MQ=255)
                                                                                                                                                                                                aatCCCTGCCTGCGGCCCGCCTAACAACTTGCCGCCGGAAAACCTCCCCAGACTGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                               <  1:621415/239‑1 (MQ=255)
                                                                                                                                                                                                        gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTg                                                                                                                                                                                                >  2:951117/1‑71 (MQ=255)
                                                                                                                                                                                                        gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGa                                                                                                                                                           >  1:432735/1‑108 (MQ=255)
                                                                                                                                                                                                        gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTTCCTAAATCAGTCACTACGGTAACATCCAGCTATTTTCCGAGTGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGTACCCCTGCATGCAGTAGTTACTCGTATGTACTTtca                         >  2:401892/1‑238 (MQ=255)
                                                                                                                                                                                                        gcctgcGGCCCGCCAAACAACTTGTCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTg                                                                                                                                                                                                <  1:951117/71‑1 (MQ=255)
                                                                                                                                                                                                         cctgcGCCCCGCCAAACACCTTGCCGCGGGAAAACCTCCCCAGACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGc                                                                                                                                                                             <  2:1561575/89‑1 (MQ=255)
                                                                                                                                                                                                          cggcgcCCCCCCAAAAAACTTTCCGCGGGAAAACCCCCCCAGCCTGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCaaa                                                                                                                                                                          <  2:380945/89‑1 (MQ=255)
                                                                                                                                                                                                           tgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGt                                                                                                                              <  1:1525784/134‑1 (MQ=255)
                                                                                                                                                                                                           tgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGt                                                                                                                              >  2:1525784/1‑134 (MQ=255)
                                                                                                                                                                                                             cGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  1:1157951/93‑1 (MQ=255)
                                                                                                                                                                                                             cGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     >  2:1157951/1‑93 (MQ=255)
                                                                                                                                                                                                             cGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTTTCGCTTTGGTg                                                              >  1:9625/1‑196 (MQ=255)
                                                                                                                                                                                                             cGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACa                   >  2:777811/1‑239 (MQ=255)
                                                                                                                                                                                                                   gCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACt                                                                                                                                                                   >  1:1561575/1‑89 (MQ=255)
                                                                                                                                                                                                                    ccTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGc                                                                                                                                                                >  1:380945/1‑91 (MQ=255)
                                                                                                                                                                                                                     cTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTa           >  2:1437472/1‑239 (MQ=255)
                                                                                                                                                                                                                              tGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGt                                                                                      >  2:436301/1‑155 (MQ=255)
                                                                                                                                                                                                                              tGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGt                                                                                      <  1:436301/155‑1 (MQ=255)
                                                                                                                                                                                                                              tGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTa                                                                         >  1:158395/1‑168 (MQ=255)
                                                                                                                                                                                                                               gTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc  >  2:72217/1‑238 (MQ=255)
                                                                                                                                                                                                                                    cGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATg                                 >  1:645736/1‑202 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                             
GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTC  >  NC_000913/3760147‑3760605

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