Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,248,4240CA25.6% 73.5 / 12.9 39V37F (GTT→TTT) lysPlysine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (25/4);  new base A (0/10);  total (25/14)
Fisher's exact test for biased strand distribution p-value = 1.57e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTAAAGCAGCCAAT  >  NC_000913/2248174‑2248659
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
aCCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGcaac                                                                                                                                                                                                                                               <  1:42456/251‑1 (MQ=255)
 ccAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGcaaca                                                                                                                                                                                                                                              <  1:323264/251‑1 (MQ=255)
         gctgCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGAcc                                                                                                                                                                                                                                      >  2:312300/1‑251 (MQ=255)
            gCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCt                                                                                                                                                                                                                                     <  1:336070/249‑1 (MQ=255)
               aCCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTAcc                                                                                                                                                                                                                                <  1:240596/251‑1 (MQ=255)
                       aaCGGCGATAGTCCCCGCCCAGTTGTACCAGTAGTTCCAGCCCAGAGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑c                                                                                                                                                                                                                      >  2:300236/1‑251 (MQ=255)
                           gCGAGAGTCACCGCCCATTTGTCCCAGTAGTCCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGTCCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGCCCTGTTCCGA‑TGGAA‑CCGcc                                                                                                                                                                                                                  <  2:53399/251‑1 (MQ=255)
                                       gCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCACCAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAatcat                                                                                                                                                                                                      >  1:329687/1‑251 (MQ=255)
                                           aGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTc                                                                                                                                                                                                  >  1:44150/1‑251 (MQ=255)
                                                             agcccagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAt                                                                                                                                                                                >  2:41920/1‑251 (MQ=255)
                                                                  agcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACAATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CAGCCAATGGCAATCATCGTCAGGTGAAGCGCCTTTAATTCACg                                                                                                                                                                           >  2:148876/1‑251 (MQ=255)
                                                                   gggcaaaGCAAAAGCTTTCTTAAACATGGTTCTGACCGAAAGGGGAAAACGACCGGGAAACA‑GGCATAAATGCAGCCAGTTCCCCGAGACTGGTAATAAGAAAGTAACCCATCAGGCCAATCAGCATATACGAGAGAAAGCCCCCCCCC‑GGCCCTGCGTGAAAAAT‑C‑GTTGGGACAGAG‑CAAAAAAAAGCCCTGTCCCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGc                                                                                                                                                                         <  1:52876/246‑1 (MQ=255)
                                                                                     cttcAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCt                                                                                                                                                        >  1:94578/1‑251 (MQ=255)
                                                                                                          aaGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAcc                                                                                                                                   >  1:191837/1‑251 (MQ=255)
                                                                                                            gggccAAACGA‑CCGGAAACCCGGAATATATGCGGCCAGTTCACCGAGCCTGGTCATCAGGAAGTAAACCATCAGCCCAATCAGAATATACGAGAGCAATCCCCCGCCC‑GACCCTGCCTGAAAAAT‑C‑GTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGA‑TGAAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCAt                                                                                                                                 <  1:27603/247‑1 (MQ=255)
                                                                                                            gTGGCAAACGACCCGGAAACC‑GTAATATAGGCAGCCATTTCCCCGAGACTGTTACTCAGGAAGTAAACCATCAGCCAAATCAGAATATACGAGAGCAGGGCCCCCCCC‑GGACCGGCATGAGAAAT‑T‑GTGGCCAGAGAGGAAACAAAAAGACCTGTACCTA‑AGAAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCAt                                                                                                                                 <  1:260575/251‑1 (MQ=255)
                                                                                                                     gAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                                                                                                                        >  1:83091/1‑251 (MQ=255)
                                                                                                                      aaCCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGc                                                                                                                       >  1:13805/1‑251 (MQ=255)
                                                                                                                         ggaaaaCC‑GGAATATAGGCAGCAAGCACACCGAGGCGGGTAATAAGGAAGAAACCCATCGGGCCAATCAGCATACACGAGAGCAACCCCCCGCCC‑GCACCTGCCTGAGAAAG‑C‑GTGGCGCGAGAGGCAAAAAAAAGCCCTGTCCCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTAt                                                                                                                    <  2:44150/248‑1 (MQ=255)
                                                                                                                            aaaCC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGACCAGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAACGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGGGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGGTTTAGTTTCGGAAACCATACAAAAATGCTATcc                                                                                                                  >  1:416461/1‑250 (MQ=255)
                                                                                                                                        tatGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGTAAACCATACAAAAATGCTATCCATCTTCAAAATg                                                                                                      >  2:8788/1‑251 (MQ=255)
                                                                                                                                                    tCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGAtttt                                                                                          >  2:98598/1‑249 (MQ=255)
                                                                                                                                                     cACCGAGCCTGGTCAACAGAAAGTAACCCACCAGCCAAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGCCCTGCCTGAAAAAT‑C‑GTTGCGCCAGAGGCAAAAAAAAGCCCTGTACCGA‑TGAAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTa                                                                                         <  1:126938/251‑1 (MQ=255)
                                                                                                                                                      aCCGAGACTGGTCATCAGGAAGTAAACAATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGt                                                                                          >  1:14514/1‑249 (MQ=255)
                                                                                                                                                                gTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAACAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGCTTGTAGCAGACGATTg                                                                              >  2:412711/1‑251 (MQ=255)
                                                                                                                                                                gTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTg                                                                              >  1:427020/1‑251 (MQ=255)
                                                                                                                                                                     caGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAAt                                                                         >  1:72020/1‑251 (MQ=255)
                                                                                                                                                                      aGGAAGTAACCCAACAGCCAAATCAGAATATACGAGAGCAACGCCCCCCCC‑GGACCTGCCGGAAAAAT‑C‑GTTGCGCCAGAGGAAAAAAAAAGCCCGGTCCCGA‑GGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATc                                                                        <  1:312300/251‑1 (MQ=255)
                                                                                                                                                                        gAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGCGGCGCGATTGTAGCAGACGATTGGCAATCCt                                                                      >  2:97898/1‑251 (MQ=255)
                                                                                                                                                                         aaGTAAACCATCAGGCTAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAG‑A‑GTTGCGCCAGAGGCACCAAAACGACCAGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGAATGTAGCAGACGCTTGGAAATCCtt                                                                     >  2:232810/1‑251 (MQ=255)
                                                                                                                                                                            aaaacaaccagacaactcaccaaataCGAGAAAAACGCCCCGCCC‑GG‑CCTGGCTTAAAAATTC‑GTTGCGCCGAAGCCAAAAAAAATCCCTGTCCCG‑‑TGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCg                                                                  <  1:28340/228‑1 (MQ=255)
                                                                                                                                                                             aaaCCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCgg                                                                 >  2:228102/1‑251 (MQ=255)
                                                                                                                                                                                 cATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGAGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTACAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGc                                                                >  2:33657/1‑248 (MQ=255)
                                                                                                                                                                                        ccagaaagagaataCGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAG‑G‑GGTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGc                                                                                                                       >  2:47114/12‑186 (MQ=255)
                                                                                                                                                                                                 atatACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTtata                                               >  1:86823/1‑249 (MQ=255)
                                                                                                                                                                                                    acggagaGCA‑TGCCCCGCCCAGGACTTGCC‑GAAAAAT‑CCGTTGCGCAGGAGCCAAAAAAAAG‑CCTTGACCGAATGAAT‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGaca                                          <  1:133627/248‑1 (MQ=255)
                                                                                                                                                                                                     aCGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGacac                                         >  1:265205/1‑251 (MQ=255)
                                                                                                                                                                                                       ggagagaaatcccccccccgaccccgcctgAAAAAG‑C‑GTGGCGACAGAGGAAAAAAAAAGCCCTGGACAGA‑GGAAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACAcc                                        <  1:97898/226‑1 (MQ=255)
                                                                                                                                                                                                              aTGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGAt                                >  1:24962/1‑251 (MQ=255)
                                                                                                                                                                                                                         cggccctgcctgAAAAAT‑C‑GTTGCGCCAGAGGCAAAAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGa                      <  1:369567/247‑1 (MQ=255)
                                                                                                                                                                                                                                                tGCGCCAGAGACAAAAAAAAGCCCTGTCCCGA‑GGAAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTAAAGCAGCCAAt  <  1:228102/250‑1 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
ACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACC‑GGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCC‑GGACCTGCCTGAGAAAT‑C‑GTTGCGCCAGAGGCAACAAAAAGACCTGTACCGA‑TGGAA‑CCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTAAAGCAGCCAAT  >  NC_000913/2248174‑2248659

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