Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 328,824 T→C 14.4% S171G (AGC→GGC)  betI ← choline‑inducible betIBA‑betT divergent operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913328,8240TC14.4% 121.6 / 2.0 58S171G (AGC→GGC) betIcholine‑inducible betIBA‑betT divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (17/31);  new base C (5/4);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCC  >  NC_000913/328598‑329055
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       
catcccgcggcCCGGGGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCGGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGTTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGGGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCGGCGGTTTGCCGCTCAGAGCcgcgcg                                                                                                                                                                                                                             <  1:155777/230‑1 (MQ=255)
cATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCcgcgcg                                                                                                                                                                                                                             <  2:1221064/239‑1 (MQ=255)
cATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCcgcgcg                                                                                                                                                                                                                             <  2:1145060/239‑1 (MQ=255)
 aTCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAGCGGGGGGTTTTTTTCGCGGGTTTTCGCTCAGGACcccgcg                                                                                                                                                                                                                             >  2:155777/1‑238 (MQ=255)
   ccTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTGAGGGGGGTTGCGCGGGGTTTTTTCAGCGGGTTTGCGCCCCCAGCCccgcgcc                                                                                                                                                                                                                           >  2:224842/1‑237 (MQ=255)
           ccGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCGCGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAAt                                                                                                                                                                                                                   <  2:312055/238‑1 (MQ=255)
            cggcggCCTGCCCGGGCGCCAGCACGTTACCGTTGGCCGGGTTAAGGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCCGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATcc                                                                                                                                                                                                                 <  1:1008465/239‑1 (MQ=255)
                     gCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGGGGGGATTAGAGCGGGGGTTTTTCAGCCGGTTTGCCCCCACGGCCCCGCGCACGCCTAATCCctcactct                                                                                                                                                                                                         >  2:472152/1‑237 (MQ=255)
                              ccAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGGGTTTTTCAGCGGGTTGCCGCTCAGAGACGCGCGCAGGCATAATCCCTCAATAAGCGCGGCCAg                                                                                                                                                                                               >  2:272473/1‑239 (MQ=255)
                                     gTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGCGGGGAATAGAGCGGGGTTTTTTCAGCGGGTTTTCTCCCCCAGCCGCGCGCCGCCCCAAACCCTCAAACAGagcgcgcacggccta                                                                                                                                                                                         >  2:150063/1‑227 (MQ=255)
                                     gTTACCGTTGGCCGGGTGAATGGTCTCGAAGGTGGGACCGCTGGTGGCGGAGGTATACCCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAg                                                                                                                                                                                        <  2:78648/239‑1 (MQ=255)
                                         ccGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATGGGGGGGTAGATGCGGGGTGATAAAGTGGGGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTGGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCg                                                                                                                                                                                    <  1:597641/239‑1 (MQ=255)
                                             tGGCCGGTTTAATGGTCTCGAAGGTGGGACCGCGGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTc                                                                                                                                                                                <  2:960712/239‑1 (MQ=255)
                                                    gTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGATTTGCGCGGGGTTTTTTCAGCCGGTTTGCCCTCCCAGCCCCGCCCAGCCCTAATTCACTAAACAACGCGGCCCGGGCCTGGCCCGcctcctcttc                                                                                                                                                                          >  2:552302/1‑234 (MQ=255)
                                                              tcGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGATTAGTGCGGGGGTTTTTCAGCGGGTTTCCGCCCCGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCCGGCCGGAGCCCCCTTCCTGTGCCCGTTcgcgcg                                                                                                                                                               >  2:425727/1‑239 (MQ=255)
                                                              tcGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGAATTAGCGGGGGGGTTTTTCTCCTGGTTTTGGCTCCGAGCCGCGCGCCGCCGAAATCCACCACACCTCGCGGCCCAGGCGCCGCCCCCCCCCTGTGCCTCGTcgcccc                                                                                                                                                               >  1:739333/1‑236 (MQ=255)
                                                                          ccGCGGGGGGCGGGGGTATAACCACATTGTATATAAGGCTGCTGTTCTGCCATTCGGGAACTCGCGTCTCCTCGGTTAATGGGGGGGTAGATGCTGGGTGATAAAGGGGCGGGTCAGGGAATTAGCGGGGGTTTTATCCGGGGGTTGGCCGCTCAGGGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGcga                                                                                                                                                    <  1:332526/238‑1 (MQ=255)
                                                                                  gggggAGGTATACCCCCCATGTATAAAAGGCGGTTGTCCGGCCATTCGGGACAGCACGTCTCCGCGTTTAATCGGTGGGGAGATGCTGGGGGATAAAGTGGGGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTc                                                                                                                                            <  1:212269/235‑1 (MQ=255)
                                                                                         gTATAACCCCCTGTTATATAAAGCGGCTGTTCTGCCATGCGGGACATCCCGTCTCCTCGGTTAAGCGGTGGGTAGATGCTGAGGGATAAGGTGGGGGGTCAGGGAATTAGCGCGGGTTTTACCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCAcc                                                                                                                                     <  1:224842/238‑1 (MQ=255)
                                                                                          tataACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTc                                                                                                                                                        <  2:405445/218‑1 (MQ=255)
                                                                                                           taAAGTTGCTGTTCGGCCTTTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGGGGCGGGTCAGGGAATTAGCGCGGGTTTTACCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAgcggc                                                                                                                  <  2:980695/239‑1 (MQ=255)
                                                                                                                   ctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCGGGGAATTAGCGCGGGGGTTTTTTACCGGGTTTTCGCCCAACGCCCCGCGCCGGCCCAATCCATCACTCACCGCCGCCCGGACGCGGCCCGCCCCCCGCGCCTGTTTCCGCCGGGCGTCTCGGCCAAACTCCCCCACCCGATTCgcacgccaccggcgcggccg                                                                                                           >  1:213569/1‑219 (MQ=255)
                                                                                                                        ctcccATTCGGGACATCCCGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACctg                                                                                                     <  2:546554/236‑1 (MQ=255)
                                                                                                                          gcgtttggggACAGCCCGTCTCCGCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGGGGGGGGTCAGGGAATTGGCGCGGGTTTTACCCAGCGGTTGGCCGCTCAGACCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACctgct                                                                                                   <  1:1105035/232‑1 (MQ=255)
                                                                                                                          gCCATTCGGGACATCACGCCTCCTCGGTTAATCGGTGGGTAGAGCCTAAGTGATAAAGTGGCGGGTAAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCGCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACctgct                                                                                                   <  2:183690/239‑1 (MQ=255)
                                                                                                                                    acttcgcGTCTCCGCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTTTCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGAt                                                                                          <  2:1146274/232‑1 (MQ=255)
                                                                                                                                     cATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCGGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATa                                                                                         <  2:932505/238‑1 (MQ=255)
                                                                                                                                         aCGTCTCCGCGGTTAATCGGGGGGTAGATGCTGAGTGATAAAGTGGGGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCa                                                                                    <  2:1096666/239‑1 (MQ=255)
                                                                                                                                           gTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGGGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCaa                                                                                                                       <  2:814405/202‑1 (MQ=255)
                                                                                                                                             tgcccgggTTAAGGGGTGGGTAGATGCTGAGGGATAAAGGGGCGGGGCGGGGAATTACGGCGGGTTTTATCCAGGGGTTTGCCGCCCCGGCCCGCGCCCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCgg                                                                                <  1:640021/233‑1 (MQ=255)
                                                                                                                                             ctccttGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCgg                                                                                <  2:341121/239‑1 (MQ=255)
                                                                                                                                             ctcctgGGTTAAGGGGGGGGTAGATGCTGAGGGATAAAGGGGCGGGTCAGGGAATTAGCGGGGGTTTTACCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCgg                                                                                <  1:425727/239‑1 (MQ=255)
                                                                                                                                             cgcctgGGTTAATGGGGGGGTATGTGCGGGGTGATAAGGGGGGGGGTCAGGAAATTAGCGGGGGTTTTATCCGGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCgg                                                                                <  1:225057/237‑1 (MQ=255)
                                                                                                                                                              gggTAGATGCTGAGGGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATActgctg                                                                <  2:129484/238‑1 (MQ=255)
                                                                                                                                                                             gATAAAGTGGCGGGTCGGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACgccagc                                                 <  2:1218811/238‑1 (MQ=255)
                                                                                                                                                                               tAAAGGGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCCCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGccagccag                                              <  2:413659/239‑1 (MQ=255)
                                                                                                                                                                                       cGGGTCAGGGAATTAGCGCGGGTTTTATCCAGGGGTTTGCCGCTCAGAGCCGCCCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCa                                       <  2:811608/238‑1 (MQ=255)
                                                                                                                                                                                        gggTCAGGGATTTAGCGCGGGTTTTATCCGGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCAt                                      <  2:1127103/238‑1 (MQ=255)
                                                                                                                                                                                              gggAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATcgccgc                                >  1:157480/1‑238 (MQ=255)
                                                                                                                                                                                                gAATTAGGGCGGTTTTTATCCGGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCAc                              <  2:263393/238‑1 (MQ=255)
                                                                                                                                                                                                gAATTAGCGGGGGTTTTAGCCAGCGTTTTCCCGCTCAGACCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATcgccgc                                <  2:909351/236‑1 (MQ=255)
                                                                                                                                                                                                  aTTAGCGCGGGTTTTACCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTg                            <  2:417557/238‑1 (MQ=255)
                                                                                                                                                                                                   ttAGCGCGGTTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCt                          <  2:50234/239‑1 (MQ=255)
                                                                                                                                                                                                   ttAGCGCGGGTTTTATCCAGCGGTTTGCCGCCCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCt                          <  2:1164548/239‑1 (MQ=255)
                                                                                                                                                                                                         cgGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCCGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCt                     >  1:214027/1‑238 (MQ=255)
                                                                                                                                                                                                         cgGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCt                     >  2:348874/1‑238 (MQ=255)
                                                                                                                                                                                                           ggTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATCGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCg                  >  1:302817/1‑239 (MQ=255)
                                                                                                                                                                                                            gTTTTACCCAGCGTTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCg                  <  2:914472/238‑1 (MQ=255)
                                                                                                                                                                                                                      gCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGCAg       >  1:408192/1‑239 (MQ=255)
                                                                                                                                                                                                                      gCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAg       >  2:5716/1‑239 (MQ=255)
                                                                                                                                                                                                                      gCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAg       >  2:1162737/1‑239 (MQ=255)
                                                                                                                                                                                                                      gCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGCTTCATCGAAg       >  1:31406/1‑239 (MQ=255)
                                                                                                                                                                                                                      gCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATAGCCGCACTGCTCACCTGCGTTTCATCGAAg       >  2:708671/1‑239 (MQ=255)
                                                                                                                                                                                                                        cGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGt      >  2:986264/1‑239 (MQ=255)
                                                                                                                                                                                                                          gTTTGCCGCCCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGtt     <  2:783383/238‑1 (MQ=255)
                                                                                                                                                                                                                             tGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTcc  >  2:1189194/1‑238 (MQ=255)
                                                                                                                                                                                                                             tGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTcc  >  2:253375/1‑238 (MQ=255)
                                                                                                                                                                                                                             tGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCCCACCTCTCACCTGCGTTTCCTCGAAGTTTcc  >  1:886092/1‑238 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       
CATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCC  >  NC_000913/328598‑329055

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