Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,488,2090GT25.0% 122.3 / 10.6 57G298C (GGC→TGC) lptGlipopolysaccharide transport system protein LptG
Reads supporting (aligned to +/- strand):  ref base G (12/30);  new base T (14/0);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.66e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATCCGAAACAGATTACCGGTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTG  >  NC_000913/4487964‑4488456
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                              
aTCCGAAACAGATTACCGGTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGAGGATGCCGATGGCGGTGTCGTTCATCTTT‑TGcccc                                                                                                                                                                                                                                                           >  1:189639/1‑250 (MQ=255)
    gAAACAGATTACCGGTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTgc                                                                                                                                                                                                                                                       <  1:441091/251‑1 (MQ=255)
              aCCGGTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTAcc                                                                                                                                                                                                                                             <  1:524031/251‑1 (MQ=255)
                 ggtgcggagcaggTGAGCGGCATCGGGAAAGCCAGCCTGACGCCGGACATCCTGGGCGCGGTGGCGCTGGAGCCGGATGCACTCTTTATCAGCGGTTTGCACAATTATGAGAAGTATCTGAAGTCGAGGGGACAGGATGCCCGACCTTCGCAGCTCAAAACGTGGAGGAAAATCTCCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GTCCCACTGCGTAGGGTACCGag                                                                                                                                                                                                                                          <  2:61543/240‑2 (MQ=255)
                 ggTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGAt                                                                                                                                                                                                                                          <  1:267129/251‑1 (MQ=255)
                 ggTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGAt                                                                                                                                                                                                                                          <  1:533782/251‑1 (MQ=255)
                  gTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGGAGATGTTGATGGCGCTTGCGTTCATATTT‑GGGCCACTGCGTAGCGTAcagc                                                                                                                                                                                                                                           >  1:587509/1‑246 (MQ=255)
                       cAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGGGGTGGTGCCTTTGGCGCTGTCTTTCCTCTTT‑TGCCCCCTgtgttgcgtcacggtggggtg                                                                                                                                                                                                                                    >  2:116114/1‑239 (MQ=255)
                                   ggCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGAGGATGCTGATGGCGCTGTCGGTCATCTTT‑GGCCCCCTGCGTAGCGGACCGATGGGCGTGCGTGTGG‑tgac                                                                                                                                                                                                                       >  1:308245/1‑248 (MQ=255)
                                     cACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTTATGATGCTGATGGCGCTGTCGTTCATCTTT‑TGCCCACCGCGTTGCGTACCGAAGGGCGGGCGTGTGG‑GCACCg                                                                                                                                                                                                                     >  1:149330/1‑251 (MQ=255)
                                             aaaCCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGGGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTGT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑Gt                                                                                                                                                                                                            <  2:117641/251‑1 (MQ=255)
                                                 caaccTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCg                                                                                                                                                                                                        <  1:236268/251‑1 (MQ=255)
                                                        ccGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGt                                                                                                                                                                                                 <  2:189639/250‑1 (MQ=255)
                                                            cGGACAAACTGGGCGTGGTGGCTCTGGACCCGGATGCATTCTCATTCAGCGGTTTGCACATCTATGTAACGTCGCTGAAGTCGACCGGTCAGGATGCCGGACGTTAGCAGCTCAACATGTGGACCAAAATCTTCCAGACGCTATCTGTGGCGGTGATGATTCTGATGGCTCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑GCACCGGTCTCA‑GTTTCGGTTTTGtcttc                                                                                                                                                                                             <  2:319148/251‑1 (MQ=255)
                                                                      gggCGTGGTGGCGCTGGCCCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACt                                                                                                                                                                                    <  2:446683/249‑1 (MQ=255)
                                                                            ggtgtcccttgACCCGGATGCCCTCTCTACCAGCGGTTTGACCTACTAGGTGAAGTATATGAAGTCAAGCGGTCAGGATGCCGGACCTTATCAGCTCAACATTTGTAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCTACTGCGTATCGTACCGATGGGCTTGTGTGTGG‑TCACAGGTATCA‑GTTTCGGTTTTTTCGTCT‑ACGTACTGGACCAGa                                                                                                                                                                            <  2:587509/241‑1 (MQ=255)
                                                                                 cgcTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGGGCCGTCGTTCATCTTT‑TGCCCACCTCGTAGCGGACCGGTGGGCGGGCGGGTTG‑GCACCGGGATCC‑GTTTTGGGTTTGtcttttcctgactgggcccaaaattt                                                                                                                                                                        >  2:558407/1‑229 (MQ=255)
                                                                                  gcTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGGTGATGGTGATGGCGCTGTCCTTTATCTTT‑TGCCCCCCGCGGAGCGTACCGGTGGGCGGGCG‑GTGGGTCACCCGGATTC‑GTTTCGGGTTTTTCTTCT‑TCGTACttggcccagtcctct                                                                                                                                                                      >  2:431583/1‑237 (MQ=255)
                                                                                               gCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCCTCTTT‑GGGCCACTGCGTAGCGTACCGATGGGCGGGCGTGTGG‑TCACCGGTATCA‑Gttt                                                                                                                                                                                                          >  1:5631/1‑203 (MQ=255)
                                                                                               gCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑Gttt                                                                                                                                                                                                          <  2:5631/203‑1 (MQ=255)
                                                                                                 aCTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGtcttct                                                                                                                                                                                            <  1:316399/215‑1 (MQ=255)
                                                                                                      ctATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGGTGATGCTGATGGCGCTGTCGTTCATCTTT‑TGCCCACTGCGTAGCGTACCGATGGGCGGGCGTGTGG‑GCACCGGTATCA‑GGTTCGGGTTTTTCTTCT‑ACCTACTGGACCAGGTCTTCGGCCCGCT‑GACGTTGGttttt                                                                                                                                                 >  1:327217/1‑249 (MQ=255)
                                                                                                        atcatcGGTTTGCACAACTATGGGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATgg                                                                                                                                               <  2:308245/251‑1 (MQ=255)
                                                                                                               gTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑TGCCCACTGCGTAGGGTACCGCTGGGCGGGCGTGTGG‑TCACCGGGATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTTTTGCATcccc                                                                                                                                        >  1:409801/1‑250 (MQ=255)
                                                                                                                   gCACAACTATGTTAAGTATCGGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCcc‑gccg                                                                                                                                    <  1:480084/250‑1 (MQ=255)
                                                                                                                       aaCTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGCAGCGTACCGATGGGCGTGCGTGTGG‑GCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATTTTCGGCCCGCT‑GACGTTGGTTTATTGCATCcc‑gcccctcat                                                                                                                               >  1:166999/1‑245 (MQ=255)
                                                                                                                         cTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGTAGATGGTGATTGCGCCGTCGTTCATCTTT‑TGCCCACTGCGTAGCGGACCCATGGGCGGGCGGGGGG‑TCCCCGGGATCC‑GTTTCGGTTTTTTCTTCTCAC‑TACTGGGCCCAATCTTCGGCCCGCT‑GTCGTTTGTTTATGGGATCcc‑gcccaaccttc                                                                                                                             >  2:540928/1‑243 (MQ=255)
                                                                                                                             gtgAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGCCGTTTATCTTT‑GGCCCACTGCGTAGCGCCCCGATGGGCGTGCGTGTGG‑TCACCGGGATCA‑GTTTCGGGCTTGTCTTCT‑TCGTCCTGGATCAGATCTTCGGCCCGCT‑GACGCTGGGTTTTGGCATCCC‑GCCGGTCACCGgcgc                                                                                                                         >  1:344503/1‑251 (MQ=255)
                                                                                                                                aaGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCCTCTTT‑GGCCCCCTGCGTAGCGTACCGATGGGCGGGCGTGTGG‑GCCCCGGTATCA‑GTTTCGGTTTTTTTTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GGCGTTGGTTTTTGGCATCCC‑GCCGATCAACGGCGCAcc                                                                                                                      >  2:295558/1‑250 (MQ=255)
                                                                                                                                         aaGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGGTGGCGCTGTCGTTCCTCTTT‑TGCCCACTGCGTAGCGTACCGATGGGCGGGCGTGGGG‑TCACCGGGATCC‑GTTTCGGGTTTGTCTTCT‑TCGGAATGGACCAGATCTTCGGCCCCCT‑GTCGTTGGGTTTTTGCATCCC‑CCCGCTCATCCGCGCACTG‑TTG‑CCa                                                                                                             >  2:336210/1‑249 (MQ=255)
                                                                                                                                             cGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAgcgc                                                                                                        <  1:246144/250‑1 (MQ=255)
                                                                                                                                             cGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGGGGCGGTGATGGTGCTGATGGCGCTGTTGTTCCTCTTT‑CGCGCACTGCGGAGCGTATCGAAGGGCGCGCGGGTGG‑TCGGCGGGATCA‑GTTTCGGGTTTGTCTTCT‑TAGGACAGGGCCCGAACTTCGGCCCGCTTGA‑GGTGGTTTTTGGCCTACC‑GCCGGTCATCTgcgcgccatttccacgcgcc                                                                                                         >  2:10358/1‑235 (MQ=255)
                                                                                                                                               aGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TCG‑CCAAGCGCCAg                                                                                                     <  1:431583/251‑1 (MQ=255)
                                                                                                                                               aGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGAAGATGCTGGTGGAGCTGTCGTTCATCTTT‑TGCCCACCGCGTAGGGTAACGATGGGGGTGCGTGTGG‑TCAACGGCATCA‑GTTTCGGTTTTGTCTTTT‑ACGTACTGGACCAGATCTTCGGCCCGCC‑TACCTTTGTTTTTGGGATCCC‑GCCCATCCACGGCGCACTG‑TTG‑CCAAACGcccc                                                                                                     >  2:437683/1‑249 (MQ=255)
                                                                                                                                                   gtcagtaggCCGGACGTTTTCAGCTCGACATGTGGAGCAAAATCTTCCGGCCGCTATCTGTGGCGGGGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGGGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTAGGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGcttc                                                                                                 <  2:344503/243‑1 (MQ=255)
                                                                                                                                                       ggATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGGTGGCGCCGTCGTTCATCTTT‑TGCCCACTGCGTAGCGGACCGATGGGCGTGCGTGGGG‑GCCCCGGGATCC‑GTTTCGGTTTTTTCTTTT‑ACGGACTGGACCAGATCTTCCGCCCCCC‑GGCGTTGGTTTTTGGCATCCC‑CCCCATCATCGGCGCACCG‑TTG‑CCTAGCGCCCGcgtcttct                                                                                             >  2:121353/1‑251 (MQ=255)
                                                                                                                                                        gATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCCCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCCCTG‑TTG‑CCCAGCGCCAGcctcttctt                                                                                            >  1:53139/1‑251 (MQ=255)
                                                                                                                                                            ccGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATc                                                                                        <  1:169166/251‑1 (MQ=255)
                                                                                                                                                            ccGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATc                                                                                        <  1:314730/251‑1 (MQ=255)
                                                                                                                                                                  nttATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGACGATGCTGGTGGAGCTCTCCTTCCACTTT‑TGCCCCCCGCGTAGCGCAACGCTGGGCGTTCGTGGGG‑TCACCGGTATCC‑GTTTCGGTTTTGTTTTCT‑ACGCACTGGGCCAGACCTTCGGCCCGCT‑GACGTTGGTTTATTGCCTCCC‑GCCCACCATCGCCGCACTC‑TTGACCAA‑CGCCCGCGTCTTCTTactccaccct                                                                                  >  2:276279/2‑242 (MQ=255)
                                                                                                                                                                       cAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTCCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCt                                                                             <  1:128534/251‑1 (MQ=255)
                                                                                                                                                                        aGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTg                                                                            <  1:295558/251‑1 (MQ=255)
                                                                                                                                                                          cTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑TGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCCGTATCA‑GTTTCGGGTTTTTCTTCT‑ACGTACTGGGCCAGATCTTCGGCCCGCC‑GACGTTGGGTTTTTGCATCCC‑CCCGATCATCGGCGCCCTG‑TTG‑CCAAGCGCCAGCTTCTTCTTTATCAGCCCGTGGCTggt                                                                          >  1:166848/1‑249 (MQ=255)
                                                                                                                                                                           tCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGAAGGCGCTGTCGTTCATCTTT‑TGCCCACCGCGTAGCGTACCGATGGGGGGGCGTGTGG‑TTACCGGGATCC‑GTTTCGGGTTTGTCTTCT‑ACGGTCTGGGCCCGATTTTCGGCCCGCT‑GACGTTGGTTTTTGGCATCCC‑GCCGCTCATCCGCGCACTG‑TTG‑CCAAGGGCCAGCTTCTTCTTAATCAGCCCGTGGCCGttt                                                                         >  2:349701/1‑250 (MQ=255)
                                                                                                                                                                                                     ccGCTATCTGTGGCGGTGATGATGCTGATGGCGCTTTTGGTCTTCTTTTGGCCC‑CTGCGGAGCGGAACGATGGGCGTGCGTGTGG‑TCACCGGTATCAAGTTT‑GGGTTTGTCTTTT‑AACTACCGGACCAGATCTTTGGCCCGCT‑GACGTTTGTTTTTGGCATCCCCGCC‑CTCATAGGCGCACTGGTTG‑CCA‑GCGCCAGCTTCTTCTTAATCAGCCTGGGGCTGGTAAATAGAAAATccgaaaccaaagaa                                                 >  2:385534/1‑236 (MQ=255)
                                                                                                                                                                                                       gCTATCTGTGGTGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGa                                              >  1:95820/1‑250 (MQ=255)
                                                                                                                                                                                                                ggCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGGGTAGCGTATCGATGGGCGTGCTTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGCCCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGTCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGa                                     <  1:10358/250‑1 (MQ=255)
                                                                                                                                                                                                                ggCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAggg                                      >  2:453844/1‑249 (MQ=255)
                                                                                                                                                                                                                 gCGGTGATGATGCTGATGGCGCTGTCGTCCATCTTT‑GCCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGTTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGCCCAGATCTCCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑ACAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGc                                   <  1:483657/251‑1 (MQ=255)
                                                                                                                                                                                                                    gTGATGATGCTGATGGCGCTGTCGTTCTTCTTT‑GGCCCACTGCGTGGCGTACCGATGGGCGTGCGTGGGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCTAATGAACGAATAAAGGGAGCGaa                                <  2:218937/251‑1 (MQ=255)
                                                                                                                                                                                                                       atgatgCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGc                             >  1:73006/1‑251 (MQ=255)
                                                                                                                                                                                                                         gatgCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCt                            <  1:549047/250‑1 (MQ=255)
                                                                                                                                                                                                                           tgCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTCCCGATTGGCGTGCGTGTGG‑TCCCCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTcc                          <  1:187456/250‑1 (MQ=255)
                                                                                                                                                                                                                            gCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTccc                         <  1:267980/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gTCGTTCATCTTT‑GGCCCCCTGCGTAGCGTACCGTTGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGtt             <  1:116114/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ttCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCa         >  2:110354/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   tt‑GGCCCACTGCGTAGCGTACCGATGGGCGGGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGtgtg  <  2:53139/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                              
ATCCGAAACAGATTACCGGTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTT‑GGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGG‑TCACCGGTATCA‑GTTTCGGTTTTGTCTTCT‑ACGTACTGGACCAGATCTTCGGCCCGCT‑GACGTTGGTTTATGGCATCCC‑GCCGATCATCGGCGCACTG‑TTG‑CCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTG  >  NC_000913/4487964‑4488456

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