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.

ATCCGAAACAGATTACCGGTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTG  >  NC_000913/4487964‑4488456
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                              
aTCCGAAACAGATTACCGGTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGAGGATGCCGATGGCGGTGTCGTTCATCTTTTGcccc                                                                                                                                                                                                                                                           >  1:189639/1‑250 (MQ=255)
    gAAACAGATTACCGGTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTgc                                                                                                                                                                                                                                                       <  1:441091/251‑1 (MQ=255)
              aCCGGTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTAcc                                                                                                                                                                                                                                             <  1:524031/251‑1 (MQ=255)
                 ggtgcggagcaggTGAGCGGCATCGGGAAAGCCAGCCTGACGCCGGACATCCTGGGCGCGGTGGCGCTGGAGCCGGATGCACTCTTTATCAGCGGTTTGCACAATTATGAGAAGTATCTGAAGTCGAGGGGACAGGATGCCCGACCTTCGCAGCTCAAAACGTGGAGGAAAATCTCCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGTCCCACTGCGTAGGGTACCGag                                                                                                                                                                                                                                          <  2:61543/240‑2 (MQ=255)
                 ggTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGAt                                                                                                                                                                                                                                          <  1:267129/251‑1 (MQ=255)
                 ggTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGAt                                                                                                                                                                                                                                          <  1:533782/251‑1 (MQ=255)
                  gTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGGAGATGTTGATGGCGCTTGCGTTCATATTTGGGCCACTGCGTAGCGTAcagc                                                                                                                                                                                                                                           >  1:587509/1‑246 (MQ=255)
                       cAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGGGGTGGTGCCTTTGGCGCTGTCTTTCCTCTTTTGCCCCCTgtgttgcgtcacggtggggtg                                                                                                                                                                                                                                    >  2:116114/1‑239 (MQ=255)
                                   ggCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGAGGATGCTGATGGCGCTGTCGGTCATCTTTGGCCCCCTGCGTAGCGGACCGATGGGCGTGCGTGTGGtgac                                                                                                                                                                                                                       >  1:308245/1‑248 (MQ=255)
                                     cACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTTATGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCACCGCGTTGCGTACCGAAGGGCGGGCGTGTGGGCACCg                                                                                                                                                                                                                     >  1:149330/1‑251 (MQ=255)
                                             aaaCCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGGGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTGTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGt                                                                                                                                                                                                            <  2:117641/251‑1 (MQ=255)
                                                 caaccTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCg                                                                                                                                                                                                        <  1:236268/251‑1 (MQ=255)
                                                        ccGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGt                                                                                                                                                                                                 <  2:189639/250‑1 (MQ=255)
                                                            cGGACAAACTGGGCGTGGTGGCTCTGGACCCGGATGCATTCTCATTCAGCGGTTTGCACATCTATGTAACGTCGCTGAAGTCGACCGGTCAGGATGCCGGACGTTAGCAGCTCAACATGTGGACCAAAATCTTCCAGACGCTATCTGTGGCGGTGATGATTCTGATGGCTCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGGCACCGGTCTCAGTTTCGGTTTTGtcttc                                                                                                                                                                                             <  2:319148/251‑1 (MQ=255)
                                                                      gggCGTGGTGGCGCTGGCCCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACt                                                                                                                                                                                    <  2:446683/249‑1 (MQ=255)
                                                                            ggtgtcccttgACCCGGATGCCCTCTCTACCAGCGGTTTGACCTACTAGGTGAAGTATATGAAGTCAAGCGGTCAGGATGCCGGACCTTATCAGCTCAACATTTGTAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCTACTGCGTATCGTACCGATGGGCTTGTGTGTGGTCACAGGTATCAGTTTCGGTTTTTTCGTCTACGTACTGGACCAGa                                                                                                                                                                            <  2:587509/241‑1 (MQ=255)
                                                                                 cgcTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGGGCCGTCGTTCATCTTTTGCCCACCTCGTAGCGGACCGGTGGGCGGGCGGGTTGGCACCGGGATCCGTTTTGGGTTTGtcttttcctgactgggcccaaaattt                                                                                                                                                                        >  2:558407/1‑229 (MQ=255)
                                                                                  gcTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGGTGATGGTGATGGCGCTGTCCTTTATCTTTTGCCCCCCGCGGAGCGTACCGGTGGGCGGGCGGTGGGTCACCCGGATTCGTTTCGGGTTTTTCTTCTTCGTACttggcccagtcctct                                                                                                                                                                      >  2:431583/1‑237 (MQ=255)
                                                                                               gCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCCTCTTTGGGCCACTGCGTAGCGTACCGATGGGCGGGCGTGTGGTCACCGGTATCAGttt                                                                                                                                                                                                          >  1:5631/1‑203 (MQ=255)
                                                                                               gCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGttt                                                                                                                                                                                                          <  2:5631/203‑1 (MQ=255)
                                                                                                 aCTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGtcttct                                                                                                                                                                                            <  1:316399/215‑1 (MQ=255)
                                                                                                      ctATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGGTGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCACTGCGTAGCGTACCGATGGGCGGGCGTGTGGGCACCGGTATCAGGTTCGGGTTTTTCTTCTACCTACTGGACCAGGTCTTCGGCCCGCTGACGTTGGttttt                                                                                                                                                 >  1:327217/1‑249 (MQ=255)
                                                                                                        atcatcGGTTTGCACAACTATGGGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATgg                                                                                                                                               <  2:308245/251‑1 (MQ=255)
                                                                                                               gTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCACTGCGTAGGGTACCGCTGGGCGGGCGTGTGGTCACCGGGATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTTTTGCATcccc                                                                                                                                        >  1:409801/1‑250 (MQ=255)
                                                                                                                   gCACAACTATGTTAAGTATCGGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCccgccg                                                                                                                                    <  1:480084/250‑1 (MQ=255)
                                                                                                                       aaCTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGCAGCGTACCGATGGGCGTGCGTGTGGGCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATTTTCGGCCCGCTGACGTTGGTTTATTGCATCccgcccctcat                                                                                                                               >  1:166999/1‑245 (MQ=255)
                                                                                                                         cTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGTAGATGGTGATTGCGCCGTCGTTCATCTTTTGCCCACTGCGTAGCGGACCCATGGGCGGGCGGGGGGTCCCCGGGATCCGTTTCGGTTTTTTCTTCTCACTACTGGGCCCAATCTTCGGCCCGCTGTCGTTTGTTTATGGGATCccgcccaaccttc                                                                                                                             >  2:540928/1‑243 (MQ=255)
                                                                                                                             gtgAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGCCGTTTATCTTTGGCCCACTGCGTAGCGCCCCGATGGGCGTGCGTGTGGTCACCGGGATCAGTTTCGGGCTTGTCTTCTTCGTCCTGGATCAGATCTTCGGCCCGCTGACGCTGGGTTTTGGCATCCCGCCGGTCACCGgcgc                                                                                                                         >  1:344503/1‑251 (MQ=255)
                                                                                                                                aaGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCCTCTTTGGCCCCCTGCGTAGCGTACCGATGGGCGGGCGTGTGGGCCCCGGTATCAGTTTCGGTTTTTTTTTCTACGTACTGGACCAGATCTTCGGCCCGCTGGCGTTGGTTTTTGGCATCCCGCCGATCAACGGCGCAcc                                                                                                                      >  2:295558/1‑250 (MQ=255)
                                                                                                                                         aaGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGGTGGCGCTGTCGTTCCTCTTTTGCCCACTGCGTAGCGTACCGATGGGCGGGCGTGGGGTCACCGGGATCCGTTTCGGGTTTGTCTTCTTCGGAATGGACCAGATCTTCGGCCCCCTGTCGTTGGGTTTTTGCATCCCCCCGCTCATCCGCGCACTGTTGCCa                                                                                                             >  2:336210/1‑249 (MQ=255)
                                                                                                                                             cGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAgcgc                                                                                                        <  1:246144/250‑1 (MQ=255)
                                                                                                                                             cGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGGGGCGGTGATGGTGCTGATGGCGCTGTTGTTCCTCTTTCGCGCACTGCGGAGCGTATCGAAGGGCGCGCGGGTGGTCGGCGGGATCAGTTTCGGGTTTGTCTTCTTAGGACAGGGCCCGAACTTCGGCCCGCTTGAGGTGGTTTTTGGCCTACCGCCGGTCATCTgcgcgccatttccacgcgcc                                                                                                         >  2:10358/1‑235 (MQ=255)
                                                                                                                                               aGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTCGCCAAGCGCCAg                                                                                                     <  1:431583/251‑1 (MQ=255)
                                                                                                                                               aGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGAAGATGCTGGTGGAGCTGTCGTTCATCTTTTGCCCACCGCGTAGGGTAACGATGGGGGTGCGTGTGGTCAACGGCATCAGTTTCGGTTTTGTCTTTTACGTACTGGACCAGATCTTCGGCCCGCCTACCTTTGTTTTTGGGATCCCGCCCATCCACGGCGCACTGTTGCCAAACGcccc                                                                                                     >  2:437683/1‑249 (MQ=255)
                                                                                                                                                   gtcagtaggCCGGACGTTTTCAGCTCGACATGTGGAGCAAAATCTTCCGGCCGCTATCTGTGGCGGGGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGGGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTAGGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGcttc                                                                                                 <  2:344503/243‑1 (MQ=255)
                                                                                                                                                       ggATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGGTGGCGCCGTCGTTCATCTTTTGCCCACTGCGTAGCGGACCGATGGGCGTGCGTGGGGGCCCCGGGATCCGTTTCGGTTTTTTCTTTTACGGACTGGACCAGATCTTCCGCCCCCCGGCGTTGGTTTTTGGCATCCCCCCCATCATCGGCGCACCGTTGCCTAGCGCCCGcgtcttct                                                                                             >  2:121353/1‑251 (MQ=255)
                                                                                                                                                        gATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCCCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCCCTGTTGCCCAGCGCCAGcctcttctt                                                                                            >  1:53139/1‑251 (MQ=255)
                                                                                                                                                            ccGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATc                                                                                        <  1:169166/251‑1 (MQ=255)
                                                                                                                                                            ccGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATc                                                                                        <  1:314730/251‑1 (MQ=255)
                                                                                                                                                                  nttATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGACGATGCTGGTGGAGCTCTCCTTCCACTTTTGCCCCCCGCGTAGCGCAACGCTGGGCGTTCGTGGGGTCACCGGTATCCGTTTCGGTTTTGTTTTCTACGCACTGGGCCAGACCTTCGGCCCGCTGACGTTGGTTTATTGCCTCCCGCCCACCATCGCCGCACTCTTGACCAACGCCCGCGTCTTCTTactccaccct                                                                                  >  2:276279/2‑242 (MQ=255)
                                                                                                                                                                       cAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTCCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCt                                                                             <  1:128534/251‑1 (MQ=255)
                                                                                                                                                                        aGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTg                                                                            <  1:295558/251‑1 (MQ=255)
                                                                                                                                                                          cTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTTGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCCGTATCAGTTTCGGGTTTTTCTTCTACGTACTGGGCCAGATCTTCGGCCCGCCGACGTTGGGTTTTTGCATCCCCCCGATCATCGGCGCCCTGTTGCCAAGCGCCAGCTTCTTCTTTATCAGCCCGTGGCTggt                                                                          >  1:166848/1‑249 (MQ=255)
                                                                                                                                                                           tCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGAAGGCGCTGTCGTTCATCTTTTGCCCACCGCGTAGCGTACCGATGGGGGGGCGTGTGGTTACCGGGATCCGTTTCGGGTTTGTCTTCTACGGTCTGGGCCCGATTTTCGGCCCGCTGACGTTGGTTTTTGGCATCCCGCCGCTCATCCGCGCACTGTTGCCAAGGGCCAGCTTCTTCTTAATCAGCCCGTGGCCGttt                                                                         >  2:349701/1‑250 (MQ=255)
                                                                                                                                                                                                     ccGCTATCTGTGGCGGTGATGATGCTGATGGCGCTTTTGGTCTTCTTTTGGCCCCTGCGGAGCGGAACGATGGGCGTGCGTGTGGTCACCGGTATCAAGTTTGGGTTTGTCTTTTAACTACCGGACCAGATCTTTGGCCCGCTGACGTTTGTTTTTGGCATCCCCGCCCTCATAGGCGCACTGGTTGCCAGCGCCAGCTTCTTCTTAATCAGCCTGGGGCTGGTAAATAGAAAATccgaaaccaaagaa                                                 >  2:385534/1‑236 (MQ=255)
                                                                                                                                                                                                       gCTATCTGTGGTGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGa                                              >  1:95820/1‑250 (MQ=255)
                                                                                                                                                                                                                ggCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGGGTAGCGTATCGATGGGCGTGCTTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGCCCAGATCTTCGGCCCGCTGACGTTGGTTTATGTCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGa                                     <  1:10358/250‑1 (MQ=255)
                                                                                                                                                                                                                ggCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAggg                                      >  2:453844/1‑249 (MQ=255)
                                                                                                                                                                                                                 gCGGTGATGATGCTGATGGCGCTGTCGTCCATCTTTGCCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGTTATCAGTTTCGGTTTTGTCTTCTACGTACTGGCCCAGATCTCCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGACAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGc                                   <  1:483657/251‑1 (MQ=255)
                                                                                                                                                                                                                    gTGATGATGCTGATGGCGCTGTCGTTCTTCTTTGGCCCACTGCGTGGCGTACCGATGGGCGTGCGTGGGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCTAATGAACGAATAAAGGGAGCGaa                                <  2:218937/251‑1 (MQ=255)
                                                                                                                                                                                                                       atgatgCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGc                             >  1:73006/1‑251 (MQ=255)
                                                                                                                                                                                                                         gatgCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCt                            <  1:549047/250‑1 (MQ=255)
                                                                                                                                                                                                                           tgCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTCCCGATTGGCGTGCGTGTGGTCCCCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTcc                          <  1:187456/250‑1 (MQ=255)
                                                                                                                                                                                                                            gCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTccc                         <  1:267980/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gTCGTTCATCTTTGGCCCCCTGCGTAGCGTACCGTTGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGtt             <  1:116114/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ttCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCa         >  2:110354/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   ttGGCCCACTGCGTAGCGTACCGATGGGCGGGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGtgtg  <  2:53139/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                              
ATCCGAAACAGATTACCGGTTCGCAGACGGTGAGCGGCACCTGGAAAACCAACCTCACGCCGGACAAACTGGGCGTGGTGGCGCTGGACCCGGATGCACTCTCTATCAGCGGTTTGCACAACTATGTGAAGTATCTGAAGTCGAGCGGTCAGGATGCCGGACGTTATCAGCTCAACATGTGGAGCAAAATCTTCCAGCCGCTATCTGTGGCGGTGATGATGCTGATGGCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTG  >  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: