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.

GAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTC  >  NC_000913/1361519‑1362010
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
gAACAGCCTGAATTGCTGGAACCTCGGGAAGCTCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTTGCCTGATCAATCATCCTTTTTCGCt                                                                                                                                                                                                                                                    >  2:424261/1‑251 (MQ=255)
gAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCt                                                                                                                                                                                                                                                    <  1:275950/251‑1 (MQ=255)
       cTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTTGCGTCCTCAATCATTCTTTTTCGCTGGGCgg                                                                                                                                                                                                                                              >  2:245531/1‑249 (MQ=255)
                       tCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGCTAGCGGCCTCAATCATCCTTTTTCGCTGGGCGTACATTGGCCCCCgggat                                                                                                                                                                                                                            >  2:134834/1‑248 (MQ=255)
                            aaGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGGCGCATGAAGTTCAGGTTGACGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGCCGTTAGCCCACGGTTGAGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAAGCTCCTCAATCCACCTTTTGCGCTGGGCGGACAGTGGGACCCGGacaggga                                                                                                                                                                                                                        >  2:192916/1‑244 (MQ=255)
                             aGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGCTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTATTGAAGCTCGCTCGCCGGCTGGCTTAATTGAGGCGGTTAGGGCCATCCATCATCCTTTTGTGACGTGAGTTCAGTGGTACCCGGAATGGAAc                                                                                                                                                                                                                       >  1:302429/1‑250 (MQ=255)
                                tCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTTGCGTCCTCAATCATCCTTTTTCCCTGGGCGTACAGTGGCACCCGGAATTGaaactg                                                                                                                                                                                                                    >  1:401856/1‑246 (MQ=255)
                                tCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCCTCAATCCTCCTTTTTCGCTGGGCGTACAGTGGCCACCGGAATGGagaagt                                                                                                                                                                                                                    >  1:73950/1‑245 (MQ=255)
                                    gAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTTGCGGCATCAATCATCCTTTTGTGCTGGGGGTTACGTGGGACCCGGGATGGAACAGTAGCgg                                                                                                                                                                                                               >  2:225992/1‑250 (MQ=255)
                                           cGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTAcgcg                                                                                                                                                                                                        <  2:393492/251‑1 (MQ=255)
                                               ggAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCtt                                                                                                                                                                                                     <  1:432208/250‑1 (MQ=255)
                                                     gCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATTGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTATCGAGTACGCGCTTTCGCGta                                                                                                                                                                                              <  2:303531/251‑1 (MQ=255)
                                                            gCACCGTCGCATTAAGTTCAGGTTGAAGAGGGGGGATTACTGGCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCGCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATTGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGt                                                                                                                                                                                        <  2:334774/250‑1 (MQ=255)
                                                                    gCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTTGCGGCATCCATACTCCTTTTGCGCTCGGGGGACAGTGGCACCCGGGATGGGACCGGAGCGAGTACGCGCTTTCGCCGATATTTTTTGAGGGCt                                                                                                                                                                               >  2:402093/1‑251 (MQ=255)
                                                                     cATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCCTCAATTCTCCTTTTTCGCTGGGCGTACAGTGGCACCCGGAATGGAACCGAAGCGGGTACGCGCTTTCGCGTATATTGTTCGAGGGCtt                                                                                                                                                                              >  1:136659/1‑251 (MQ=255)
                                                                          aGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGGGGCATCAATCACCCTTTTTCGCTGGGCGTACAGTGGGACCCCGAATGGGAAAGGATCGGGGACGCGCTTTCGAGTATATTGTTTGAGGGGTTTATc                                                                                                                                                                          >  1:201677/1‑250 (MQ=255)
                                                                           gTTCAGGTTGAAGAGGGGTGAGTACTGTCTGCTTTGTTCCCTGAATGTAGCAACTTTTGGGTAAACTCTCTACAGGGACAAGGGGCGAAGGTCGTTAGCCTGCGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGTGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGGGGCTCCCGGAATGGAACAGTATCGAGTACGCGCTTTCGCGGATATTGTTCGAGGGCTTTATctc                                                                                                                                                                        <  2:281523/251‑3 (MQ=255)
                                                                                ggTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCTCCGCtt                                                                                                                                                                   <  1:52480/251‑1 (MQ=255)
                                                                                ggTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCtt                                                                                                                                                                   <  2:177468/251‑1 (MQ=255)
                                                                                 gTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATTAATCATCCTTTTGCGCTGGGGGTAACGTGGCCCCCCGAATGGAACAGTAGCGAGGACGCGCTTTCGCGTGTATTGTTCGAGGGGTTTTTTACCGCttt                                                                                                                                                                  >  2:404410/1‑250 (MQ=255)
                                                                                           ggggATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCat                                                                                                                                                        <  1:87303/251‑1 (MQ=255)
                                                                                                       cTGCGTTGTTAACTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGCCAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTTCGCTGGGCGGACAGTGGCACCCGGAATGGAAACGTAGCGAGGAAGCGCTTTCGCGTATATTGTTCGAGGGATTTATCACCGCTTGTccgaaccctctagctgaaaa                                                                                                                                             >  1:412038/1‑231 (MQ=255)
                                                                                                              gttACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGAc                                                                                                                                     >  1:89604/1‑251 (MQ=255)
                                                                                                                tACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTTCGCTGGGCGGACAGCGGGACCCGGAATGGGACAGTAGCGAGTACGCGCTTTCGCGGATATTTTTTGCGGGCTTTATTACCCCTTGTCAGCACCCTATCGCTGAAAAACAGCGActc                                                                                                                                   >  1:413248/1‑251 (MQ=255)
                                                                                                                     gAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAACGTCATCAATCCTCCTTTTGGGCTGGGGGGACAAAGGGCCCCGGGATGGAACACGAACGCGTACGGGCTTTCGAGGATATTGTTCGCGGGCTTTATCACCGCTTATCCGCACCATATCGCTGgaaaaaaacggctcttacc                                                                                                                              >  1:219081/1‑237 (MQ=255)
                                                                                                                      aaTGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCAGCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCa                                                                                                                             <  1:173569/251‑1 (MQ=255)
                                                                                                                           aGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACa                                                                                                                        <  1:254104/251‑1 (MQ=255)
                                                                                                                           aGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGGCATCAATCATCCTTTTGCGCTGGGGGTACAGGGGCACCCGGAATGGAACAGTAGCGGGTACGCGCCTTCGCGGATATTGTTCGAGGGCTTTATTACCGCTTGTTAGCCCCATATCGCTGAAAAACAACGCCTCTGACCACTAcc                                                                                                                        >  1:389346/1‑250 (MQ=255)
                                                                                                                               aCTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGTAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGttt                                                                                                                    <  2:307100/251‑1 (MQ=255)
                                                                                                                                     gggTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGa                                                                                                               <  2:10151/250‑1 (MQ=255)
                                                                                                                                       gTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCTCCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAAt                                                                                                            <  1:465859/251‑1 (MQ=255)
                                                                                                                                         aaaCTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATg                                                                                                           <  1:136394/250‑1 (MQ=255)
                                                                                                                                         aaaCTCGCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGc                                                                                                          <  2:389346/251‑1 (MQ=255)
                                                                                                                                                ctACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGTTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTGCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAata                                                                                                     <  1:237902/249‑1 (MQ=255)
                                                                                                                                                ctACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTGGCGTCCTCAATCCACCTTTTTCGCTGGGCGTACCGTTGCCCCCGGAATAGGACCGTAACACGTCCACGCTTTGGCCGCTATTGTTCGCGGGCTTTATCACCGCTTGTCAGCCCCATATCGCTGCAAAACAGCGGCCCTGACCCCTACCCTTTAAGGAAATTCACATATg                                                                                                   >  1:273602/1‑251 (MQ=255)
                                                                                                                                                ctACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCTCCATATCGCTTAAAAACAGCGACTCGGACCACTACAGTTTAAGGAAATGCAAat                                                                                                      <  2:319531/248‑1 (MQ=255)
                                                                                                                                                  aCATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATg                                                                                                   <  2:161335/249‑1 (MQ=255)
                                                                                                                                                       gACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGtgatg                                                                                            <  2:201677/251‑1 (MQ=255)
                                                                                                                                                        aCAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGGGGCATCAAACATCCTTTTGCGCTGGGCGTACAGTGGCACCCCGGATGGAACAGTAGCGAGGACGCGCTTTTGGGGTTATTGTTCGAGGGCGTTATCCCCCCTTTTCAGCCCCCTATTCCTGAAAAACAAAGGCCCTTACCACTACAGTTTAAAGAAATGCAAATATGAggggtgg                                                                                           >  2:272078/1‑245 (MQ=255)
                                                                                                                                                          aaGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTTCGCTGGGCGGACAGTGGCCCCCGGAAAGGAACAgg                                                                                                                                                                                                                    >  1:262577/1‑127 (MQ=255)
                                                                                                                                                          aaGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTTCGCTGGGCGGACAGTGGCACCCGGAATGGAACAGGAGCGAGTACGCGCTTTCGCGGATATTGTTCGAGGGCTTTATTACCGCTTTTCAGCACCCTTTCGCTGGAAAACAACGGCTCTGACCCCTACAGTTTAAAGAAATGCAAATATGAGtggatagg                                                                                         >  1:5015/1‑245 (MQ=255)
                                                                                                                                                          aaGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGt                                                                                                                                                                                                                    <  2:262577/128‑1 (MQ=255)
                                                                                                                                                           aGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTTCGCTGGGGGGACAGGGGGCCCCGGAATGGAAAAGTAGCGGGTACGCGCTTTCTCGGATATTGGTTGAGGGCTTTATCACCGCTTGGCAGCACCCAATTGCTGGAAAACAAGGGCTTTTACCACTACAGGTTAAGGAAATTCAAATAAGAGTGATGAggg                                                                                        >  1:206995/1‑251 (MQ=255)
                                                                                                                                                             gggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGAc                                                                                      <  1:326500/251‑1 (MQ=255)
                                                                                                                                                              ggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACt                                                                                     <  1:369201/251‑1 (MQ=255)
                                                                                                                                                              ggCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATAGCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACt                                                                                     <  2:177463/251‑1 (MQ=255)
                                                                                                                                                                 gAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACGCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGgc                                                                                  <  2:73950/251‑1 (MQ=255)
                                                                                                                                                                       cGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAACGTCATCAATCATCCTTTTTCGCTGGGCGGACAGTGGCCCCCGGAATGGAACAGTAGCGAGGACGCGCTTTCGCGGATATTTTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAACGACTCTGGCCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGcccgg                                                                            >  1:374860/1‑248 (MQ=255)
                                                                                                                                                                          tAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCAGCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGaa                                                                          <  2:219081/250‑1 (MQ=255)
                                                                                                                                                                           aGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGaaaa                                                                        <  2:89604/251‑1 (MQ=255)
                                                                                                                                                                                            tGAAGGTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTTGAGGTAATACTCATCCTTTTTAGCTGGGGGGAACGTGGGACCAGGGATGGAACAAGATCGAGTACGCGCTTTTGCGTATATTTTTCGGGGGCTTTCTCACCCATTTGCAGCACAATATCGCTGGAAAACAGAGCCTCTGACCAACACCGGTTAAGGAAATGTAAAAATAGGTGATGAGGGACAGGCGGCATGAGAAAGCTTGTCggtataacg                                                       >  2:489165/1‑244 (MQ=255)
                                                                                                                                                                                             gAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCg                                                       <  2:136659/250‑1 (MQ=255)
                                                                                                                                                                                                       cgcCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAgg                                            <  2:401856/251‑1 (MQ=255)
                                                                                                                                                                                                            gATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCtt                                       <  1:70788/251‑1 (MQ=255)
                                                                                                                                                                                                             aTGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCttt                                      <  1:272078/251‑1 (MQ=255)
                                                                                                                                                                                                                   tAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTAATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACaa                                 <  1:77443/250‑1 (MQ=255)
                                                                                                                                                                                                                     gTTGAGGCGGTTAGCGTCATCAATCAGCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTGTAAGGAAATGCAAATATGAGTGAGGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAAcgt                              <  2:169668/251‑1 (MQ=255)
                                                                                                                                                                                                                        gAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAAcgtcgt                           <  2:417473/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGa                   >  2:185679/1‑250 (MQ=255)
                                                                                                                                                                                                                                           caatcaTCCTTTTTCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCgg          <  2:413248/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 tcCTTTTGCGCTGGGCGTACAGTGGCACCCGGCATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACAGCTTGTCAGCACCATATCGCTGACAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGAc    >  1:434347/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  ccTTGTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATTAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACt   <  2:246807/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     tttGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTc  <  1:402093/249‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
GAACAGCCTGAATTGCTGGAACATCGGGAAGATCCCGAACTGCCGGTGGAACAGCAATATGCACCGTCGCATGAAGTTCAGGTTGAAGAGGGGGGATTACTGTCTGCGTTGTTACCTGAATGTAGCAACTTTTGGGTAAACTCTCTACATGGACAAGGGGCGAAGGTCGTTAGCCCACGGTTGCGTGTTGAAGCTCGCTCGCCGGATGGCTTAGTTGAGGCGGTTAGCGTCATCAATCATCCTTTTGCGCTGGGCGTACAGTGGCACCCGGAATGGAACAGTAGCGAGTACGCGCTTTCGCGTATATTGTTCGAGGGCTTTATCACCGCTTGTCAGCACCATATCGCTGAAAAACAGCGACTCTGACCACTACAGTTTAAGGAAATGCAAATATGAGTGATGAGGGACTGGCGCCAGGAAAACGCTTGTCGGAAATCCGCCAGCAGCAGGGGCTTTCACAACGTCGTGCCGCCGAACTCTCCGGGCTGACTC  >  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: