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

CATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCC  >  NC_000913/328598‑329055
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       
catcccgcggcCCGGGGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCGGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGTTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGGGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCGG‑CGGTTTGCCGCTCAGAGCcgcgcg                                                                                                                                                                                                                             <  1:155777/230‑1 (MQ=255)
cATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCcgcgcg                                                                                                                                                                                                                             <  2:1221064/239‑1 (MQ=255)
cATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCcgcgcg                                                                                                                                                                                                                             <  2:1145060/239‑1 (MQ=255)
 aTCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAGCGGGGGGTTTTTTTCG‑CGGGTTTTCGCTCAGGACcccgcg                                                                                                                                                                                                                             >  2:155777/1‑238 (MQ=255)
   ccTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTGAGGGGGGTTGCGCGGGGTTTTTTCAG‑CGGGTTTGCGCCCCCAGCCccgcgcc                                                                                                                                                                                                                           >  2:224842/1‑237 (MQ=255)
           ccGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCGCGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAAt                                                                                                                                                                                                                   <  2:312055/238‑1 (MQ=255)
            cggcggCCTGCCCGGGCGCCAGCACGTTACCGTTGGCCGGGTTAAGGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCCG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATcc                                                                                                                                                                                                                 <  1:1008465/239‑1 (MQ=255)
                     gCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGGGGGGATTAGAGCGGGGGTTTTTCAGCCGGTTTGCC‑CCCACGGCCCCGCGCACGCCTAATCCctcactct                                                                                                                                                                                                         >  2:472152/1‑237 (MQ=255)
                              ccAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGGGTTTTTCAG‑CGGGTTGCCGCTCAGAGACGCGCGCAGGCATAATCCCTCAATAAGCGCGGCCAg                                                                                                                                                                                               >  2:272473/1‑239 (MQ=255)
                                     gTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGCGGGGAATAGAGCGGGGTTTTTTCAG‑CGGGTTTTCTCCCCCAGCCGCGCGCCGCCCCAAACCCTCAAACAGagcgcgcacggccta                                                                                                                                                                                         >  2:150063/1‑227 (MQ=255)
                                     gTTACCGTTGGCCGGGTGAATGGTCTCGAAGGTGGGACCGCTGGTGGCGGAGGTATACCCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAg                                                                                                                                                                                        <  2:78648/239‑1 (MQ=255)
                                         ccGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATGGGGGGGTAGATGCGGGGTGATAAAGTGGGGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTGGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCg                                                                                                                                                                                    <  1:597641/239‑1 (MQ=255)
                                             tGGCCGGTTTAATGGTCTCGAAGGTGGGACCGCGGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTc                                                                                                                                                                                <  2:960712/239‑1 (MQ=255)
                                                    gTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGATTTGCGCGGGGTTTTTTCAGCCGGTTTGCC‑CTCCCAGCCCCGCCCAGCCCTAATTCACTAAACAACGCGGCCCGGGCCTGGCCCGcctcctcttc                                                                                                                                                                          >  2:552302/1‑234 (MQ=255)
                                                              tcGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGATTAGTGCGGGGGTTTTTCAG‑CGGGTTTCCGCCCCGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCCGGCCGGAGCCCCCTTCCTGTGCCCGTTcgcgcg                                                                                                                                                               >  2:425727/1‑239 (MQ=255)
                                                              tcGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGAATTAGCGGGGGGGTTTTTCTC‑CTGGTTTTGGCTCCGAGCCGCGCGCCGCCGAAATCCACCACACCTCGCGGCCCAGGCGCCGCCCCCCCCCTGTGCCTCGTcgcccc                                                                                                                                                               >  1:739333/1‑236 (MQ=255)
                                                                          ccGCGGGGGGCGGGGGTATAACCACATTGTATATAAGGCTGCTGTTCTGCCATTCGGGAACTCGCGTCTCCTCGGTTAATGGGGGGGTAGATGCTGGGTGATAAAGGGGCGGGTCAGGGAATTAGCGGGGGTTTTATCCGG‑GGGTTGGCCGCTCAGGGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGcga                                                                                                                                                    <  1:332526/238‑1 (MQ=255)
                                                                                  gggggAGGTATACCCCCCATGTATAAAAGGCGGTTGTCCGGCCATTCGGGACAGCACGTCTCCGCGTTTAATCGGTGGGGAGATGCTGGGGGATAAAGTGGGGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTc                                                                                                                                            <  1:212269/235‑1 (MQ=255)
                                                                                         gTATAACCCCCTGTTATATAAAGCGGCTGTTCTGCCATGCGGGACATCCCGTCTCCTCGGTTAAGCGGTGGGTAGATGCTGAGGGATAAGGTGGGGGGTCAGGGAATTAGCGCGGGTTTTACCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCAcc                                                                                                                                     <  1:224842/238‑1 (MQ=255)
                                                                                          tataACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTc                                                                                                                                                        <  2:405445/218‑1 (MQ=255)
                                                                                                           taAAGTTGCTGTTCGGCCTTTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGGGGCGGGTCAGGGAATTAGCGCGGGTTTTACCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAgcggc                                                                                                                  <  2:980695/239‑1 (MQ=255)
                                                                                                                   ctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCGGGGAATTAGCGCGGGGGTTTTTTAC‑CGGGTTTTCGCCCAACGCCCCGCGCCGGCCCAATCCATCACTCACCGCCGCCCGGACGCGGCCCGCCCCCCGCGCCTGTTTCCGCCGGGCGTCTCGGCCAAACTCCCCCACCCGATTCgcacgccaccggcgcggccg                                                                                                           >  1:213569/1‑219 (MQ=255)
                                                                                                                        ctcccATTCGGGACATCCCGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACctg                                                                                                     <  2:546554/236‑1 (MQ=255)
                                                                                                                          gcgtttggggACAGCCCGTCTCCGCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGGGGGGGGTCAGGGAATTGGCGCGGGTTTTACCCAG‑CGGTTGGCCGCTCAGACCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACctgct                                                                                                   <  1:1105035/232‑1 (MQ=255)
                                                                                                                          gCCATTCGGGACATCACGCCTCCTCGGTTAATCGGTGGGTAGAGCCTAAGTGATAAAGTGGCGGGTAAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCGCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACctgct                                                                                                   <  2:183690/239‑1 (MQ=255)
                                                                                                                                    acttcgcGTCTCCGCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTTTCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGAt                                                                                          <  2:1146274/232‑1 (MQ=255)
                                                                                                                                     cATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCGGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATa                                                                                         <  2:932505/238‑1 (MQ=255)
                                                                                                                                         aCGTCTCCGCGGTTAATCGGGGGGTAGATGCTGAGTGATAAAGTGGGGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCa                                                                                    <  2:1096666/239‑1 (MQ=255)
                                                                                                                                           gTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGGGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCaa                                                                                                                       <  2:814405/202‑1 (MQ=255)
                                                                                                                                             tgcccgggTTAAGGGGTGGGTAGATGCTGAGGGATAAAGGGGCGGGGCGGGGAATTACGGCGGGTTTTATCCAG‑GGGTTTGCCGCCCCGGCCCGCGCCCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCgg                                                                                <  1:640021/233‑1 (MQ=255)
                                                                                                                                             ctccttGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCgg                                                                                <  2:341121/239‑1 (MQ=255)
                                                                                                                                             ctcctgGGTTAAGGGGGGGGTAGATGCTGAGGGATAAAGGGGCGGGTCAGGGAATTAGCGGGGGTTTTACCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCgg                                                                                <  1:425727/239‑1 (MQ=255)
                                                                                                                                             cgcctgGGTTAATGGGGGGGTATGTGCGGGGTGATAAGGGGGGGGGTCAGGAAATTAGCGGGGGTTTTATCCGG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCgg                                                                                <  1:225057/237‑1 (MQ=255)
                                                                                                                                                              gggTAGATGCTGAGGGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATActgctg                                                                <  2:129484/238‑1 (MQ=255)
                                                                                                                                                                             gATAAAGTGGCGGGTCGGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACgccagc                                                 <  2:1218811/238‑1 (MQ=255)
                                                                                                                                                                               tAAAGGGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCCCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGccagccag                                              <  2:413659/239‑1 (MQ=255)
                                                                                                                                                                                       cGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑GGGTTTGCCGCTCAGAGCCGCCCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCa                                       <  2:811608/238‑1 (MQ=255)
                                                                                                                                                                                        gggTCAGGGATTTAGCGCGGGTTTTATCCGG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCAt                                      <  2:1127103/238‑1 (MQ=255)
                                                                                                                                                                                              gggAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATcgccgc                                >  1:157480/1‑238 (MQ=255)
                                                                                                                                                                                                gAATTAGGGCGGTTTTTATCCGG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCAc                              <  2:263393/238‑1 (MQ=255)
                                                                                                                                                                                                gAATTAGCGGGGGTTTTAGCCAG‑CGTTTTCCCGCTCAGACCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATcgccgc                                <  2:909351/236‑1 (MQ=255)
                                                                                                                                                                                                  aTTAGCGCGGGTTTTACCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTg                            <  2:417557/238‑1 (MQ=255)
                                                                                                                                                                                                   ttAGCGCGGTTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCt                          <  2:50234/239‑1 (MQ=255)
                                                                                                                                                                                                   ttAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCCCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCt                          <  2:1164548/239‑1 (MQ=255)
                                                                                                                                                                                                         cgGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCCGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCt                     >  1:214027/1‑238 (MQ=255)
                                                                                                                                                                                                         cgGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCt                     >  2:348874/1‑238 (MQ=255)
                                                                                                                                                                                                           ggTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATCGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCg                  >  1:302817/1‑239 (MQ=255)
                                                                                                                                                                                                            gTTTTACCCAG‑CGTTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCg                  <  2:914472/238‑1 (MQ=255)
                                                                                                                                                                                                                      g‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGCAg       >  1:408192/1‑239 (MQ=255)
                                                                                                                                                                                                                      g‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAg       >  2:5716/1‑239 (MQ=255)
                                                                                                                                                                                                                      g‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAg       >  2:1162737/1‑239 (MQ=255)
                                                                                                                                                                                                                      g‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGCTTCATCGAAg       >  1:31406/1‑239 (MQ=255)
                                                                                                                                                                                                                      g‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATAGCCGCACTGCTCACCTGCGTTTCATCGAAg       >  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)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       
CATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAG‑CGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCC  >  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: ‑