Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,828,5060TG30.8% 84.5 / 10.9 69V257G (GTG→GGG) srlDsorbitol‑6‑phosphate 2‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (11/34);  new base G (18/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.05e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACGGCATTACGGTGCATTCACTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGG‑CAAATTTCTCGTTTTAACCGTGCCTTCGACACACTATGCCAGCAAGGGCGGGTTGGCGTGGGCCGTTC  >  NC_000913/2828260‑2828754
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                
aCGGCATTACGGTGCATTCACTGATGCTCGGTAACCTTCTGAAATCGCCGATGTTCCAGTCACTGGTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATg                                                                                                                                                                                                                                                              <  2:200272/251‑1 (MQ=255)
aCGGCATTACGGTGCATTCACTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATg                                                                                                                                                                                                                                                              <  1:93086/251‑1 (MQ=255)
aCGGCATTACGGTGCATTCACTGAGGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGAt                                                                                                                                                                                                                                                               <  2:862508/250‑1 (MQ=255)
  ggCATTACGGTGCATTCACTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTTCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGtt                                                                                                                                                                                                                                                            >  1:542445/1‑251 (MQ=255)
     aTTACGGTGCATTCACGGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑g                                                                                                                                                                                                                                                        <  1:700417/251‑1 (MQ=255)
        aCGGTGCATTCACTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GAt                                                                                                                                                                                                                                                      <  1:791663/250‑1 (MQ=255)
        aCGGTGCATTCACTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATc                                                                                                                                                                                                                                                     <  1:802087/251‑1 (MQ=255)
         cGGTGCATTCACTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCT‑TATc                                                                                                                                                                                                                                                     >  1:717421/1‑250 (MQ=255)
             gCATTCACTGAGGCTTGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGTGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACa                                                                                                                                                                                                                                                 <  2:558806/250‑1 (MQ=255)
                  cACTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCT‑GGTCAACAGCg‑ggg                                                                                                                                                                                                                                          >  1:628957/1‑249 (MQ=255)
                  cACTGATGCCTGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑ga                                                                                                                                                                                                                                           <  1:810057/250‑1 (MQ=255)
                   aCTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGCTCAATGTAACCGGCGGGCAGGGGATGTTCT‑GATCaaacgcgg                                                                                                                                                                                                                                             >  2:493026/1‑242 (MQ=255)
                   aCTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGGTGTTCT‑GATCaaacaccg                                                                                                                                                                                                                                             >  1:322310/1‑242 (MQ=255)
                         ggccgggagccTGCTAAAATCGCCGATGTTCCAGTCACGGGTGCCACAATCCGCGACCAAGCTGGGTATCAAACCGGCTCAAGTCGAGCAGTATTACATCGACAAAGTTCCGCTCAAACGCGGCCGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCAt                                                                                                                                                                                                                                   <  2:720653/241‑1 (MQ=255)
                         gCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGAACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGTCGGTCGGGTGATGTTCC‑TATCAACAGCG‑GGGA‑Tcccat                                                                                                                                                                                                                                  >  1:170416/1‑248 (MQ=255)
                            cGGTATCCTGCTGAAATCGACGATGTTCCAGTCACTGTTGCCACAAGACGCGACCAAGCTGGGGATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCTATTATCAAGATGTGCTGAATATGCTGCTGTTCTACTCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑GCCATTatg                                                                                                                                                                                                                               <  2:610348/251‑3 (MQ=255)
                                aaCCTGCTTAAATCGCCGATGTTCCAGTCCCTGTGTCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGAt                                                                                                                                                                                                                            <  2:721546/250‑1 (MQ=255)
                                     gctgAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCg                                                                                                                                                                                                                       <  1:316514/250‑1 (MQ=255)
                                            tCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAAAAGCG‑GAGC‑TCCATTAAGGGTCTCCGTGAGA‑Ct                                                                                                                                                                                                               >  1:590793/1‑250 (MQ=255)
                                                cGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGCCAGGAGAGGTTCT‑GATCAAAAGCG‑GAGA‑TCCATTAAgggtgtccgcgggaatat                                                                                                                                                                                                              >  1:46000/1‑232 (MQ=255)
                                                cGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTAACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGa                                                                                                                                                                                                          <  1:493026/251‑1 (MQ=255)
                                                  atgtgcccggaaCTGTTGCCACAATACGCGATCAAGCTGGGTATCAAACCGGATCAAGGCGAGCAGTAGTACAGCGACAAAGTACCGCTCAAACGCGGCTGCTATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGGCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAAt                                                                                                                                                                                                        <  2:459427/240‑1 (MQ=255)
                                                  aTGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGCTTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGAAGTTCT‑GATCAACAACG‑GAGA‑ACCATTAAGGAGCTCCGTGGGA‑CCATAAAAt                                                                                                                                                                                                        >  1:557782/1‑251 (MQ=255)
                                                  aTGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATTCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGGTGGTCT‑TATCAACCACG‑GGGA‑TCCATTTAGGGTCTCCGTGGGA‑CTATAGaa                                                                                                                                                                                                         >  1:573963/1‑250 (MQ=255)
                                                               ggTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCTCCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTa                                                                                                                                                                                           <  2:397388/250‑1 (MQ=255)
                                                                gtttccACAATACGCGACCATGCAGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGGTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTTAGA‑CTATAGAATGCCTGATGCGCt                                                                                                                                                                                            <  2:451883/245‑1 (MQ=255)
                                                                gTTGCCACAATACGCGAGCAAGCTGGGGATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTTT‑GGTCAACAACG‑GGGA‑TCCATTAAGGAACTCCGTGAGC‑CTATAGAATTCCTGCTTCGCTAc                                                                                                                                                                                          >  1:702637/1‑251 (MQ=255)
                                                                   gCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCt                                                                                                                                                                                       <  1:354518/251‑1 (MQ=255)
                                                                      cacaatcggcgacAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCtcat                                                                                                                                                                                    <  2:307316/239‑1 (MQ=255)
                                                                                aCCAAGCTGGTTATCAAACCGGATCAAGTCGAGCCGTATTACATCGACAAAGCTCCGCTCAAACGCGGCTTCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTATGGCGTCGTACTTCTCCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACa                                                                                                                                                                          <  2:7048/251‑1 (MQ=255)
                                                                                 ccAAGCTGGGTATCAAACCGGAGCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCGCATCAGGCATACAg                                                                                                                                                                         <  1:135407/251‑1 (MQ=255)
                                                                                        ggttaacAAACCGGAGCATGTCGAGCAGTATTACATCGACAAAGTACCACTAAAAAGCTGCTGCGATGTGCCAGATGTGATGAATACGGTGATGTTCTAGTCCAGACCTAAGGCGTCGAACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGTTGATGCTCT‑GATCAACAGCG‑TAGA‑TCCATTAAGGATCTTCGTGAGA‑CTATAGAATGCCTGATGCGCAACGTGCATCAGTCATACAGG‑ACTTct                                                                                                                                                                 <  2:171186/245‑2 (MQ=255)
                                                                                        gggTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGG‑ACTTcc                                                                                                                                                                 <  1:637454/251‑1 (MQ=255)
                                                                                            ataaaaCCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCGGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑GCCATTAAGGATCGCCGTGAGA‑CTATAGAATGCCGGATGCGCTACGCTCATCAGGCATACAGG‑ACTTCCGCCa                                                                                                                                                             <  2:532199/248‑1 (MQ=255)
                                                                                                        cAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCT‑GATCAACAGCG‑GAGA‑TCCCTTAAGGATCTCCGTGAGA‑CTATAGAATGCCTTATGCGCTACGCTCCTCAGGCATACAGG‑ACTTCCGCCCCTACAT‑TAAgg                                                                                                                                                 >  1:31780/1‑250 (MQ=255)
                                                                                                             cGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGGCCGGTGGTGTTCT‑TGTCAACAGCG‑GAGG‑TCCATTAAAGATCTCCGTGAGA‑CTATAGGATGCCTGATGCGCTACCCTCCTCAGGCATACAGGCACT‑CCGCCCCCACAT‑TAAGGaaaaat                                                                                                                                           >  2:654355/1‑249 (MQ=255)
                                                                                                                  aGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGGGA‑CTATAGAATGCCTGATGCGCCACGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGGTATgg                                                                                                                                      >  1:624121/1‑251 (MQ=255)
                                                                                                                  aGTATTACAACGACAAAGTCCCGCACAACCGCGGCTGTGTTTATCAAGATGTGCTAAATATGCTGCTGTTATACTCCCGTCCTAAGGCGGCGTACTGCACGGGAAAGTGGATAAATGCCCCCGGCGGTCAGGGGTTGTTGT‑GATCAACAGGG‑GGGA‑CCCATTAGGGATTTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCTTTTAAGGCAAACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATgg                                                                                                                                      <  2:188956/251‑1 (MQ=255)
                                                                                                                  aGTAATACATCGACAAAGTACCGCTCAACCGCGGCTGCTATTATCAAGATGTTCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCCCCGTACAGTCGATCAATTTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑CCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCTCTTCGCTCACCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATgg                                                                                                                                      <  1:88349/251‑1 (MQ=255)
                                                                                                                   gTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCTGGGGATTGTCT‑GATCCACAGCG‑GTGA‑TCCATTAAGGATCTCCCTGGGAACTA‑AGAAAGACTGTGGCGCATCGCTCCTCCGGCATACCCG‑GCTTCCCCCCCCACAT‑TTAGGAAAAGTTATGGt                                                                                                                                     >  2:81546/1‑251 (MQ=255)
                                                                                                                        aCATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGGTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGAGATGATCTTGA‑CAACAACG‑G‑GAGTACATTTAGGATCTCCGTGAGA‑CTATAGAAAGCCAGAAGCGCTACGCTCCTCAGGCGTACAAG‑ACTTCCGCCACTACAT‑TTAGGAAAAGTTATGGTAtacc                                                                                                                                >  1:451741/1‑248 (MQ=255)
                                                                                                                                 aaGTACCGCCCAAACGCGGCGGGGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCATACCTAAGGCGTCGGACTGCACCGGACAGTCAATAAAGGTCACCGGCGGTCAGGTGATGTTCG‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCTCATCAGGCATAAAGG‑TCTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACtcatca                                                                                                                      <  1:811932/251‑1 (MQ=255)
                                                                                                                                      ccGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGCCCTAAGGGTTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCGGTGAGA‑CTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTc                                                                                                                 <  1:834603/251‑1 (MQ=255)
                                                                                                                                            aaaCGCGGCTGCGATTCTCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCT‑GATCAACAGAG‑GAGA‑TCCCTTAAGGATCTCCGTGAGA‑CTATAAAATGCCTGAAGCGCTACGCTCATCAGGACTACAGG‑ACTTCCGCCACTACAT‑TAAGGGAAAGTTATTGTATCC‑GCACTCATCCCCgctgcccat                                                                                                           >  1:724549/1‑243 (MQ=255)
                                                                                                                                               cgcgGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATTTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATc                                                                                                        <  1:770180/251‑1 (MQ=255)
                                                                                                                                                gcgGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCT‑GTTAAACACCG‑GGGA‑TCCCATAAGGGTCTCCCTGAGC‑CTATAAAATGACTCCTGCGCTCCGCTCATCAAGCGTACAAG‑GCTTTCGCCCCTA‑ATATTAGGGAAAGTTTTGGTATCC‑GCACtactcacccccccccttatac                                                                                                       >  2:784000/1‑231 (MQ=255)
                                                                                                                                                    cTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCT‑GATCCACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTg                                                                                                   >  1:590312/1‑251 (MQ=255)
                                                                                                                                                     tGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCT‑GATCAACAACG‑GAGA‑TCCATTAAAGATCTCCGTGAGA‑CTATAAAATGCCTGATGCGCTACGCTCATCAGGCATACAAG‑ACTTCCGCCACTACAT‑TTAGGAAAAATTATGGTATTA‑GCACTCATCACCGTCGCCGTTATCGCCTgg                                                                                                  >  1:289451/1‑251 (MQ=255)
                                                                                                                                                         aTTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGCCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGG‑CCTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGtgtg                                                                                              <  1:40349/251‑1 (MQ=255)
                                                                                                                                                          ttATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGAGATGTGCT‑GAACAACAGCG‑GAGG‑TCCATACAAGATCTCCGTGAGA‑CTATAGGATGCCTGCGGCGCTACGCACAACAAGCATACCGG‑ACCTCCGCCCCCCCAT‑TAAGGGACAGTTAATGAAACC‑GCACCCATCACCGGCGCCCTATTCGCCTGgtgagc                                                                                             >  2:457952/1‑248 (MQ=255)
                                                                                                                                                                   tgtgTTGAATAGGCTGCGGGTCTACGCCAGTCCTAGGGAGTCGTACTGCCTCGTACAGGCGATCAATGTCACCGGCGTTCAGTTGATTTTCT‑GAATAACGGCG‑GAGA‑TCCATAAAGGATCTCCGTGAGA‑ATATAGATTGCCTGTTGCGCTACGCTAATTAGGAATACAGG‑ATTGCCGCCACTACAT‑TAATGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTgg                                                                                     <  1:762227/250‑1 (MQ=255)
                                                                                                                                                                         gAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCT‑GATCAACCGCG‑GAGA‑TCCATTAAGGATCGCCGCGGGA‑CTATAGAAAGCCCGCTGCGCTCCGCCCACCAGGCCTAAAGG‑GCTTCCGCCACCAACT‑TAAGGGAAAGTTTTGGTATCC‑CCACTCATCACCCGCGCCGTTATTGCCTgggggtggcaaatggcccttgg                                                                              >  2:670771/1‑231 (MQ=255)
                                                                                                                                                                             atGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTTCGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTTGCCTTAagcggc                                                                          >  1:521800/1‑251 (MQ=255)
                                                                                                                                                                             atGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAAAGCCTGATGCGCTACGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggcggc                                                                          >  1:6908/1‑251 (MQ=255)
                                                                                                                                                                             atGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCT‑GATCAACAGCG‑GAGA‑AACCTTAAGGGTCTCCGTGAGA‑ATATAGAATGCCTGATGCGCCACGCTCCTACCGCATACAGG‑CATTCCGCCACTACAT‑TTAGGAAAAGTTAAGGAATCC‑GCCCTCATCACCGTCCCCCTTCTCGCCTGGTGTGCGCAACtcgcattcagcggc                                                                          >  2:295104/1‑238 (MQ=255)
                                                                                                                                                                                      ttCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GGTCAACAGCG‑GAGA‑TCCATTTAGGATCTCCGGGGGA‑CTATAGAATGCCTGATGCGCTACGCCCATCCGGCATACAGG‑GCTTTCGCCCCCACAT‑TTAGGCAAAGTTTTGGTATCC‑CCCCTCCTCACCCTCCCCGTTTTCGCCTGGTGTTGGCCACCGGCCTTTGGCCGCTGGGCAAAt                                                                 >  2:639804/1‑251 (MQ=255)
                                                                                                                                                                                              cggtccaagggcttcgtACTCCGCCGGACAGACGATCAATGTACCCGGCGGTCGTGTAATTTTTT‑GACAAACAGCG‑GACA‑TATATTAAGGTTCTCAGTGAGA‑CTATAGAATGCCTTTTGCTCTACTCTCACCAGGCATACAGC‑CCTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGG‑CAAATTTCTCGttt                                                        <  2:682458/238‑1 (MQ=255)
                                                                                                                                                                                                           gcgtACTGCACCGGAAAGTCGATCAATGTCACGGGCGGTCAGGGGATGTTCA‑GATCACCTGCG‑GAAA‑TCCATTAAGGATCTCGGTGAGT‑CTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGG‑CAAATTTCTCGTTTTAACCGTGCCtgc                                           <  1:719929/250‑3 (MQ=255)
                                                                                                                                                                                                               aCTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGGGGAGTTCC‑GGTCCACCGCG‑GAGT‑TCCCTTTAGGGTCTCCGTAAGA‑ATATAGAATGCCTGGTGCGCTTCGCCCATCTGGCGTACAGG‑ACTTCCGCCACCACAT‑TAAGGAAAAGTTATTGTATCCTGC‑CTCATCACCGCCCCCGTTCTCGCCTTGTTTGCCCAACCGGCCTTAGGGGGCTGG‑CGAATTTCTCGGTTTTACTGTGCCTTccacc                                       >  2:12668/1‑247 (MQ=255)
                                                                                                                                                                                                                  ccACCGAACAGGCGATCAACGTCTCCGGCGTGCAGGTGCTGTTCT‑AACCAACAGCG‑TAGA‑TCAATAAAGGAGCTTGGTGAGA‑TGATAAATTCCCGGATGCGTTCCCCGCAATAGGCATACAGC‑TCTTCCGCCACGACAT‑TAATGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGTTGTGCGCATCTGGCCTGAGGCGGCTGG‑CAAATTTCTCGTTTTAACCGTGCCTTCGacacac                                    <  1:82029/250‑1 (MQ=255)
                                                                                                                                                                                                                      cGGACAGTCGATCAATGTCACCGGCGGTCAGGTGCTGTTCT‑GATCAACCGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATATAAAGCCCGATGTGCTAAGCACATCAGGCATACAGG‑ACTTCCGCCCCTCCAT‑TAAGGAAAAGATATGGTATCC‑TCACTCATCCCCCACGCCCTTATCGCATGGGGTTCGCACCCGGCCTTAAGGTGCTGG‑CCAATTTCTCGtttttaacccgcccttgacacaccatt                                >  2:18686/1‑228 (MQ=255)
                                                                                                                                                                                                                       ggACAGTCGATCAATGTCACCGGCGGTCAGGAGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑ATATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TTAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGG‑CAAATTTCTCGTTTTAAACGTGCCTTCGACACACttagc                               >  1:294825/1‑247 (MQ=255)
                                                                                                                                                                                                                          aaGTCGAACAATGCCACGGGCGGTCGGGTGTTTTCCT‑GATCAACA‑CGTGAGA‑TCCATTAAGGGTCCCGGGGAGC‑CTATAGAATCCCTGATCCGCTCCCCTCATCAGGCAAACGGC‑CCTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGAGGCTGG‑CAAATTTCACGTTTTAACCGTGCCTTCGACACACTATGCCAg                            <  1:234312/250‑1 (MQ=255)
                                                                                                                                                                                                                               gATCAATGTCACCGGCGGTCAGGGGATTTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTTTAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATGGGATCC‑GCACTCATCACCGTCGCCGTTTTCGCCTGGGGTGCGCAACTGGCCTTAGGCGGCCGG‑CAAATTTCTCGTTTTAACCGTGCCTTTGACACAATATTCCAGCAAgg                       >  1:862009/1‑251 (MQ=255)
                                                                                                                                                                                                                                aTCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCCACGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATGGGATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGGGCGCAACTGGCCTTAGGCGGCTGG‑CAAATTTCTCGTTTTAACCGTGCCTTCGACACACTATGCCAGCAAggg                      >  1:694710/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGG‑ACTTCCGCCCCTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGG‑CAAATTTCTCGTTTTAAACGTGCCTTCGACACACTATGCCAGCAAGGGCg                    >  1:817111/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cACCGGCGGTAAGGTGATGTTCT‑GATCAACGGCG‑GAGA‑CCCATTATGTATCCCCGTGAGA‑CTATAGAATGCCTGATCCGCTACGCACATCAGGCATACTGG‑TCTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGG‑CAAATTTCTCGTTTTAACCGTGCCTTCGACACACTATGCCAGCAAGGGCGGGTTg               <  1:184417/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 tCAGGTGATGTTCT‑AATCAACAGCG‑GGGA‑TCCATTAAGGACCTCCGTGAGA‑TTAAAGAAGGCCGGATGCGCTCCGCTCATCAGGCATACAGG‑CCTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGG‑CAAATTTCTCGTTTTAACCGTGCCTTCGACACACTATGCCAGCAAGGGCGGGTTGGCGTGGGcc      <  1:644679/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    ggTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGG‑CAAATTTCTCGTTTTAACCGTGCCTTCGACACACTATGCCAGCAAGGGCGGGTTGGCGTGGGCCGTTc  >  1:312452/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                
ACGGCATTACGGTGCATTCACTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAACGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCT‑GATCAACAGCG‑GAGA‑TCCATTAAGGATCTCCGTGAGA‑CTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGG‑ACTTCCGCCACTACAT‑TAAGGAAAAGTTATGGTATCC‑GCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGG‑CAAATTTCTCGTTTTAACCGTGCCTTCGACACACTATGCCAGCAAGGGCGGGTTGGCGTGGGCCGTTC  >  NC_000913/2828260‑2828754

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