Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913954,3880CT31.5% 60.2 / 17.6 56A27T (GCA→ACA) focAformate channel
Reads supporting (aligned to +/- strand):  ref base C (29/8);  new base T (0/17);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 2.29e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.52e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAAT  >  NC_000913/954152‑954615
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                    
gATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCt                                                                                                                                                                                                                                     <  2:687401/238‑1 (MQ=255)
gATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCt                                                                                                                                                                                                                                     <  2:886346/238‑1 (MQ=255)
   gcgcCCACTCGTCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCCGATAGAAAGGCTTAAGCGG‑ATGTTTCGTTGCTTtat                                                                                                                                                                                                                                 >  1:979170/1‑239 (MQ=255)
      cccACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTTTAg                                                                                                                                                                                                                               >  1:818836/1‑238 (MQ=255)
       ccACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCt                                                                                                                                                                                                                                     <  1:841478/231‑1 (MQ=255)
                 tAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCCGAAAGAAAAGCTTAAACGG‑ATGGTTCGTTGCTTTATAGa‑aac                                                                                                                                                                                                                          >  2:841478/1‑231 (MQ=255)
                                 cacGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGc                                                                                                                                                                                                    >  1:241721/1‑237 (MQ=255)
                                    ggTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCTCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGCGATGAAA‑ACACCGGCGGTAAACGCCAGATAGAAAGGCTTAAGCGG‑ATGTTTCGTTGCTTTTTAGA‑CACCCGCCCCTTCTGTCACTTTGGCCAtt                                                                                                                                                                                                >  1:697022/1‑238 (MQ=255)
                                    ggTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑AACCCGCCCTTTCGGGCACCTTGGCCAtt                                                                                                                                                                                                >  1:608561/1‑238 (MQ=255)
                                     gTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGAAAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTgc                                                                                                                                                                                              >  2:725721/1‑239 (MQ=255)
                                      tGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTTGCCCTTgca                                                                                                                                                                                             >  1:45795/1‑239 (MQ=255)
                                         aaGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGGAATCGCCAGATAGAAAGTCTTAAATGG‑ATGTTTTGTTGGTTTATAGA‑CACCCGCCTCTTCTGCCCCTTTGGGCATTgccgtc                                                                                                                                                                                          >  1:69207/1‑235 (MQ=255)
                                           taaaaGAGCGCGCCCCCGAAAACAACACAAAGAATCAGCCCCAGAGAAAAAAAAATG‑CCGCCAACCATTTTTG‑CCATGCCGGGGGACATGGT‑CCTGGGCCAGGTGGTTGCGGTAAAAAAGAAGACGAATGCAATGGGAATAAAA‑CCACCGGCGGTAATCGCCAAAAAAAAAGTCTTAGGGGA‑ATTTTTCGTTGTTTTAGAAA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAgg                                                                                                                                                                                        <  2:611732/237‑1 (MQ=255)
                                              aaGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAACGG‑GTGTTTCGTTGCTTTATAGA‑CACCCCCCCCTTCGGCCACTTTGGCCCTTGCAGCAGGAAg                                                                                                                                                                                     >  1:44759/1‑239 (MQ=255)
                                                gagaTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGGCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCCCTTCGGCCACTTTGGGCATTGCAGGAGGAAGTa                                                                                                                                                                                   >  2:528974/1‑239 (MQ=255)
                                                        ccccGCAAACAACAAAAAGAACCAGCCCAGGGAAGAAG‑AAATGCCCGCCAACCTGTTTCG‑CCTAGCCGAGGGGATTTGTCCCTGCCCCA‑GTGGTTGTTGAAATAAAAAAGACAAATGCGATTGAGAAAAAA‑CCCCCGGCGTAAACCGCCAGAAAAAAATTTTAAAGGGG‑TTTTTTCGTTGTTAAAAAAC‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaaaa                                                                                                                                                                                <  2:710434/232‑1 (MQ=255)
                                                             cAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATGGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCCGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCTGCCACTTTGGCCCTTGCAGCAGGAAGTAAAAGa                                                                                                                                                                              >  1:938153/1‑231 (MQ=255)
                                                                              cAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATTCCCAGATNGAAAGGCTTAAGCGG‑CTGGTTCGTTGCTTTAAAGA‑CACCCGcc                                                                                                                                                                                                                     >  1:477354/1‑175 (MQ=255)
                                                                                 ccccAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTTCAGCAGGAAGTAAAAGATAAAAAGGGTTGTCAGATTTCACACTaa                                                                                                                                                  >  2:635524/1‑239 (MQ=255)
                                                                                  cccAGAGAGAAGCAAATC‑CCGCAAACCAGTTTTG‑CCATGCGGAAGGCCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑CCCCCGGCGGTAATCGCCAGATAGAAAGTCTAAAGCGG‑ATGTTTCTTTGTTTTATAGC‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCacac                                                                                                                                                     <  1:840157/235‑1 (MQ=255)
                                                                                  ccaagagagAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATAAAA‑CCCCCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCttt                                                                                                                                                          <  1:821601/227‑1 (MQ=255)
                                                                                           aagcagccg‑ccgccAA‑CTGTTTTGCCCATGCAGAGGGCGATTGGGCCGGGGCCG‑GGGGTGGTTGTAAGAAAAAAGACGAAGGCGATTGAGAGAAAA‑ACCCGGGCGTTAACCGCCAAAAAAAAAGTCTAAGGCGG‑TTTTTTCTTTGTTTTAAAGA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCtt                                                                                                                                         <  2:632438/230‑1 (MQ=255)
                                                                                           aaGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTa                                                                                                                                                   >  1:381533/1‑228 (MQ=255)
                                                                                           aaGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGGAAAAGATCAAAAAGGTTGTCAGCTTTCACACTa                                                                                                                                                   >  1:580065/1‑228 (MQ=255)
                                                                                            aGCAAATG‑CCGCCACCCAGTTTGC‑CAAGGCGGAAGGGCTTGGGGCTTGTGCAG‑GGGGTGGCTGTAAAAAAAAAGACGAAGGCGATGGAGATAAAA‑CCCCCGGCGGTAATCCCCAGATAGAAATTCTTAAGGGG‑ATTTTTCGTTGTTTTATAGA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctc                                                                                                                                             <  2:608561/233‑1 (MQ=255)
                                                                                              cAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATGTAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGGAATCGCCAAGTAAAAAATCTTAAACGG‑AAGTTTCCTTGCTTTATAGA‑CACACGCCTCTTCGGCCACTTTGGTCATTGCCGCAGGAAGAAAAAGATCAAAAAGGGTGTCAGGttttcaacttaatc                                                                                                                                               >  2:926595/1‑218 (MQ=255)
                                                                                                                gtttccccAGGCCAAAGGGAATTGTCCCTGTGCCG‑GTGTTTGCTGTAATAAAAAAGCCGAATGCGTTTGAGATGAAA‑ACCCCGGCGGTAACCGCCAGATAAAAAGTCTTAAGCGG‑TTGTTTCGTTGTTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGAt                                                                                                                     <  1:560731/233‑1 (MQ=255)
                                                                                                                       cATGCCGAGGGGCATTGGGCCTGGCCAA‑GGGGTGGGTGTAATAAAGAAGAGGAAGGCGATGGAGAAAAAA‑CCCCCGGCGATAACCGCAAAGAAAAAAGTCTTAGGCGG‑TTGTTTCGTTGCTTTATAGA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                                              <  2:241721/239‑1 (MQ=255)
                                                                                                                       cATGCCGAAGGTCATTGGGCCTGTGCCA‑GTGGTTGCTGTGATAAAAAAGACGAATGCGTTTGAGATAAAA‑CCCCCGGCGTTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTTTTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                                              <  1:528974/239‑1 (MQ=255)
                                                                                                                           acGAAGGACTTTGTGCCGGTGAGA‑GGTGTTTGTATGATATAAAAAAAGAATGCTATTGAGATAAAA‑CCCCCGGCTGTAACCCCCAGATAGAAAGTCTTAAGGGG‑ATTTTGCTTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                                              <  2:39132/234‑1 (MQ=255)
                                                                                                                                     ttGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGACAGAAAGTCTTAAGCGG‑ATTTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGGCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATg                                                                                                >  1:707012/1‑239 (MQ=255)
                                                                                                                                          gcttgtgcagttGGTTGTGGTAATATAAAAGAGGAAGGGGATGGAGATAAAA‑CCACCGGCGGTAACCGCAAGAAAAAAAGTCTAAGGCGG‑ATTTTTCGTTGTTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATgaga                                                                                             <  2:979170/227‑1 (MQ=255)
                                                                                                                                           cggtgccA‑GTGGTTGCTGTGATATAGAAGACGAAGGCGATTGAGATGAAA‑CCCCCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑TTTTTTCTTTGTTTTATAGA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGaaa                                                                                           <  1:725721/236‑1 (MQ=255)
                                                                                                                                                      gtttgttttaatatAGATGACAAATGCGATTGAG‑TGAAAGTCCCCGGCGGTATTCCCCAGATGAAAAGTCTTAGGCGG‑ATTTTTGTTTGTTTAATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGAt                                                                                <  1:898266/229‑1 (MQ=255)
                                                                                                                                                       tttgcgtggATATAAAAGAGGAGTGAGATTGAGATAAAA‑ACACCGGCGGTAATCGCAAGATAAAAATTTTAAAGGGG‑ATGTTTCGTTGTTTAATAGA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGAt                                                                                <  2:485530/230‑1 (MQ=255)
                                                                                                                                                          gCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGGGG‑GTGTTTTCTTGCTTTTTAGA‑CACCCGCCCCTTCCGCCCATTTGGCCATTGCAGCCGGAAGTAAAAGATCAAAAGGGGTTTCAGCTTTTCCACTTACTCTCTCTTTATTAATTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATCAATCata                                                                            >  2:973476/1‑239 (MQ=255)
                                                                                                                                                           ttgagATATAGAAAACAAATGCGATTGAAATGAAA‑ACCCCGGGGGAAATCGCCAGAAAAAAAGTCTAAAGGGT‑TTTTTTGGTTTCTTTAAAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAAct                                                                                                                                                <  2:119996/166‑1 (MQ=255)
                                                                                                                                                            tgtgATAAAGAAGACGAATGCGATTGAAATAAAA‑CCCCCGGCGGTAATCGCAAGATAGAAAGTTTTAAGGGT‑ATGTTTCTTTGTTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCatat                                                                           <  1:802953/238‑1 (MQ=255)
                                                                                                                                                            ggagATATAGAAGTAAAGTGCGATGGAGAAAAAA‑ACCCCGGGGGAAATCGCAAGAAAGAAACTGTAAAGGGG‑TTGTTGTGTTGTTGAAAAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCatat                                                                           <  2:32595/235‑1 (MQ=255)
                                                                                                                                                              taatatAGAAGACGAATGCGATTGAGATGAAA‑CCCCGGGCGGTAATCGCAAGATAGAAAGTCTAAGGCGG‑ATGTTTGGTTGTTTTATAGA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                                              <  1:618540/199‑1 (MQ=255)
                                                                                                                                                                   aaaaagaGGATTGCGTTTGGAATAAAA‑ACCCCGCCGGTAAGCGCCAGATAAAAAGTTTTAAGGGG‑TTGTTTCTTTGTTTTAAAGA‑ACCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATga                                                                                               <  2:812286/208‑1 (MQ=255)
                                                                                                                                                                     gaagaCGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTttat                                                                                                                                      >  1:119996/1‑170 (MQ=255)
                                                                                                                                                                        gaCGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCCGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCg                                                                       >  1:32595/1‑230 (MQ=255)
                                                                                                                                                                        gaCGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTta                                                              >  1:401390/1‑239 (MQ=255)
                                                                                                                                                                        gaCGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGGCTTAAGCGG‑AAGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTCACACTTACTCTCTCTTTATTAAATCGGCGACGAGATTCTAACAAAACATtata                                                                                                    >  2:618540/1‑201 (MQ=255)
                                                                                                                                                                             aTGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                                     >  1:812286/1‑211 (MQ=255)
                                                                                                                                                                                  aTTGAGATGAAT‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTTGCCTGGCTTATAGGCCCGTa                                                    >  1:997440/1‑239 (MQ=255)
                                                                                                                                                                                                       tgTAATCGCAAGATAG‑AAGTGATAAGGGGAAGGTTTCG‑TGTTTTTAAGACCACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCCCACTAACTCTCTCTTTATTAAGGCGGCGACGAGATACTAACAAAGCAGTATAGATGAGAAATTGATATAGAGCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCAAGGTTTTTATGc                                 <  1:973476/237‑1 (MQ=255)
                                                                                                                                                                                                           aTCCCCAGTTAGAAAGTCTTAAGCGG‑ATTTTTGTTTCTTTTAAAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa                            <  1:926595/239‑1 (MQ=255)
                                                                                                                                                                                                                         aGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATtata                                                                                                    >  1:669990/1‑153 (MQ=255)
                                                                                                                                                                                                                            ttttaagggtttGTTTTTTTGATTTAAACA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGa            <  2:580065/228‑1 (MQ=255)
                                                                                                                                                                                                                             ttAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAatat                >  1:175818/1‑233 (MQ=255)
                                                                                                                                                                                                                             ttAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGAtt          >  1:923712/1‑239 (MQ=255)
                                                                                                                                                                                                                                      aTGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCAt                                            >  2:160704/1‑197 (MQ=255)
                                                                                                                                                                                                                                      aTGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCAt                                            <  1:160704/197‑1 (MQ=255)
                                                                                                                                                                                                                                      aTGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCaaat  >  2:529454/1‑239 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                    
GATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAA‑ACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGG‑ATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAAT  >  NC_000913/954152‑954615

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