Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,062,576 G→A 9.8% L94F (CTT→TTT)  cbpM ← modulator of CbpA co‑chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,062,5760GA 9.8% 190.9 / 4.7 74L94F (CTT→TTT) cbpMmodulator of CbpA co‑chaperone
Reads supporting (aligned to +/- strand):  ref base G (29/36);  new base A (6/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

GGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTCAG  >  NC_000913/1062344‑1062797
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                               
ggCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGAAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACCCATCACG‑GATAGGCTCAAAACCGGGAAAGCCGcg                                                                                                                                                                                                                            >  1:501705/1‑237 (MQ=255)
                 tGCTGGTGTTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGAACGGTTAAAATGCGTCATCGCATCCGGCAAAAACACATCACG‑GGAGAGCTAAAAACCGGGAAAACCGCCGGGGAAGCAGg‑gggggtt                                                                                                                                                                                                         >  2:443354/1‑235 (MQ=255)
                 ggCTGGTGTTAGTTGAGTGCGACCCCCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGtttt                                                                                                                                                                                                         <  1:443354/238‑1 (MQ=255)
                  gCTGGTGTTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCCCAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGGCATCGCCTCCGGCAAACACCCATCACG‑GATGAGCTACAAAACGGGAAACACGCCGGCGCCGCAGG‑CGGgttt                                                                                                                                                                                                         >  1:275847/1‑238 (MQ=255)
                      gtgtTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTTTAGACCGGTTAAGATGCGTCATCGCATCCGGCCAAAACACATCACC‑GATGCGCTAAAAACCGGGAAAACCCCTGGCGCAGCAGG‑CGGTTTTccc                                                                                                                                                                                                      >  1:509379/1‑236 (MQ=255)
                       tgtTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGTTGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGAAGCGTCATCGCATCCGGCAAAAACCCATCACG‑GATGGGCTACAAACCGGGAAAGCCAGCGGCCCACCAGG‑CGGTTTTccccc                                                                                                                                                                                                    >  2:791036/1‑235 (MQ=255)
                       tgtTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGGCCAGCAGC‑CGGTTTTCCTg                                                                                                                                                                                                     >  1:530862/1‑237 (MQ=255)
                       tgtTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGGTAAGATGCGTCATCGCATCCGGCCAACACACATAACG‑GATGAGTCACAAACCTGAAAAACCGCCAGCCAAACAGG‑CGGGTTTccgga                                                                                                                                                                                                    >  2:317274/1‑235 (MQ=255)
                       tgtTAGCGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTCCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCCTTATGCCGCAGTTGCCGTATGCGGCGTGAACGTCTTATCCGGCCAACAGGCACTGGACTCTTTGTACACCGGTTAAAATGCGGCCTACCAACCGGCAAACACACATCAAC‑GATGAGCTACAACCCCGGAAAGACGCTGGCTCACACGG‑CGGtttttgtt                                                                                                                                                                                                     >  2:404269/1‑233 (MQ=255)
                        gtTAGTGGAGTGCGACCCCCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGtttt                                                                                                                                                                                                         <  1:822259/232‑1 (MQ=255)
                           tggggggTGCGACCCTCAAAACAGATAACGTCGTTTGGGCGCCTGAAAAGACGGGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCATG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAg                                                                                                                                                                                               <  2:316601/233‑1 (MQ=255)
                           aGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCCCG‑GATGAGCTCCAAACCGGGCAAACCGCTGGCGCAGCAGG‑CGGGTTTCCTGGTTCAg                                                                                                                                                                                               >  1:190812/1‑239 (MQ=255)
                           aGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAg                                                                                                                                                                                               <  1:11227/239‑1 (MQ=255)
                           aGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAAGCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACCCACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCCGCTTTAg                                                                                                                                                                                               >  1:821263/1‑239 (MQ=255)
                                 ggtcGACCACCAAAACAGATAACGTCGTTTGTGCGCCTTAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTgcg                                                                                                                                                                                          <  1:632117/235‑1 (MQ=255)
                                 gTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACAAACATCACG‑GCTGGGCTACAAACCCGGAAAACCGCCGGGCCAGCAAG‑CGGGTTTTCTGCCTcgcgggcac                                                                                                                                                                                         >  2:827665/1‑231 (MQ=255)
                                 gTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACCGGAACTTGAATCTTTGTAGACCGGTTAAGCTGCGCCACCCAATCCGGCCACCACCCA‑CTCGCGAGGACCACCAAACCGGGAAAACCCCGGGACCAACAGGAGGGTTTTCCTGCTCCCGGAgcg                                                                                                                                                                                          >  2:686864/1‑239 (MQ=255)
                                       ccaccaAAAAAGATAACGTTGTTTGTGCCCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATatc                                                                                                                                                                                   <  1:234502/239‑1 (MQ=255)
                                         accaAAACAGATAACGTCGTTTGGGCGCCTGAAAAGGCGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATatcat                                                                                                                                                                                 <  1:544048/239‑1 (MQ=255)
                                                    tAACGTCGTTTGTGCTCCTTAAAAGACGCGTTTTGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAACCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGt                                                                                                                                                                       <  2:190812/238‑1 (MQ=255)
                                                            tttGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCCGCTTCAGGGGCGCAATATCATCCCTTaaagccacgccca                                                                                                                                                               >  1:584791/1‑227 (MQ=255)
                                                              ttgcccccTAAAAAG‑CGCTTTTAGCGTCCCATCAGGCATTAGTGCTCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCGTTGTAGCCCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACcgc                                                                                                                                                           <  2:710006/233‑1 (MQ=255)
                                                               ttgCGCTTGAAAAGCGGCGTTTAGCGTCGCTTCAGGCATTTGGGCGCAGTTTCCGGATGGGGCGGGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGGAGGCCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACcgc                                                                                                                                                           <  2:276045/238‑1 (MQ=255)
                                                               tgtcggcTGAAAAGGCGCGTTTAGCGTCGCATCAGGCATGATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGGCAGCGCCACcgc                                                                                                                                                           <  1:352317/233‑1 (MQ=255)
                                                                   gcCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATTCCG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGGGCGCAATATAATCCATTAACGTCAGCGCCCCCCCGAt                                                                                                                                                        >  1:833270/1‑238 (MQ=255)
                                                                          aaGGCGCGTTTAGCGTCCCATCAGGCATTAGGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGcc                                                                                                                                                <  1:791036/239‑1 (MQ=255)
                                                                                            cactaGGGATTATGGGCCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTGTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCCACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCAt                                                                                                                              <  1:761354/235‑1 (MQ=255)
                                                                                                   cATTAGGGGGCGGTTGCCGGATGGGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACg                                                                                                                          <  1:335139/236‑1 (MQ=255)
                                                                                                   cATTAGGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTAACCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACg                                                                                                                          <  1:7149/236‑1 (MQ=255)
                                                                                                      tagggCGCAGTTGCTGGATGCGGCGTTAACGTCTTATCCTGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGATCTACAAACCGGGAAAGCCGCTGGCGCTGCAGG‑CGGTTTTCCTGCTTCAg                                                                                                                                                                                               <  1:886645/161‑1 (MQ=255)
                                                                                                           cgcAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTgcg                                                                                                                                                                                          >  2:83402/1‑164 (MQ=255)
                                                                                                           cgcAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTgcg                                                                                                                                                                                          <  1:83402/164‑1 (MQ=255)
                                                                                                              aGTTGCCGGATGCGGCGTGAACGGCTTATCCGGCCCACAGGAACTGGAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACcgcg                                                                                                             <  1:815822/238‑1 (MQ=255)
                                                                                                                ttGCCGGAGGGGGGGTTAACGTTTTAACCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACcgcgcg                                                                                                           <  1:343135/238‑1 (MQ=255)
                                                                                                                        tGCGGCGTGACCTGTTTTTCCGGCCCACGGGAACTGTAATCTTTGGAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGcacca                                                                                                   <  2:899558/238‑1 (MQ=255)
                                                                                                                             cgTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGGCAGCGCCACCCCGATCCCCGGCCAGTCCAGAGCCAGGTCCTTACCCCAGCGTACCGCCCGGTGCACCACaaag                                                                                              >  1:229652/1‑236 (MQ=255)
                                                                                                                                  aCGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCa                                                                                        <  1:317274/239‑1 (MQ=255)
                                                                                                                                   cGTCTTTCCCGGCCCACAGGAACTGTAATTTTTGGAGACCGGTTAAGATGCGGCATCGCATCCGGAAAACACCCAGCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCa                                                                                        <  2:229652/238‑1 (MQ=255)
                                                                                                                                     ttttttCCGGCCCACAGGAATGGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATgg                                                                                     <  2:81838/234‑1 (MQ=255)
                                                                                                                                            cGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCGAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCCCCACAATGGCGGCATGGACGT‑Caa                                                                             >  1:808048/1‑239 (MQ=255)
                                                                                                                                              gCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCCGAGCCCGTTCCTGACGCCAGCGTACCGCGCGTTGCACCACAATGGCGGCCTGGTCGG‑CAAAt                                                                           >  1:310705/1‑239 (MQ=255)
                                                                                                                                                     ggcaCTGTTATTTTTTTAGCCCGGTTAAGTTGCGTCATCGGTTCCGGCAAACCCAACTCCCG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAgg                                                                    <  2:420687/236‑1 (MQ=255)
                                                                                                                                                        aCTGTAATTGTTTTAGGCCGGTTAATATGGCGCTCCCCACCCGGCAAACACACATCACG‑GATGAGTCACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGTCATGGTCGT‑CAAATCCCCAGgttg                                                                 <  2:821263/239‑1 (MQ=255)
                                                                                                                                                         cTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCACCGCCACCGCGATCCCCGGCCAGTCCATAGCCAGTTCATTACGCCAGCGTACCGCGCGTTGCACCACAATGGCGGACTGGGCGT‑CAAATACCCCGgctgt                                                                >  1:436317/1‑239 (MQ=255)
                                                                                                                                                         atttAATTTTTGTAGACCGGTTAAGATGCGTCATCGCCTCCGGCAAACACACCTCACG‑GATGAGCTAAAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGgttg                                                                 <  2:310705/235‑1 (MQ=255)
                                                                                                                                                          tGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACcgc                                                                                                                                                           <  1:240733/148‑1 (MQ=255)
                                                                                                                                                            tAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAAAGCCAGTTCATGACGCAGGCGGACCGCGCGGTGCACCACAATGGCGGCATGGTCGC‑CAAATACCCAGGTGGttt                                                              >  1:847079/1‑238 (MQ=255)
                                                                                                                                                            tAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGGTTTCCTGCTTCAGGTGCGCAATATCATCCCTTAACGTCAGCGCCACCGCGATCCCCGGCCCGTCCACACCCCGTTCATCACCCGGGAGTCCCCCCGGTTGCCCCACAAAGGGCGCAGGGGGGC‑CCAAAACCCCGGTGGttt                                                              >  2:757653/1‑238 (MQ=255)
                                                                                                                                                              aTCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCCTGACCCAGGCGTACCGCCCGTTGCACCCCAATGCCGCCAAGGCCGT‑CAAAAACCCAGGTTGTTTccc                                                           >  1:420687/1‑238 (MQ=255)
                                                                                                                                                              aTCTTTGGAGACCGGTTAAAAAGCCTCAACGCATCCGGCAAACCCACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCGCCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCt                                                           <  1:397875/239‑1 (MQ=255)
                                                                                                                                                                 tttGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCCCGCGTTGCACCCCAATGCGGGCATGGTCGT‑CAAAAACCCAGGTTGTTTCCTgg                                                         >  1:671723/1‑238 (MQ=255)
                                                                                                                                                                        aacGGTTAAACTGGCTCATCGACTCCGGAAAACACGAGTCCAG‑GATGGGCTACAAACCGGGAAAGCCGCTTGCGGAGTAGG‑CGGTTTGCCTGCTTCAGGTGCCTAATATCATCCCTTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCAcg                                                 <  2:254561/237‑1 (MQ=255)
                                                                                                                                                                         cggtttAAGATGCGTCATCGCCTACGGCAAACACACATCAGG‑GATGAGCCACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑GGGTTTTCCTGATTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCAcg                                                 <  2:609465/234‑1 (MQ=255)
                                                                                                                                                                         ccGGTTAAGATGCGTCATCGCATCCGGCAAACACACAACACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTCCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCAcg                                                 <  1:130404/238‑1 (MQ=255)
                                                                                                                                                                         ccGGTTAAGAGGCGTCAGCGCATCCGGCAAACCCACATCCCG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCAcg                                                 <  1:396185/238‑1 (MQ=255)
                                                                                                                                                                          cGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATg                                                                                      >  1:509730/1‑201 (MQ=255)
                                                                                                                                                                          cGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCGGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGGGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAAGGCCAGTTCATGAACCAGCCGTACCGCGCGTTTCCCCCACATGGCGGCCACGTCCT‑CAAAAACCAACGGTGTTTCCTGGACATCAccgg                                               >  2:130404/1‑236 (MQ=255)
                                                                                                                                                                           ggTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CCAATACCCAGGGTGTTTCCTGGATCTCACCCgg                                              >  1:99277/1‑239 (MQ=255)
                                                                                                                                                                             ttAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCCGGTTGTTTCCTGGATTTCACCCgggt                                            >  1:606517/1‑237 (MQ=255)
                                                                                                                                                                             ttAAGATGCGTAATCGCATCCGGCAAACACCCATACCG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACcgcg                                                                                                             <  2:489844/175‑1 (MQ=255)
                                                                                                                                                                                      gTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCCTGACGCAGGCGTACCCCCCCTTGCACCCCAATCGCGGCCGGGTCTG‑CCAAAAACCAGGATGGTTCCCGGAACCCACCCGGGTCAaccaac                                    >  2:505526/1‑236 (MQ=255)
                                                                                                                                                                                       tCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGGGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGcaccc                                                                                                   >  1:489844/1‑174 (MQ=255)
                                                                                                                                                                                       tCATCGCATCCAGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCCCTGGCGCAGCAGG‑CGGATTTCCTCCTTAAAGTGCGCACTATCATCCCTTAACGTCAGCGCCAACACGTTCCCCAGCCCGTCCAGAGCCAGTTCCTGACGCCGGCGTACCGCGCGCTGCACAACAATGCGTGCATCAGCGT‑CAAACACCCAGGTTGTTTCCGGGAGAACAAGCCGGTCAACCAcccc                                  >  2:587908/1‑239 (MQ=255)
                                                                                                                                                                                               tcccgacaaCACACCTCACG‑GATGGGCTACAAACCGGGAAAGCCGCGGGCGCAGCCGG‑CGGGTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGa                          <  2:99277/232‑1 (MQ=255)
                                                                                                                                                                                               tCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGGTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCCCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGGTTCCTGGATCTTACCCGGTTCCACCCCCCCCAAACCGa                          >  1:286184/1‑239 (MQ=255)
                                                                                                                                                                                                  ggCAAACCCACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACaa                       <  1:815690/239‑1 (MQ=255)
                                                                                                                                                                                                          cacgTCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATttcattc               <  1:587908/239‑1 (MQ=255)
                                                                                                                                                                                                             tccccG‑GGTGAGCTACAAGCCGGGAAAACCGCTGGGCCAGCAGG‑CGGGTTTCCTGCTTCAGGGGCGCAATATCACCCATTACCGTCAGCGCCCCCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAAc            <  2:671723/235‑1 (MQ=255)
                                                                                                                                                                                                              tCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGGTAGCGCCACCCCGATCCCCGGCCCGTCCAAAGCCAGTTACTGACCCAGGCGTACCGCGCGGTGCACCACAATGGCGGCAGTGGCGTCCAAACACCAGGGTTGTTTCC‑GGATCTCACGCGGTGCAACCACCCCCAAACACAAAAATTCATTaaaa            >  2:404599/1‑237 (MQ=255)
                                                                                                                                                                                                              ccgcG‑GATGGGCCACAAGCCGGGAAAGCCGGTGGCGCAGAGGG‑CGGTTTTCCCGCTCCAGGTGCGCAATAACACCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAAct           <  2:597529/236‑1 (MQ=255)
                                                                                                                                                                                                                      tGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCCCGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCCCTTc    >  1:597529/1‑239 (MQ=255)
                                                                                                                                                                                                                        aGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCACGCGGGTCAACCACCCCCAAACCGACAATTTCATTCAACTCCCCTTCag  >  1:805849/1‑239 (MQ=255)
                                                                                                                                                                                                                        aGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCACCCGGTTCAACCAACCCCAAACCCACAAATTCCTTCAAACCCCCTTctc  >  1:43781/1‑237 (MQ=255)
                                                                                                                                                                                                                        aGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCCCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAAGCGTACCGCGCGTTGCACCACAATGGCGGCAAGGTCCT‑CCAATAACCCGGTTGTTTCCTGGATCTCACCCCGGTCAACCCCCCCCCAACCCACAATTTCATTTAACACCCCTTCag  >  2:333529/1‑239 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                               
GGCACTAAGCCAGTTATTGCTGGTGTTAGTGGAGTGCGACCACCAAAACAGATAACGTCGTTTGTGCGCCTGAAAAGACGCGTTTAGCGTCGCATCAGGCATTATGGCGCAGTTGCCGGATGCGGCGTGAACGTCTTATCCGGCCCACAGGAACTGTAATCTTTGTAGACCGGTTAAGATGCGTCATCGCATCCGGCAAACACACATCACG‑GATGAGCTACAAACCGGGAAAGCCGCTGGCGCAGCAGG‑CGGTTTTCCTGCTTCAGGTGCGCAATATCATCCATTAACGTCAGCGCCACCGCGATCCCCGGCCAGTCCAGAGCCAGTTCATGACGCAGGCGTACCGCGCGTTGCACCACAATGGCGGCATGGTCGT‑CAAATACCCAGGTTGTTTCCTGGATCTCACGCGGTTCAACCACCCCCAAACCGACAATTTCATTCAACTCCTCTTCAG  >  NC_000913/1062344‑1062797

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