Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 900,209 C→A 8.5% D123Y (GAT→TAT)  artJ ← arginine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913900,2090CA 8.5% 156.6 / 4.9 51D123Y (GAT→TAT) artJarginine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (20/26);  new base A (2/3);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

CGTTATTCAGTTTTTCCAGCAGGGCTTTGTTATCCGGACGTACAGCGATGCCCAGGCCGGTGCCAAAATATTGCGGATCGGTCACTTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCAAAGCCGACAATCTCATTATTAGCACCTATAGATTCAAAGGGTGGATAGGTGGCTGAAACGCCAAAATT  >  NC_000913/899909‑900506
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                         
cgTTATTCAGTTTTTCCAGCAGGGCTTTGTTATCCGGACGTACAGCGATGCCCAGGCCGGTGCCACAATATTGCGGATCGGTCACTTTCTCAGTAGCAACACCCAGTTGTGGATTGGGTTTCCGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATCCGACCATTTTTCAGCTCGCTAAAGGCATTCTGATCACTGGCCTACGAGCCAGTTTTCACTTCCGGGTGCTCATCCTGCATATATATCTGGAGCATAGAACCGTTTTCCATCCCAATCCGTATCCCTTTCAGATc                                                                                                                                                                                                                                                                                                           >  1:429892/1‑301 (MQ=255)
     ttCAGTTTTTCCAGCAGGGCTTTGTTATCCGGACGTACAGCGATGCCCAGGCCGGTGCCAAAATATTGCGGATCGGTCAATTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCACATACCACATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATAACTTTCTGGTGCTTAGTACCGTTTTCCATCCCAAAACGTTACCCTTCCAAACAGGCaa                                                                                                                                                                                                                                                                                                      >  1:259496/1‑301 (MQ=255)
             ttCCAGCAGGGCTTTGTTATCCGGACGTACAGCGATGCCCAGGCCGGTGCCAAAATATTGCGGATCGGTCACTTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAAAACGTTTGCCGATCTGAAAGGCAAACGTATTg                                                                                                                                                                                                                                                                                              <  2:584753/301‑1 (MQ=255)
                cagcagGGTTTTTTTATCCGGTCGTTCAGCGATGCCCAGGCCGGTGCCAAAATATTGCGGATCGGTCACTTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAAAACGTTTGCCGATCTGAAAGGCAAACGTATTggga                                                                                                                                                                                                                                                                                           <  1:88590/301‑4 (MQ=255)
                                       gTACAGCGATGCCCAGGCCGGTGCCAAAATATTGCGGATCGGTCACTTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCacgac                                                                                                                                                                                                                                                                    <  1:66176/301‑1 (MQ=255)
                                        tACAGCGATGCCCAGGCCGGTGCCAAAATATTGCGGATCGGTCACTTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTATTGTAGGTATATGTTTTGGCAATCacgacg                                                                                                                                                                                                                                                                   >  2:368489/1‑301 (MQ=255)
                                            gCGATGCCCAGGCCGGTGCCAAAATATTGCGGATCGGTCACTTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTg                                                                                                                                                                                                                                                               <  2:526317/301‑1 (MQ=255)
                                            gCGATGCCCAGGCCGGTGCCAAAATATTGCGGATCGGTCACTTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTg                                                                                                                                                                                                                                                               >  2:625153/1‑301 (MQ=255)
                                                  cccAGGCCGGTGCCAAAATATTGCGGATCGGTCACTTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTTGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTAAGtttt                                                                                                                                                                                                                                                         >  1:591075/1‑301 (MQ=255)
                                                             gCCAAAATATTGCGGATCGGTCACTTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAAGGc                                                                                                                                                                                                                                              >  1:460965/1‑301 (MQ=255)
                                                               aaaatgtttGCGGTTTGGTCATTTTCTCAGTAGCAACACTCAGTTGTTGATTGGTTTTTAGCCATTCGTTTACCACCGCTGTGTCACCAAGTACCCCATCAATACGACCATGTGTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGt                                                                                                                                                                                                                                            <  1:588893/294‑1 (MQ=255)
                                                                aaaaTATGTCGGATCGGTCACTTTCTCAGTATCAACCCCAAGTTGCGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGGTTTCACTTCCGGGTGCTGATCCTGAATATATTGCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTg                                                                                                                                                                                                                                           <  1:576443/301‑1 (MQ=255)
                                                                       tGCGGATCGGTCACTTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACg                                                                                                                                                                                                                                    <  2:461529/301‑1 (MQ=255)
                                                                         cGGATCGGTCACTTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGAt                                                                                                                                                                                                                                  <  2:591836/301‑1 (MQ=255)
                                                                                          cAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAAAGATACCTGTTTGCTACGccc                                                                                                                                                                                                                 >  1:42012/1‑299 (MQ=255)
                                                                                                                    tttttTAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTa                                                                                                                                                                                       <  2:605069/300‑1 (MQ=255)
                                                                                                                             cATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCaa                                                                                                                                                                              >  2:231187/1‑300 (MQ=255)
                                                                                                                                    ttACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGACAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGATCCGGGGTGATATCAAAACCGGAGATTACGGCGTCATATTTTCt                                                                                                                                                                       >  1:509863/1‑301 (MQ=255)
                                                                                                                                     ggcttttccgatctCACCAAATACCCCATCTATACGACCTTTTTGTAGATCGATAAAGGCATTCTGATATCTGTCATAAGAGACAGTTTTCACTTCCGGTTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTg                                                                                                                                                                      <  1:370479/288‑1 (MQ=255)
                                                                                                                                      cccaccGCGGTGTCACCAAATACCCCATCAATACGACCTTTTTTCAGTTCGATAAAGGCATTCTGATAACTGTCATAAGGGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTGCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGa                                                                                                                                                                     <  1:624294/300‑1 (MQ=255)
                                                                                                                                             gtttgtCACCAAATACCCCATCAATACGACCATTTTTCAGTTCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTGGCCTTTCAGATCGGCAAACGTGTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAg                                                                                                                                                              <  1:465291/298‑1 (MQ=255)
                                                                                                                                             ctttgtacacacaTACCCCCTCAATACGACCATTTTTCAGGTCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTGTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAg                                                                                                                                                              <  1:368005/289‑1 (MQ=255)
                                                                                                                                              tttgtCACCAAATACCCCATTAATTCGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATATGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAgg                                                                                                                                                             <  1:230944/299‑1 (MQ=255)
                                                                                                                                                   cACCAAATACCCCATCAATACGTCCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACgg                                                                                                                                                        <  2:42065/301‑1 (MQ=255)
                                                                                                                                                   cACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGagatcggaagagca                                                                                                                                                        >  2:370975/1‑288 (MQ=255)
                                                                                                                                                                       cGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAAcgcgg                                                                                                                                    >  2:465563/1‑300 (MQ=255)
                                                                                                                                                                        tttccgatctCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTg                                                                                                                                   <  2:195310/292‑1 (MQ=255)
                                                                                                                                                                                  cAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGCGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGagatcggaag                                                                                                                         >  1:195067/1‑292 (MQ=255)
                                                                                                                                                                                    gATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTAATTAGTAAAAGt                                                                                                                       >  1:562876/1‑301 (MQ=255)
                                                                                                                                                                                            aaGGCATTCTGATAACTGTCATAAGGGATAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTctg                                                                                                               <  1:464999/301‑1 (MQ=255)
                                                                                                                                                                                            aaGGCATTCTGATAACTGTCATAAGAGACAGTTTTCCCTTCCGGGTGCTGATCCTGACTATATTTCTGGTGCGTAGTACCGTTTTCCCTCACAATACGTTTGCCTTCCAGATCAGCCACCGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTAAGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGATGTCGACCGCGGGGTTAGTAAAAgataatttct                                                                                                               >  1:326526/1‑292 (MQ=255)
                                                                                                                                                                                               gCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTctgcct                                                                                                            >  2:463861/1‑301 (MQ=255)
                                                                                                                                                                                                          attGTCATAAGTGACAGTTTTCTCTTCCGGGTGCTGATCCTGAATATATTGCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCAtttgtttg                                                                                                 <  2:556323/299‑1 (MQ=255)
                                                                                                                                                                                                          attGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCAtttgtttg                                                                                                 <  2:259764/299‑1 (MQ=255)
                                                                                                                                                                                                          aCTTTCATAAGAGACAGTTTTCACTTCCGGTTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATGTCAGGGACGGGATCAGGCTGTCGAACGCGTGATGAGTAAAAGTACATTCTGCCTGCAtttgtttg                                                                                                 <  2:326930/301‑1 (MQ=255)
                                                                                                                                                                                                           cTGTCATAAGGGACAGGTTTCACTTCCGGGTGCTGTTCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGc                                                                                                <  1:55637/301‑1 (MQ=255)
                                                                                                                                                                                                                       aCAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTg                                                                                    <  2:640519/301‑1 (MQ=255)
                                                                                                                                                                                                                                     cGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCa                                                                      <  2:367125/301‑1 (MQ=255)
                                                                                                                                                                                                                                                  tGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCCGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATAGATATCAAAGCCGACAAtcc                                                         >  2:331550/1‑300 (MQ=255)
                                                                                                                                                                                                                                                        tatTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATAAAAGCCGACAATCTCattaa                                                   >  2:208793/1‑300 (MQ=255)
                                                                                                                                                                                                                                                          tggagttcagacgtgtgctcttccgaTCTCAAAACGTTTGCCGATCTGAAAGGCAAACGTATTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCAAAGCCGACAATCTCAttatta                                                 <  2:170027/276‑1 (MQ=255)
                                                                                                                                                                                                                                                                gTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCTATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCAAAGCCGACAATCTCATTAATAGCACCt                                           >  1:419698/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                   cttagtaCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCTATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCAAAGCCGACAATCTCATTATTAGCACCtata                                        <  2:420261/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                      agtaCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGTAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGCTATCAAAGCCGACAATCTCATTATTAGCACCTATAGAt                                     >  1:374724/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                       gtaCCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCAAAGCCGACAATACCATTATTAGCACCTAAAGaat                                    >  1:187073/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                         tccGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCAAAGCCGACAATCTCATTATTAGCACCTATAGATTCa                                  <  2:187315/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                ttcttCCCAATTCGTTTGCCTTTCAGTTCGTCATACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCAAAGCCGACAATCTCATTATTAGCACCTATAGATTCAAAGGGTg                           <  1:208550/297‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       cAAAACGTTTGCCGATCTGAAAGGCAAACGTATTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCAAAGCCGACAATCTCAttattaagatcggaagagcgtcgtgtagggaaaga                    >  1:169830/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                                            tgTTTGCCTTTCAGATCGGCAAACGTTTTGTATGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTATACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCAAAGCCGACAATCTCATTATTAGCACCTATAGATTCAAAGGGTGGATAGGTGGCTg               <  2:430456/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                     ttCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCAAAGCCGACAATCTCATTATTAGCACCTATAGATTCAAAGGGTGGATAGGAGGCTGAAACGCCaa      >  1:153494/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                         gATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCAAAGCCGACAATCTCATTATTAGCACCTATAGATTCAAAGGGTGGATAGGTGGCTGAAACGCCAAAAtt  <  2:563590/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                         gATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATAAAAGCCGACAATCTCATTATTAGCACCTAACAATTCAAAGGGAGGAAAGTAGGCAAAAAAGCAAAAAta  >  2:149136/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                         
