Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,361,7650GT25.4% 98.7 / 10.5 63A216S (GCG→TCG) puuDgamma‑glutamyl‑gamma‑aminobutyrate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (11/36);  new base T (15/1);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 7.95e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTC  >  NC_000913/1361519‑1362010
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
gAACAGCCTGAATTGCTGGAACCTCGGGAAGCTCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTT‑G‑CCTGATCAATCATCCTTTTTCGCt                                                                                                                                                                                                                                                    >  2:424261/1‑251 (MQ=255)
gAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCt                                                                                                                                                                                                                                                    <  1:275950/251‑1 (MQ=255)
       cTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTT‑G‑CGTCCTCAATCATTCTTTTTCGCTGGGCgg                                                                                                                                                                                                                                              >  2:245531/1‑249 (MQ=255)
                       tCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGCTA‑G‑CGGCCTCAATCATCCTTTTTCGCTGGGCGT‑ACATTGGCCCCCgggat                                                                                                                                                                                                                            >  2:134834/1‑248 (MQ=255)
                            aaGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGGCGCATGAAGTTCAGGTTGACGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGCCGTTAGCCCACGGTTGAGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAAG‑C‑TCCTCAATCCACCTTTTGCGCTGGGCGG‑ACAGTGGGACCCGGacaggga                                                                                                                                                                                                                        >  2:192916/1‑244 (MQ=255)
                             aGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGCTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTATTGAAGCTCGCTCGCCGGCTGGCTTAATTGAGGCGGTTA‑G‑GGCCATCCATCATCCTTTTGTGACGTGAGT‑TCAGTGGTACCCGGAATGGAAc                                                                                                                                                                                                                       >  1:302429/1‑250 (MQ=255)
                                tCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTT‑G‑CGTCCTCAATCATCCTTTTTCCCTGGGCGT‑ACAGTGGCACCCGGAATTGaaactg                                                                                                                                                                                                                    >  1:401856/1‑246 (MQ=255)
                                tCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCCTCAATCCTCCTTTTTCGCTGGGCGT‑ACAGTGGCCACCGGAATGGagaagt                                                                                                                                                                                                                    >  1:73950/1‑245 (MQ=255)
                                    gAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTT‑G‑CGGCATCAATCATCCTTTTGTGCTGGGGGTTAC‑GTGGGACCCGGGATGGAACAGTAGCgg                                                                                                                                                                                                               >  2:225992/1‑250 (MQ=255)
                                           cGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTAcgcg                                                                                                                                                                                                        <  2:393492/251‑1 (MQ=255)
                                               ggAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCtt                                                                                                                                                                                                     <  1:432208/250‑1 (MQ=255)
                                                     gCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATTGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTATCGAGTACGCGCTTTCGCGta                                                                                                                                                                                              <  2:303531/251‑1 (MQ=255)
                                                            gCACCGTCGCATTAAGTTCAGGTTGAAGAGGGGGGATTACTGGCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCGCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATTGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGt                                                                                                                                                                                        <  2:334774/250‑1 (MQ=255)
                                                                    gCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTT‑G‑CGGCATCCATACTCCTTTTGCGCTCGGGGG‑ACAGTGGCACCCGGGATGGGACCGGAGCGAGTACGCGCTTTCGCCGATATTTTTTGAGGGCt                                                                                                                                                                               >  2:402093/1‑251 (MQ=255)
                                                                     cATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCCTCAATTCTCCTTTTTCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACCGAAGCGGGTACGCGCTTTCGCGTATATTGTTCGAGGGCtt                                                                                                                                                                              >  1:136659/1‑251 (MQ=255)
                                                                          aGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑GGGCATCAATCACCCTTTTTCGCTGGGCGT‑ACAGTGGGACCCCGAATGGGAAAGGATCGGGGACGCGCTTTCGAGTATATTGTTTGAGGGGTTTATc                                                                                                                                                                          >  1:201677/1‑250 (MQ=255)
                                                                           gTTCAGGTTGAAGAGGGGTGAGTACTGTCTGCTTTGTTCCCTGAATGTAGCAACTTTTGGGTAAACTCTCTACAGGGACAAGGGGCGAAGGTCGTTAGCCTGCGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGTGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGGGGCTCCCGGAATGGAACAGTATCGAGTACGCGCTTTCGCGGATATTGTTCGAGGGCTTTATctc                                                                                                                                                                        <  2:281523/251‑3 (MQ=255)
                                                                                ggTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCTCCGCtt                                                                                                                                                                   <  1:52480/251‑1 (MQ=255)
                                                                                ggTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCtt                                                                                                                                                                   <  2:177468/251‑1 (MQ=255)
                                                                                 gTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATTAATCATCCTTTTGCGCTGGGGGT‑AACGTGGCCCCCCGAATGGAACAGTAGCGAGGACGCGCTTTCGCGTGTATTGTTCGAGGGGTTTTTTACCGCttt                                                                                                                                                                  >  2:404410/1‑250 (MQ=255)
                                                                                           ggggATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCat                                                                                                                                                        <  1:87303/251‑1 (MQ=255)
                                                                                                       cTGCGTTGTTAACTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGCCAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTTCGCTGGGCGG‑ACAGTGGCACCCGGAATGGAAACGTAGCGAGGAAGCGCTTTCGCGTATATTGTTCGAGGGATTTATCACCGCTTGTccgaaccctctagctgaaaa                                                                                                                                             >  1:412038/1‑231 (MQ=255)
                                                                                                              gttACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGAc                                                                                                                                     >  1:89604/1‑251 (MQ=255)
                                                                                                                tACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTTCGCTGGGCGG‑ACAGCGGGACCCGGAATGGGACAGTAGCGAGTACGCGCTTTCGCGGATATTTTTTGCGGGCTTTATTACCCCTTGTCAGCACCCTATCGCTGAAAAACAGCGActc                                                                                                                                   >  1:413248/1‑251 (MQ=255)
                                                                                                                     gAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑A‑CGTCATCAATCCTCCTTTTGGGCTGGGGGG‑ACAAAGGGCCCCGGGATGGAACACGAACGCGTACGGGCTTTCGAGGATATTGTTCGCGGGCTTTATCACCGCTTATCCGCACCATATCGCTGgaaaaaaacggctcttacc                                                                                                                              >  1:219081/1‑237 (MQ=255)
                                                                                                                      aaTGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCAGCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCa                                                                                                                             <  1:173569/251‑1 (MQ=255)
                                                                                                                           aGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACa                                                                                                                        <  1:254104/251‑1 (MQ=255)
                                                                                                                           aGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGGCATCAATCATCCTTTTGCGCTGGGGGT‑ACAGGGGCACCCGGAATGGAACAGTAGCGGGTACGCGCCTTCGCGGATATTGTTCGAGGGCTTTATTACCGCTTGTTAGCCCCATATCGCTGAAAAACAACGCCTCTGACCACTAcc                                                                                                                        >  1:389346/1‑250 (MQ=255)
                                                                                                                               aCTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGTAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGttt                                                                                                                    <  2:307100/251‑1 (MQ=255)
                                                                                                                                     gggTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGa                                                                                                               <  2:10151/250‑1 (MQ=255)
                                                                                                                                       gTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCTCCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAAt                                                                                                            <  1:465859/251‑1 (MQ=255)
                                                                                                                                         aaaCTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATg                                                                                                           <  1:136394/250‑1 (MQ=255)
                                                                                                                                         aaaCTCGCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGc                                                                                                          <  2:389346/251‑1 (MQ=255)
                                                                                                                                                ctACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGTTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTGCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAata                                                                                                     <  1:237902/249‑1 (MQ=255)
                                                                                                                                                ctACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTG‑G‑CGTCCTCAATCCACCTTTTTCGCTGGGCGT‑ACCGTTGCCCCCGGAATAGGACCGTAACACGTCCACGCTTTGGCCGCTATTGTTCGCGGGCTTTATCACCGCTTGTCAGCCCCATATCGCTGCAAAACAGCGGCCCTGACCCCTACCCTTTAAGGAAATTCACATATg                                                                                                   >  1:273602/1‑251 (MQ=255)
                                                                                                                                                ctACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCTCCATATCGCTTAAAAACAGCGACTCGGACCACTACAGTTTAAGGAAATGCAAat                                                                                                      <  2:319531/248‑1 (MQ=255)
                                                                                                                                                  aCATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATg                                                                                                   <  2:161335/249‑1 (MQ=255)
                                                                                                                                                       gACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGtgatg                                                                                            <  2:201677/251‑1 (MQ=255)
                                                                                                                                                        aCAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑GGGCATCAAACATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCCGGATGGAACAGTAGCGAGGACGCGCTTTTGGGGTTATTGTTCGAGGGCGTTATCCCCCCTTTTCAGCCCCCTATTCCTGAAAAACAAAGGCCCTTACCACTACAGTTTAAAGAAATGCAAATATGAggggtgg                                                                                           >  2:272078/1‑245 (MQ=255)
                                                                                                                                                          aaGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTTCGCTGGGCGG‑ACAGTGGCCCCCGGAAAGGAACAgg                                                                                                                                                                                                                    >  1:262577/1‑127 (MQ=255)
                                                                                                                                                          aaGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTTCGCTGGGCGG‑ACAGTGGCACCCGGAATGGAACAGGAGCGAGTACGCGCTTTCGCGGATATTGTTCGAGGGCTTTATTACCGCTTTTCAGCACCCTTTCGCTGGAAAACAACGGCTCTGACCCCTACAGTTTAAAGAAATGCAAATATGAGtggatagg                                                                                         >  1:5015/1‑245 (MQ=255)
                                                                                                                                                          aaGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGt                                                                                                                                                                                                                    <  2:262577/128‑1 (MQ=255)
                                                                                                                                                           aGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTTCGCTGGGGGG‑ACAGGGGGCCCCGGAATGGAAAAGTAGCGGGTACGCGCTTTCTCGGATATTGGTTGAGGGCTTTATCACCGCTTGGCAGCACCCAATTGCTGGAAAACAAGGGCTTTTACCACTACAGGTTAAGGAAATTCAAATAAGAGTGATGAggg                                                                                        >  1:206995/1‑251 (MQ=255)
                                                                                                                                                             gggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGAc                                                                                      <  1:326500/251‑1 (MQ=255)
                                                                                                                                                              ggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACt                                                                                     <  1:369201/251‑1 (MQ=255)
                                                                                                                                                              ggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATAGCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACt                                                                                     <  2:177463/251‑1 (MQ=255)
                                                                                                                                                                 gAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACGCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGgc                                                                                  <  2:73950/251‑1 (MQ=255)
                                                                                                                                                                       cGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑A‑CGTCATCAATCATCCTTTTTCGCTGGGCGG‑ACAGTGGCCCCCGGAATGGAACAGTAGCGAGGACGCGCTTTCGCGGATATTTTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAACGACTCTGGCCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGcccgg                                                                            >  1:374860/1‑248 (MQ=255)
                                                                                                                                                                          tAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCAGCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGaa                                                                          <  2:219081/250‑1 (MQ=255)
                                                                                                                                                                           aGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGaaaa                                                                        <  2:89604/251‑1 (MQ=255)
                                                                                                                                                                                            tGAAGGTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTT‑GAGGTAAT‑ACTCATCCTTTTTAGCTGGGGGG‑AACGTGGGACCAGGGATGGAACAAGATCGAGTACGCGCTTTTGCGTATATTTTTCGGGGGCTTTCTCACCCATTTGCAGCACAATATCGCTGGAAAACAGAGCCTCTGACCAACACCGGTTAAGGAAATGTAAAAATAGGTGATGAGGGACAGGCGGCATGAGAAAGCTTGTCggtataacg                                                       >  2:489165/1‑244 (MQ=255)
                                                                                                                                                                                             gAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCg                                                       <  2:136659/250‑1 (MQ=255)
                                                                                                                                                                                                       cgcCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAgg                                            <  2:401856/251‑1 (MQ=255)
                                                                                                                                                                                                            gATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCtt                                       <  1:70788/251‑1 (MQ=255)
                                                                                                                                                                                                             aTGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCttt                                      <  1:272078/251‑1 (MQ=255)
                                                                                                                                                                                                                   tAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTAATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACaa                                 <  1:77443/250‑1 (MQ=255)
                                                                                                                                                                                                                     gTTGAGGCGGTTA‑G‑CGTCATCAATCAGCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTGTAAGGAAATGCAAATATGAGTGAGGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAAcgt                              <  2:169668/251‑1 (MQ=255)
                                                                                                                                                                                                                        gAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAAcgtcgt                           <  2:417473/251‑1 (MQ=255)
                                                                                                                                                                                                                                 a‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGa                   >  2:185679/1‑250 (MQ=255)
                                                                                                                                                                                                                                           caatcaTCCTTTTTCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCgg          <  2:413248/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 tcCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGCATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACAGCTTGTCAGCACCATATCGCTGACAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGAc    >  1:434347/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  ccTTGTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATTAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACt   <  2:246807/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     tttGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTc  <  1:402093/249‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
GAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTA‑G‑CGTCATCAATCATCCTTTTGCGCTGGGCGT‑ACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTC  >  NC_000913/1361519‑1362010

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