Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,965,5740CA25.4% 89.7 / 40.4 55K19N (AAG→AATrhlBATP‑dependent RNA helicase RhlB
Reads supporting (aligned to +/- strand):  ref base C (17/24);  new base A (3/11);  total (20/35)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACG  >  NC_000913/3965325‑3965824
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          
ggCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCt                                                                                                                                                                                                                                                           >  2:53240/1‑251 (MQ=255)
      gATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGAt                                                                                                                                                                                                                                                     >  1:200642/1‑251 (MQ=255)
              gggCAAGTTCCCGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCg                                                                                                                                                                                                                                             <  1:192015/250‑1 (MQ=255)
              cGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCg                                                                                                                                                                                                                                             >  1:137673/1‑251 (MQ=255)
                   aGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCAGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGACTGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCTAAGTc                                                                                                                                                                                                                                        >  2:373550/1‑251 (MQ=255)
                     ttCACGCGGCGGTGCCATAATAAAGGCAAGCGGCGTATTCACCTTGCGATCGTCAATCGCAGGATGAGAGGGAAGATAATGAAACGTGGACGTAAGAAACGCCATCGTTTTCCCGGTCCCGGTTTGCGCCGCCCGGGTTACGCCACGCCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCg                                                                                                                                                                                                                                       <  2:12477/250‑1 (MQ=255)
                      gcACGCGTCGGCGCAATAATTAAGGACCGCGGCTGATTCATCTTGCGAGCGGCAATCGCAGGATGAGAGAGACGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGCCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCgg                                                                                                                                                                                                                                      >  2:153675/2‑250 (MQ=255)
                       cACGCGTCGGTGCCAAAATTAAGGCACGCGGTTGATCCCGCTTGCGAGCGGCAATAGCAGGATGAGAGAGAAGATAATGAACCGGTGACGTAAGAACCGCCATCGTTTCCCCGGCACCGGTTTGCGCCGGCCCGGCTACGTCCCGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACATTCGGATGCAGGGCGAAGTCGGaa                                                                                                                                                                                                                                    <  2:47207/251‑1 (MQ=255)
                           ggTCGGTGCCATAATTAAGGCACGCGGCTGATTCCCCTTGCGATCGGCAACCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAACCGCCATCTTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACt                                                                                                                                                                                                                                <  2:200642/250‑1 (MQ=255)
                              cGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACATTCGGATGCAgg                                                                                                                                                                                                                                                >  1:92572/1‑232 (MQ=255)
                              cGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACATTCGGATGCAGGGCGAAGTCGGAAAACttct                                                                                                                                                                                                                             <  1:79726/251‑1 (MQ=255)
                              cGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAAGCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACATTCGGATGCAgg                                                                                                                                                                                                                                                <  2:92572/232‑1 (MQ=255)
                                                  ggggcggATTCCCCTTGCGACCGGCAACGGCAGGATGGGAAGGAAGATAAGGAAACGTTGACGTAAGAAACGCCACCGTTTTCCCGGTACCGGTTTGCCCCTGCCCGACTACGTCCCCACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACATTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCGGTTAAATGTGTTTTg                                                                                                                                                                                                         <  1:305939/245‑1 (MQ=255)
                                                        gATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAATCGTTGAGGTAAGAACCGCCATCTTTTCCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACATTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCAt                                                                                                                                                                                                    <  1:110849/250‑1 (MQ=255)
                                                           tCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAgt                                                                                                                                                                                                 <  1:267611/250‑1 (MQ=255)
                                                                   cGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGa                                                                                                                                                                                         >  1:189928/1‑250 (MQ=255)
                                                                    gATCGGCAATCGCAGGCTGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCCAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAACACTTCTGTTCTGTTCAATGTGTTTTGCTCATAGTGTGGTAgcat                                                                                                                                                                                       >  1:341815/1‑248 (MQ=255)
                                                                       cGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTCCCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACATTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATc                                                                                                                                                                                    <  1:282260/251‑1 (MQ=255)
                                                                         gcattCGCAGGATGAGAGAGAAGATAATTAACCGTTGACGTAAGCACCGCCATCGTTTCCCCGGTACCGGTTTACCCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATTTGTTTTGCTCATAGTGTGGTAGAATATCAg                                                                                                                                                                                  <  1:34279/247‑1 (MQ=255)
                                                                                     tGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGAGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACCACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCt                                                                                                                                                                      >  2:225798/1‑251 (MQ=255)
                                                                                      gagagagaAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACATTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCtt                                                                                                                                                                     <  1:102407/251‑1 (MQ=255)
                                                                                                               cGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGCACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGCTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGGGAAATa                                                                                                                                             >  1:157351/1‑250 (MQ=255)
                                                                                                                gtaaaaaccGCCATGGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTGCGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACATTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATaaa                                                                                                                                           <  1:282762/243‑1 (MQ=255)
                                                                                                                gTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTTTGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACt                                                                                                                                                                            <  1:228710/218‑1 (MQ=255)
                                                                                                                gTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTTTGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACt                                                                                                                                                                            >  2:228710/1‑218 (MQ=255)
                                                                                                                  aaGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGt                                                                                                                                         <  1:173557/251‑1 (MQ=255)
                                                                                                                   aGAAACCCCTTCGTTTCCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTc                                                                                                                                        <  2:382161/251‑1 (MQ=255)
                                                                                                                    gAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTACCCGGCTACGTCACGACTCGCCCGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCCAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGACGGTAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTCGAATATCCGCTTACTATTGCTTTAAGAAAGCGTATCCGGTGAAATAAAGTc                                                                                                                                        >  2:103299/1‑250 (MQ=255)
                                                                                                                            ctctttgtCCCGGTACCGGTTTGCTCCTGCCCGCCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAtgt                                                                                                                                                                                                               <  2:325983/164‑1 (MQ=255)
                                                                                                                            aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAtgt                                                                                                                                                                                                               >  1:325983/1‑171 (MQ=255)
                                                                                                                            aTCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCttt                                                                                                                                 >  1:338364/1‑249 (MQ=255)
                                                                                                                                ttttCCCGTTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTg                                                                                                                            <  1:131799/250‑1 (MQ=255)
                                                                                                                                 tttCCCGGGACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGt                                                                                                                          >  1:400615/1‑251 (MQ=255)
                                                                                                                                  ttCCCGGTACCTGTTGGCGCCTGCCCGGCTACTTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGt                                                                                                                          <  2:133175/250‑1 (MQ=255)
                                                                                                                                     ccGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACATTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAAt                                                                                                                      >  1:300222/1‑251 (MQ=255)
                                                                                                                                      gtgTACCGGTGTGCGCCGGCCCGGGTACGGCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATg                                                                                                                     <  2:341815/249‑1 (MQ=255)
                                                                                                                                      cGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATg                                                                                                                     <  1:331147/251‑1 (MQ=255)
                                                                                                                                            cGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTacac                                                                                                               <  1:141828/251‑1 (MQ=255)
                                                                                                                                                            ggCTAAGTCACGACCCGCCAGCGTCAGCGGTCGGGCCAGTGCCTGAATGGGCGTACAGTTATGCAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAac                                                                                               >  2:113425/1‑251 (MQ=255)
                                                                                                                                                                         aCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTa                                                                                   <  1:373550/251‑1 (MQ=255)
                                                                                                                                                                               cAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCtg                                                                             <  1:10391/251‑1 (MQ=255)
                                                                                                                                                                                      aGCGGAAGTGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACATCCTTCGGATTCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGtt                                                                      <  1:157772/251‑1 (MQ=255)
                                                                                                                                                                                          gAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTtatat                                                                  <  1:415044/251‑1 (MQ=255)
                                                                                                                                                                                          gAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACATTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTtatat                                                                  <  2:331153/251‑1 (MQ=255)
                                                                                                                                                                                              ggCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCACCCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATg                                                               >  2:413257/1‑250 (MQ=255)
                                                                                                                                                                                                ccAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGc                                                            <  2:26124/251‑1 (MQ=255)
                                                                                                                                                                                                     gCCTGAGTGGGGTTACAGTTATGAAACCCTGTTTTTTCAAGGGCTTCTTCAACCTTCGGATGCGGGCCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATTGTGTGGTAGAAGATCAGCTTACAAGTGCTTTACGATAGCGTATCCGGTAAAATAAAGTCAGCCTTGAGTTGGGTAAGGTTACGCCAACACCGAAATCAACACGCTTGGCTTATTCCTGTGGAGTTATATATGAGCGATaa                                                       <  1:373582/251‑1 (MQ=255)
                                                                                                                                                                                                     gCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACATTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATaa                                                       <  1:405860/251‑1 (MQ=255)
                                                                                                                                                                                                        ggAATGGGGGTACAGTTATGAACCCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTAGCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAat                                                    <  2:137673/250‑1 (MQ=255)
                                                                                                                                                                                                          aaTGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAatt                                                   <  1:28943/250‑1 (MQ=255)
                                                                                                                                                                                                               gCGTACAGTTATGAAACCCTTGTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTGCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCACCCTTTAGTTGGTTAATGTTACACCAACAAGGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGTTAAAATTATTccc                                             <  2:147985/251‑3 (MQ=255)
                                                                                                                                                                                                                 gtaaaGTTAAGAAACCCCTTTTTTTCAAGGGCTTTTTAAACATTCGGATGCAGGGGGAAGTCGGAAAACTTCTTGTCTGGTAAATGGGTTTTGCTCAAAGTGTGGTAGAATATCAGCTTACTTTTGCTTTACGAAAGCGTAACCGGTTAAAAAAAATCAACCTTTAGTTGGGTAATTGTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGTAGTTATATATGAGCGATAAAATTATTCACct                                           <  1:311068/247‑1 (MQ=255)
                                                                                                                                                                                                                         aTGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTCCTATTGCTTTACGAACGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCtgactgac                                    <  2:291042/250‑1 (MQ=255)
                                                                                                                                                                                                                                              ggCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTAc              >  1:63351/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  tCTACAACATTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCa           >  1:216317/1‑250 (MQ=255)
                                                                                                                                                                                                                                                          cTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACg  >  2:380297/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          cTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCAATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACg  >  2:179535/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          
GGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAATAAAGTCAACCTTTAGTTGGTTAATGTTACACCAACAACGAAACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACG  >  NC_000913/3965325‑3965824

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