CGTTATTCAGTTTTTCCAGCAGGGCTTTGTTATCCGGACGTACAGCGATGCCCAGGCCGGTGCCAAAATATTGCGGATCGGTCACTTTCTCAGTAGCAACACCAAGTTGTGGATTGGTTTTCAGCCATTCGTTTACCACCGCTGTGTCACCAAATACCCCATCAATACGACCATTTTTCAGATCGATAAAGGCATTCTGATAACTGTCATAAGAGACAGTTTTCACTTCCGGGTGCTGATCCTGAATATATTTCTGGTGCGTAGTACCGTTTTCCATCCCAATACGTTTGCCTTTCAGATCGGCAAACGTTTTGTAGGTATCTTTTTTGGCAATCACGACGGCTGAGTTTTCATAGTAGGGCGTGGTAAACGATACCTGTTTGCTACGCTCCGGGGTGATATCCATACCGGAGATTACGGCGTCATATTTTCTGAATTTCAGGGACGGGATCAGGCTGTCGAACGCGTGATTAGTAAAAGTACATTCTGCCTGCATTTGTTTGCACAAGGCTTTTGCCAGATCGATATCAAAGCCGACAATCTCATTATTAGCACCTATAGATTCAAAGGGTGGATAGGTGGCTGAAACGCCAAAATT  >  NC_000913/899909‑900506

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: