Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913954,3880CT31.3% 93.1 / 20.8 68A27T (GCA→ACA) focAformate channel
Reads supporting (aligned to +/- strand):  ref base C (42/4);  new base T (0/21);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 9.75e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATA  >  NC_000913/954155‑954600
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                    
gcgcCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTtat                                                                                                                                                                                                                 >  1:530652/1‑239 (MQ=255)
         cgacttgacaaaaacaaTCAACAGGGGGAAAGTAAAGAGATCGGCTCCGCAAACACCAAAAAGAACCAGCCCCAGAGAGAAGCAAATGCCGCCACCCATTTTTGCCAGGCCGAAGGCCATTGTGCCTGCGCCAGTGTTTGCT‑GTGATAT‑AGA‑AGACGAATGCAATTGAGATAAAACCC‑C‑CGGCG‑TTAA‑TCGCCAG‑AAAGAAAGTTTTAG‑GCGGTTGTTTCGTTGTTTTATAGACCCCCg                                                                                                                                                                                                        <  2:639039/231‑1 (MQ=255)
          gCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GGAA‑TCGCCAG‑ATAGAAAGACTTAA‑ACGGATGTTTCGTTGCTTTATAGACACCCg                                                                                                                                                                                                        >  1:1127124/1‑238 (MQ=255)
                  aacaacaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCCCTTCg                                                                                                                                                                                                >  1:2293276/1‑238 (MQ=255)
                      acaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCCGcc                                                                                                                                                                                             >  1:791151/1‑237 (MQ=255)
                      acaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCCCCCCTTTGGCCaa                                                                                                                                                                                           >  1:782548/1‑238 (MQ=255)
                            aacggggTGAAAGAAAAGAGACCGGCCCCGCAAAAAACAAAAAGAACCAGCCCCAGAGAGAAGCAAGTGCCGCCAACAGGTTTGGCCATGCCGAGGGGATTTGGGCCGGTGCCAGTGGTTGTG‑GAAATAA‑AAA‑AGAGGAAGGCGATTGAGAGAAAAACC‑C‑GGGGG‑GTAA‑CCCCCAG‑AGAGAAAGTTTAAA‑GGGTTTGTTTCGTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGc                                                                                                                                                                                     <  2:791151/234‑1 (MQ=255)
                                  gTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGGTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTgc                                                                                                                                                                               >  1:1181642/1‑239 (MQ=255)
                                                       ccGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGACGGTTCGTTGCTTTATAGACACCCCCCTCTTCGGCCACTTTGGCCATTGCCGCAGGAAGTAAAAGATCaaa                                                                                                                                                          >  1:618582/1‑239 (MQ=255)
                                                        cGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCCG‑ATAGAAAGGCTTAA‑GAGGATGGTTCGTTGCTTTATTGACACCCGCCTCTTCGGCCACTTTTGCCATTGCAGCAGGAAGTAAAAGATCaaaa                                                                                                                                                         >  1:1046412/1‑239 (MQ=255)
                                                        cGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTTCAGCAGGAAGTAAAAAATCaaaa                                                                                                                                                         >  1:829768/1‑239 (MQ=255)
                                                               aacaCAAAGACTCACCCCCAGAGAGAAGCAAGTGCCGCCAACCGTTTTTCCGATGCCGAGGGCAATTGTGCTTGTGCCGGTGTTTGCG‑GTAATAA‑AAA‑AGGCGAAGGCGATTGAGAAAAAACCC‑C‑CGGCG‑TAAT‑TCCCCAG‑AAAAAAAGTCTAAA‑GCGGTTGTTTTGTTGTTTAAAAGCCCCCCGCCTCTTCGGCCACTTTGGc                                                                                                                                                                                     <  2:1081561/204‑1 (MQ=255)
                                                               aacaCAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCCCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAAATCAAAAgg                                                                                                                                                       >  1:819328/1‑234 (MQ=255)
                                                                         aTCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTgcagca                                                                                                                                                                           >  1:1081561/1‑204 (MQ=255)
                                                                                ccAGGGAGCAGAAAATGCCACCAACCTGTTTGGCCAAGCCGAAGGGCTTTGTGCTGGGGCCAGTGGTGGGT‑GTAATAA‑AAA‑AGACGAATGCGATTGAGATAAAACCG‑C‑CGCCG‑GAAA‑TCCCCAG‑ATGAAAGGTATAAG‑GCGATTGTTTCTTTGTTTTATAAAAACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAAct                                                                                                                                 <  2:782548/239‑1 (MQ=255)
                                                                                        aaGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAAGGTTGTCAGCTTTCACACTa                                                                                                                                    >  1:246055/1‑228 (MQ=255)
                                                                                         aGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCCCTTCGGCCACTTTGGCCATTTCAGCAGGAAGTAAAAAATCAAAAGGGTTTTCAGCTTTCCCACCaaa                                                                                                                                  >  1:1987108/1‑228 (MQ=255)
                                                                                                     cAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGAAACCCGCCTCTTCGGCCCCTTTGGCCCTTGCATCAGGAAGTAAAAAATCAACAGGGGTGTCCGCCTTTACACTAACTCTCCCTTTATTAATTCGGcg                                                                                                             >  2:1989857/1‑238 (MQ=255)
                                                                                                      caCCAGTTTTGGCAAGCCGAAGGGCATTGTCTCTGT‑CCAGGGGGTGGTGGTGAAGT‑AGAGAGA‑GAAGTAGGTTG‑GA‑GTAAACA‑CGCGGGG‑TTAAGTCGCCAT‑ATAAATAGTTTTAG‑TGGTTTGTTTTGTTGTTTAATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCACTTTATTAAGGCGGcgac                                                                                                           <  2:1265641/238‑1 (MQ=255)
                                                                                                            ttttGCCATGCGGAGGGGCGTTTGGCTTGGGCGGGTGTTTGCT‑GAAATAA‑AAA‑AGACGAATGGGATTGAAAAAAAAACC‑C‑GGGGG‑AAAA‑TCGCAAG‑ATGAAAATTCTAAA‑GGGGTTTTTTGTTTGCTTTAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcgac                                                                                                           <  2:530652/233‑1 (MQ=255)
                                                                                                                   aTGCGGAGGGACATTGGGCTGGTGCAGGTGGTTGTT‑GTAAAAA‑AAA‑GGAGGAAGGCGATGGAGA‑AAAAACC‑C‑CGGGC‑GGAA‑ACCCCAGAAGAGAAAGTCTAAA‑GGGGATTTTTTGTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACaa                                                                                              <  2:829768/239‑1 (MQ=255)
                                                                                                                   aTGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAAGCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAAATCTCTCTTTATTAAGTCGGCGACGAGATACTAACaa                                                                                              >  1:2328779/1‑239 (MQ=255)
                                                                                                                       cgaggggttttGGGCCTGTGCCGGTGTTGTTT‑GAGATAT‑AAA‑AGACAAATGGTATTGAGAGAAAACCC‑C‑CGGCG‑TTAC‑GCCCCAG‑ATAAAAATTTTTAA‑GCGGGTTTTTCGTTGTTTTAGAGCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCa                                                                                          <  2:819328/230‑1 (MQ=255)
                                                                                                                                ttGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCCTTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAg                                                                                                                   >  1:1063574/1‑205 (MQ=255)
                                                                                                                                   tgcctgtgccAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTTTAGATgaga                                                                              >  1:815127/1‑239 (MQ=255)
                                                                                                                                   tgcctgggccGGTGTTTGTT‑GTAATAA‑AGA‑AGACAAAGGCGTTTGAGAAAAAACCC‑C‑GGGCG‑TAAA‑CCGCCAA‑AAAAAAATTTTAAA‑GGGGTTTTTTCTTTGTTTTAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctc                                                                                                                              <  2:2103586/191‑1 (MQ=255)
                                                                                                                                     cctgtgccGGTGGTTGGT‑GTG‑TATAAGA‑AGAAAAAGGATAGTGAGAGAAAAACC‑C‑CGGTG‑GAAA‑TC‑CCAGAAAAGAAGGTTTAAA‑GTGATTGTTTTGTTGTTTAAGAACCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGaa                                                                             <  2:2383201/238‑1 (MQ=255)
                                                                                                                                        gtgccAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTg                                                                         >  2:380082/1‑239 (MQ=255)
                                                                                                                                             aGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGGAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTaa                                                                                                                    >  1:2103586/1‑191 (MQ=255)
                                                                                                                                                   tGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTTTTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAAATCat                                                              >  1:1229869/1‑239 (MQ=255)
                                                                                                                                                      t‑gtaATAA‑AGA‑AGACAAATGCGATTGAGATAAAAACC‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAAAAAGTCTAAA‑GGGGTTTTTTTGTTTTTTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTgcagca                                                                                                                                                                           <  2:102285/127‑1 (MQ=255)
                                                                                                                                                          taaaaaaaa‑agaCAAAGGCATTTGAGAGAAAACCC‑C‑CGGCG‑GTAA‑CCGCCAG‑AAAAAAAGTCTAAA‑GGGGTTGTTTTTTTGTTTTATAGCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGatat                                                                     <  2:904755/219‑1 (MQ=255)
                                                                                                                                                            tat‑aGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTCACACTAACTCTCTCTTTATTTAGTCGGCGACGAGATACTaa                                                                                                 >  1:1679303/1‑196 (MQ=255)
                                                                                                                                                             taaaga‑agaCGAATGCGATTGAGATAAAAACC‑C‑CGGGG‑TTAA‑TCGCCAG‑ATAGAAATTCTTAG‑GCGGTTTTTTCGTTGTTTTAAAAACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctc                                                                                                                              <  1:2069763/164‑1 (MQ=255)
                                                                                                                                                              t‑aGA‑AGAAGAATGCGATGGAGATGAAAACC‑C‑CGGCG‑GTAA‑TCGCAAG‑AAAAAAAGTCTTAA‑GGGTTTGTTTCGTTGTTTTAAAGCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcgacga                                                                                                         <  2:313269/186‑1 (MQ=255)
                                                                                                                                                               aata‑agaCGAATGCGTTTGAGATGAAACCC‑C‑CGGCG‑TTAA‑TCGCCAG‑AAAAAAATTCTAAA‑GCGGTTTTTTCTTTGTTTTAAAAACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                               <  2:1369328/193‑1 (MQ=255)
                                                                                                                                                                  a‑agaCGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaaaa                                                                                                                                                                 >  1:102285/1‑127 (MQ=255)
                                                                                                                                                                     gaCGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATAtc                                                           >  1:904755/1‑227 (MQ=255)
                                                                                                                                                                        gAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGc                                                                                                                                                                                     >  1:449402/1‑102 (MQ=255)
                                                                                                                                                                        gAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCCGCAGGAAGTAAAAGAACAAAAGGGTTTTCAGCTTTCACACTAACTCTCTCTTTTTTaa                                                                                                                    >  2:2069763/1‑167 (MQ=255)
                                                                                                                                                                        gAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GAAA‑TTGGCGG‑AGAGAAAGGCGTAA‑GCGGGGGTTTCGTTGCTTTATAGTAAACCGCCCCTTCGGCCACTTTGGCCAATTCAGTATGAAGTAAAAGATAAAAAGGGGTGGCAGATTTTACACTAAATCTCTCTTTATTAAGTCGGCGTAGAGATACAAACAAAACATTATAGATGAGAAATTGATATAAATAATAATCCGCCTGGCTTAAAg                                            >  2:1097401/1‑239 (MQ=255)
                                                                                                                                                                        aagagtgcggtagaTGTAAACACC‑CGGGGTGTAA‑TCGCAAG‑ATAAAATGTGAAAA‑GCGGATGTTGCTTGGGTTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACAAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCGCGCCTGGCTtat                                              <  2:772820/228‑1 (MQ=255)
                                                                                                                                                                         aaTGAGATTGAGATGAAAACA‑C‑GGGCG‑GTAA‑GCGCCAG‑AAAAAATGTTTAAA‑GCGTTTGTTTCTTTGTTTTAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTtata                                             <  1:1052885/237‑1 (MQ=255)
                                                                                                                                                                          aagCGTTTGAGATAAAAACC‑C‑CGGCT‑GTAA‑CCG‑CAA‑TTAAAAAATTTTAAGGGGGTTTTTTCGTTTTTTTAAGAACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCgg                                                                                                               <  2:1364669/168‑1 (MQ=255)
                                                                                                                                                                          aTGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCCG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                               >  1:313269/1‑186 (MQ=255)
                                                                                                                                                                          aTGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATtata                                                                                     >  1:1369328/1‑196 (MQ=255)
                                                                                                                                                                          aTGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGACTTAA‑GCGGGTGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCCGGAAGTAAAAGATCAAAAGGGGTGTTAGCTTTCACACTAACTCTCTCTTTATTAAGTTGGCGGCGAGATACTAACAAAGCATTATAGATGATGAATTGATTTTGAACATATTTCGCCTGGCtttttgg                                           >  2:1302993/1‑233 (MQ=255)
                                                                                                                                                                                   agaTGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTTACACTCACTCTCTCTTTATTAAGTCGGCGACGAGATACTAAAAAAGCATTATATATGAGAAATTGATATAGATCATATATCGCCCGGCTTATAGGCCCGTAAc                                   >  2:1052885/1‑237 (MQ=255)
                                                                                                                                                                                   aaaaaaaaaCC‑C‑CGGCG‑GAAA‑CCCCCAG‑AAAGAAAGTCTTAA‑GCGGTTTTTTCGTTGTTTTAAAACCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcgac                                                                                                           <  2:373683/160‑1 (MQ=255)
                                                                                                                                                                                       gAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCCCTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCg                                >  1:1424042/1‑236 (MQ=255)
                                                                                                                                                                                        aaaaCA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTTTAGACAACCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACCAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGa                                                                        >  1:752346/1‑195 (MQ=255)
                                                                                                                                                                                        aaaaCA‑C‑CGGCG‑GTAA‑TCCCCAG‑ATAGAAAATCTTAA‑ACGGATGTTTCGTTGCTTtata                                                                                                                                                                                                                >  2:817723/1‑59 (MQ=255)
                                                                                                                                                                                          aaCA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGGCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCCCTTCGGCCACTTTGGc                                                                                                                                                                                     >  1:1840083/1‑84 (MQ=255)
                                                                                                                                                                                             a‑c‑CGGCG‑TTAA‑ACGCCAG‑ATAAAAAGTCTAAA‑GCGGGTTTTTCGTTGCTTAAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTttatt                                                                                                                      <  2:145223/144‑1 (MQ=255)
                                                                                                                                                                                             a‑c‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTaa                                                                                                 >  1:373683/1‑165 (MQ=255)
                                                                                                                                                                                                 cggcg‑gTAA‑TCGCCCG‑ATAGAAAGTCTTAA‑GCGGGTGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCCGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATTGATGAGAAATTGATATAGATCATATTCCGCCTGGCTTATAGGCCCGTAAATCGCATGGGTTTTa                     >  2:666432/1‑239 (MQ=255)
                                                                                                                                                                                                          a‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcga                                                                                                            >  1:145223/1‑144 (MQ=255)
                                                                                                                                                                                                               ccAG‑ATAGAAAGTCTTAA‑GCGGTTTTTTCTTTGCTTTAAACACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGa         <  1:1302993/239‑1 (MQ=255)
                                                                                                                                                                                                                cAG‑ATAGAAAGTCTAAG‑GGGGTTGTTTCGTTGTTTTATAGACACCCGCCTCTTCGGCCACTTTg                                                                                                                                                                                       <  1:2183571/64‑1 (MQ=255)
                                                                                                                                                                                                                  aaagaGAAATTCTAAA‑GAGGTTGTTTGGTTGTTTTAACAACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACGCTCTCTTTATTAAGTCGGCGAGGAGATGCTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGGGt       <  1:1201591/235‑1 (MQ=255)
                                                                                                                                                                                                                       aaaatttttAA‑GTGGTTGTTTCGTTGTTTAATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAata  <  2:246055/232‑1 (MQ=255)
                                                                                                                                                                                                                           gTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTgcag                                                                                                                                                                             >  2:2183571/1‑64 (MQ=255)
                                                                                                                                                                                                                            tCTTAG‑GCGGATGTTTCGTTGTTTTATACACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTg                                                    <  2:1593865/184‑1 (MQ=255)
                                                                                                                                                                                                                               tAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTcaca                                                                                                                                       >  1:1176191/1‑98 (MQ=255)
                                                                                                                                                                                                                                       aTGTTTCTTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTg                                                    <  1:1756959/174‑1 (MQ=255)
                                                                                                                                                                                                                                       aTGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGc                                                       >  1:1838303/1‑171 (MQ=255)
                                                                                                                                                                                                                                       aTGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTg                                                    >  2:1756959/1‑174 (MQ=255)
                                                                                                                                                                                                                                       aTGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGc                                          >  1:1593865/1‑184 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                    
GCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCT‑GTGATAT‑AGA‑AGACGAATGCGATTGAGATGAAAACA‑C‑CGGCG‑GTAA‑TCGCCAG‑ATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATA  >  NC_000913/954155‑954600

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