Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,299 A→C 19.0% V418G (GTG→GGG)  aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913120,2990AC19.0% 155.7 / 2.6 69V418G (GTG→GGG) aroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base A (45/8);  new base C (2/14);  total (47/22)
Fisher's exact test for biased strand distribution p-value = 1.85e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTG‑TACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGACAC  >  NC_000913/120077‑120537
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTG‑TACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAACCGAGACCTAACACGATCAGCCATACCGGGATCAGGGATACCGAAATCGCCCTTCCTGGGGTCATCAG‑CATAATACACAGAACCGCCGCCaa                                                                                                                                                                                                                                 >  1:757541/1‑238 (MQ=255)
                      gATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTG‑TACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAACCGACAACTAACACGATCAGCCATACCGTGATCAGGTATACCGAAATAGCCATTACCGGGGGCATCAG‑CATAATACACAGTACCCCCCCCCAAAAcaaccagcgaagccagtta                                                                                                                                                                                                           >  1:325121/1‑222 (MQ=255)
                         tGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTG‑TACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCCGGTATACCGAAATCGCCATTCCTTGGGTCCTCAG‑CATAATCACCAGTTCCGCCGCCATAAATAGCAAGCAGATCCAGGTAcca                                                                                                                                                                                                        >  2:303724/1‑238 (MQ=255)
                          gctctactGACGATTTGCGGGGCGTTGGTGTAAAGATTATTGCCCTCCCCCTG‑TACGGGTGAGGGCGTAGAGAGATTAATGCTCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATCCCGGGATCAGGTATACCGAAATCGCCATCCCTGGGGTAATCAG‑CATAATCCCCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCa                                                                                                                                                                                                       <  2:251402/236‑1 (MQ=255)
                          gCACTAGTGGCGTTTGACGGGGCGTGGGTGTAAAGATTATTGCCCTCACCGTG‑TACGGGTGAGGGCGTAGAGATATTAAGGCGCTTTTCCG‑GTTTTGGCGGTTTTCTTTTTAAACAGATAGCCGATCCCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAAGCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCa                                                                                                                                                                                                       <  2:645278/239‑1 (MQ=255)
                          ccACTACTGAGGTTTTGCGGGGCTTTGGTGTAAAGTTTATTGCCCTCACCTTG‑TAGGGGTGAGGGCGTAGAGAGATTAATGCGTTTTTACG‑GCTTGGGCGGTTTTCTCTTTAAACAGATAGCCGATCCCTAACACGATCAGCCATCCCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCa                                                                                                                                                                                                       <  2:757541/238‑1 (MQ=255)
                            actactGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTG‑TACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGTCGATACCTAACACGATCAGCCATACCGGGATCAGGGATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGAAGATCCAGTTACCCAGc                                                                                                                                                                                                     >  2:30986/1‑239 (MQ=255)
                                               ttaaaagaaacaggATTATTGCCCT‑TCCCTGTTACGGGTGAGGGCGTAGAGAAATTAAGGCGTTTTTACG‑GTTTGGGCGGTTTTCCTTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATCCCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCCCCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑Ga                                                                                                                                                                                 <  2:651495/226‑1 (MQ=255)
                                                 gTTGGTGTAAAGATTATTGCCCTCACCCTG‑TACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCCTAAACAGCAGGCAGATCCAGTTACCCACCGGATAAAGCAGAGCAGG‑Gaa                                                                                                                                                                                >  2:206606/1‑238 (MQ=255)
                                                                  tGCCCTCACCCTG‑TACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTAccg                                                                                                                                                                                                                                        <  2:27244/166‑1 (MQ=255)
                                                                       tCACCCTG‑TACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CCTAATCACCCGAACCGCCGCCCTAAACAACCGGCAGATCCAGTTAACCAGCCGATAAAACAGAGCAGG‑GAAAGGAGGTACCCCggcccttaaac                                                                                                                                                         >  1:651495/1‑229 (MQ=255)
                                                                         aCCCTG‑TACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAATTACCCAGCGGATAAAGCAAGGAAGG‑GAAGGGAGTTACACCGCCTTGGTCCt                                                                                                                                                         >  2:757456/1‑237 (MQ=255)
                                                                              g‑TACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATAACGAAATCGCCATTCCCGGGGTCCTCAA‑CATAATAACCAGTAACGCCGCCCTAAAAAGCAGGCAGAACACGTTTCCCCACGGGTAAAACACAACAGG‑GAAACGAGTTACACCGCCttttccacgacttgg                                                                                                                                                  >  1:492984/1‑226 (MQ=255)
                                                                              g‑TACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAAAAGGTAGGCGATACCTAACACGATCAGCCATACCGGGATCAGGGATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCCTAAACAGCAGCCAGATCCAGTTACCCAGCGGCTAAAGCAGAAACAG‑GAAGAGAGTTACCACGCCATGTTCCCGCCTGgc                                                                                                                                                  >  1:100436/1‑239 (MQ=255)
                                                                                tACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCCGATCCAGGTACCCCGCGGATAAAGCAAAGCAGG‑GAAGCGAGTTACCACGCCCTggttccgactggg                                                                                                                                                  >  2:242908/1‑226 (MQ=255)
                                                                                tACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAACCGATACCTAAAACGATCAGCCATACCGGGATCAGGGATACCGAAAACGCCATTCCTGGGGTCATACG‑CCTAATCACCAGTACCCCCCCCCCAAACAGCAGGCAGAGCCAGTTACCCAGCGGAAAAAGCAAGGCGGG‑GAAACGAGGTACCACGCCTTTTTCCCGCTTgggg                                                                                                                                                 >  1:932764/1‑237 (MQ=255)
                                                                                     gtggggggttAGAGATATTAATGCTTTTTTCCG‑GCTTTGGGGTTTTTCCTTTAAAACAGATGGCCGATACTTAACACGATCAGCCATACCGGGATCAGTTATCCCGAAATCGCCATCCCGGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACg                                                                                                                                                                   <  2:575705/207‑1 (MQ=255)
                                                                                          ggCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTTTCTTTAAACAAATAGCCGATACCTAACACGATCAGCCATAACGGGATCAGGTATACCGAAATCGCCATTCCTGGGGGCATCAG‑CCTAATCACCAGTACCCCCGCCATAAACAGCAGGCCGATCCAGTTACCCCCCGGATACACCCAAGCAGG‑GAAACGCAGTACCAACCCTTTTTCCCGCTTGGggcgcagagaat                                                                                                                                       >  1:286062/1‑231 (MQ=255)
                                                                                          ggCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAag                                                                                                                                                                                             >  2:201220/1‑186 (MQ=255)
                                                                                                 ggggtattaTGCCCTTTTTCGTGTTTGGGCGTTTTCCTTTTAAAACAGATCGCGGACCCCAAACACGACCACACATCCGTGGAACAGGTATCCCGAATCCGCCATTCGTGGGGTCATCAG‑CATAATCCCCAGTCCCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCata                                                                                                                                  <  1:826819/230‑1 (MQ=255)
                                                                                                   aaaTTAATGGTTTTTTCCG‑GTTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATCCCGACCACGACCAGCCAACCGGGGAACAGGTAGCCCAAAATCGCCTTTCCGGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATg                                                                                                                                <  2:492984/235‑1 (MQ=255)
                                                                                                     aTTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGGCATCAG‑CATAATCAACAGTACCGCCGCCATAAACAGAAAGCAGATCCATGTACCCAGCGGGAAAAGCAGAGCCGG‑GGAACGAGATACCACGCCTTGTTCCTGCTTGGCGCGCCGGAATTTCATATgaccc                                                                                                                            >  1:555294/1‑235 (MQ=255)
                                                                                                          tGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCCTATGCGCCCGGCt                                                                                                                       >  2:127285/1‑239 (MQ=255)
                                                                                                             gcTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCGGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGAGAGTTACCACGACTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAat                                                                                                                    >  2:124866/1‑239 (MQ=255)
                                                                                                                     g‑gCTTTGGCGTTTTTCTTTTAAAACAGATAGCCGATACTTAACCCGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCTTTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGccc                                                                                                            <  2:100436/239‑1 (MQ=255)
                                                                                                                      tataaagtctcttttCTTTTTAAACAGATCCCCTATCCCAAACGCGATCCGCCACCCCGGGAACGGGTATCCCGAACTCCCCATCCCGGGGATCAACAG‑CATAACCCCAAGTCCCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGccc                                                                                                            <  1:436591/228‑1 (MQ=255)
                                                                                                                      ggttttGGCGGTTTTCTCTTAAAACAGATAGCCGATCCCTAAACCGATCAGCCATCCCGGGATAGGGTAGACCGAACTCGCAATTCCGGGGGTCATCAG‑CATAATCCCCAGTCCCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGccc                                                                                                            <  2:286062/236‑1 (MQ=255)
                                                                                                                          ttGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGCAGTACCACGCCCTGTTCCTGCTTTGCCCGACGGAATTTCCAATGCGCCAGGCTAATCATCACCCCGtt                                                                                                        >  2:501665/1‑239 (MQ=255)
                                                                                                                          ttGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTACTGGGGTCATCAG‑CATAATCACCCGTACCGCCGCCATAAACAGCAGGCAGATCCAGCTAACCAACGGGTAAAGCAGCACAGGCAAAGCG‑GTTACCACACCTTGTTCCTGCTTATGGCGACCGAATTTCAAAGGGGCCCGCCCAAACACCCCCCAGtt                                                                                                        >  1:182079/1‑239 (MQ=255)
                                                                                                                            ggCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAag                                                                                                                                                                                             >  2:502604/1‑153 (MQ=255)
                                                                                                                                gTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATaaa                                                                                                                                                                                                                              >  2:427632/1‑116 (MQ=255)
                                                                                                                                 ttttttttttAAACAGATTGCCTAACCCAAAACCGACCCCCCAACCCGTGATCGGGTAAACCAAACCCGCCTTTCGGGGGGTCATCAG‑AATAACCACAAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTAcc                                                                                                 <  1:30986/239‑1 (MQ=255)
                                                                                                                                             aCAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTAc                                                                                                  >  2:436591/1‑226 (MQ=255)
                                                                                                                                                         tACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGTTTACCCGGGCAGATACCACCACCGCCAtt                                                                         >  2:932555/1‑239 (MQ=255)
                                                                                                                                                               acacGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCCTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAg                                                                    >  2:469805/1‑238 (MQ=255)
                                                                                                                                                                     tCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGaaa                                                             >  2:537578/1‑239 (MQ=255)
                                                                                                                                                                     tCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCAGCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCCTTACCAGTCCGcaa                                                             >  2:976805/1‑239 (MQ=255)
                                                                                                                                                                        gCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGGGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCg                                                          >  2:51677/1‑239 (MQ=255)
                                                                                                                                                                         cattACCGGAACCAGGTGACCCGAAACCGCCA‑TCCCGGGGTTAACAGAAATAATCCCCAGTCCCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCg                                                          <  2:113857/235‑1 (MQ=255)
                                                                                                                                                                             aCCGGGCTCAGGTATACCGAAATCGCAATTCCTGGGGTCATCCG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATACAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGAGCGACGGAATTTCATATGCGCCAGGACAATCACCGCCCCGTTGATTACCAGTGACATTACAACCAGCGCCATTACAAGGCCGAAAGAggcccc                                                     >  1:724519/1‑235 (MQ=255)
                                                                                                                                                                             aCCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGActc                                                     >  2:614817/1‑239 (MQ=255)
                                                                                                                                                                             aCCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTCCAGTTACCACCCGCGCCATTAACAGTCCGAAAACGGAccc                                                     >  1:861533/1‑237 (MQ=255)
                                                                                                                                                                                gggATCAGGTATCCCGAAACCCCCATTCCGGGGGTCATCAG‑CATAACCCCCAGTCCCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAag                                                                                                                                                                                             <  1:358138/101‑1 (MQ=255)
                                                                                                                                                                                   aTCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCTGATCCAGTTACCCAGCGGATAAAGCAGCGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATCTGCCAGCCTAATCATCGCCCCGTTGATTCCCAGTGCCATTACAACCAGCGCCCTTAACacaccaaaacagatcccctggg                                                >  1:113857/1‑217 (MQ=255)
                                                                                                                                                                                   aTCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACAAGGGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACCCTGGGGc                                               >  2:756620/1‑239 (MQ=255)
                                                                                                                                                                                   aTCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCCAGTTACCAAGCCTTGTTCCTGCTTGGCGCCACGGAATTTCCTATGAGCCAGGCTAATCATCGCCCAGTTGCTTACCAATGCCGCAACCACCCGCCCCCATAACAATACGAAAACGGCCACTGGGGc                                               >  1:70826/1‑239 (MQ=255)
                                                                                                                                                                                   aTCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAAGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTAACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGGGCCAGGCTAATAATCGCCAAGTTGATTACCAGGGCAGTAACACCCATGCCCATTAACAGACCGAAAGGGAGCCCTGGGGc                                               >  1:771901/1‑239 (MQ=255)
                                                                                                                                                                                     cAGGTATCCCAAACCCGCCATTCTGGGGGCAATCAG‑AATAACCCCCATCACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCaa                                             <  2:555294/239‑1 (MQ=255)
                                                                                                                                                                                      aGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCCGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAACGCCATTAACAGTCCGAAAGCGGACCCTGTGGCAAg                                            >  2:353119/1‑239 (MQ=255)
                                                                                                                                                                                      aGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGAGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAACGCCCTTAACAACCCGAAAAGGGACCCTGGGGaaaa                                            >  1:309909/1‑238 (MQ=255)
                                                                                                                                                                                        gtaccccGAAATCGCCATCCCTGGGGCCCTCAG‑CAAAATCCCCAGTCCCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGa                                                                                                                                              <  1:288974/134‑1 (MQ=255)
                                                                                                                                                                                         tataCCGAAATCGCCATCCCGGGGGTCATCAG‑CAAAACCCCCAGTCCCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAgg                                           <  2:325121/237‑1 (MQ=255)
                                                                                                                                                                                          ataCCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAgg‑g                                                                                                                                                                                  >  2:358138/1‑101 (MQ=255)
                                                                                                                                                                                           taCCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGt                                       >  2:365389/1‑239 (MQ=255)
                                                                                                                                                                                            aCCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTAcc                                                                                                                                                                                                         >  2:822257/1‑77 (MQ=255)
                                                                                                                                                                                              caaaacccccATTCCGGGGTTCACAAC‑CAAAACCCCAAGTCCCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTaa                                    <  2:932764/231‑1 (MQ=255)
                                                                                                                                                                                                aaaTCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAata                                                                                                                   >  2:5278/1‑157 (MQ=255)
                                                                                                                                                                                                aaaTCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGAt                                                                                        >  2:762274/1‑185 (MQ=255)
                                                                                                                                                                                                aaaTCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGAAGGTaacc                                  >  1:149425/1‑237 (MQ=255)
                                                                                                                                                                                                aaaTCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACAACCAGCGCCATGAACAGTCCGAAAGCGGACTCTGGGGCAAGGGAGTTAATc                                  >  2:845504/1‑239 (MQ=255)
                                                                                                                                                                                                  aTCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGa                                                                                                                                                                                                  >  2:966259/1‑78 (MQ=255)
                                                                                                                                                                                                  aTCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCa                                                                                                                                    >  2:288974/1‑139 (MQ=255)
                                                                                                                                                                                                     gCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTAc                                                                                                  >  2:997217/1‑170 (MQ=255)
                                                                                                                                                                                                          tcgggggcactcgaG‑CATAACCCCCAGTCCCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCTGATACCGCCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGcaca                        <  2:182079/226‑1 (MQ=255)
                                                                                                                                                                                                          tCCTGGGGTCATCAG‑CATAACCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGcaca                        <  2:861533/239‑1 (MQ=255)
                                                                                                                                                                                                          gccgggggACATCAA‑CAAAACCCCCAGTCCCGCCGCCATAAACAGCAAGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGcaca                        <  2:771901/235‑1 (MQ=255)
                                                                                                                                                                                                          cccTGGGGACCACAG‑CAAATCCCCCGGTACCGCCGCCATAAACAGCAGGCAAGTCCAGTTACCCAGCGTATAAAGCAGAGCAGG‑GAAGCGAGTGACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGcaca                        <  2:452704/238‑1 (MQ=255)
                                                                                                                                                                                                                    atcaG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GCAGCGGGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGTAGATACACCCCGCGCCATTATCAGTCCGAAACCGACCACTGGCGAAAGGCCGTTAAACATTAAGccaaaccacgtaac              >  1:383301/1‑226 (MQ=255)
                                                                                                                                                                                                                    atcaG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAACCGGACCCTgggg                                                >  2:645106/1‑205 (MQ=255)
                                                                                                                                                                                                                                 aCCAGACCCGCGCCCATAAACAGCGGGCAGATCCAGTTGCCCAGCGGATAAAGAAGAGCAGG‑GAAGGGATTTGCCCCGCCTTGTTCCTGGTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTCCCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGTACTCTGGGGCAAGTTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGacac  <  1:303724/239‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
CCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTG‑TACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACG‑GCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAG‑CATAATCACCAGTACCGCCGCCATAAACAGCAGGCAGATCCAGTTACCCAGCGGATAAAGCAGAGCAGG‑GAAGCGAGTTACCACGCCTTGTTCCTGCTTGGCGCGACGGAATTTCATATGCGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGACAC  >  NC_000913/120077‑120537

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