Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,474 T→C Q162Q (CAA→CAG yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,4740TC98.0% 149.1 / NA 51Q162Q (CAA→CAGyaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (30/20);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGC  >  NC_000913/7232‑7717
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
ccaAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTcc                                                                                                                                                                                                                                             <  1:57250/251‑1 (MQ=255)
ccaAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTcc                                                                                                                                                                                                                                             <  1:417433/251‑1 (MQ=255)
   aaaGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTGCGCCTGAATTCCGCt                                                                                                                                                                                                                                          <  1:167891/251‑1 (MQ=255)
                       ccAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGccgta                                                                                                                                                                                                                           <  1:496671/246‑4 (MQ=255)
                       ccAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTAAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGCGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTAACCAGCGGGGGAAAATCAAAAGAAAAACTAAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGccgtaagaac                                                                                                                                                                                                                      >  2:496671/1‑243 (MQ=255)
                             gACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCa                                                                                                                                                                                                                <  1:265600/251‑1 (MQ=255)
                              aCGTGGGGAAGTTGACCGATATTCATTACGCAACTCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGAGACGCCATGAAGACCGCGAGTGATCGCAAGCAGAGCAAAGAAAGCGCGAATAATACCTGTCACAAGCGGGGGAAAATCAAACAAAACACTCAGGGCTCGCGCAACGGCGTTAGCCTGAATTCAGCAGCAAATTAttccgaaggcactgagca                                                                                                                                                                                                                >  2:334784/1‑233 (MQ=255)
                                   gggAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGa                                                                                                                                                                                                          <  1:153053/251‑1 (MQ=255)
                                   gggAAGTTGCCCGAGAGTCAATACGCAAATCCCCAGGCTAGTCAGTACCCAGAATAGCGCCACCAACCGGACAAACCCCTGCATGAGCCGGGCGACGCCAAGAAGACCGCGAGTGATCGCCAGCAAAGCAAAGTCAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACGCACTGATCGCGAAACGGCGTCCGCCTGAAGTACGCTGAAAAGTATGCCGTAGGGGTTGAGCAAAAAGa                                                                                                                                                                                                          <  2:55406/251‑1 (MQ=255)
                                         ttGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCg                                                                                                                                                                                                     >  1:450525/1‑250 (MQ=255)
                                          tGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTAGCCTGAATTCCGCTGCAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGa                                                                                                                                                                                                    >  1:628505/1‑250 (MQ=255)
                                                                          ggTCAGTACCCAGATGATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGccc                                                                                                                                                                     <  1:349658/249‑1 (MQ=255)
                                                                                  cccAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCat                                                                                                                                                            >  1:672594/1‑250 (MQ=255)
                                                                                    cAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCctat                                                                                                                                                          >  2:603899/1‑250 (MQ=255)
                                                                                      aattctCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATAcc                                                                                                                                                       <  2:191002/246‑1 (MQ=255)
                                                                                          aTCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGc                                                                                                                                                   >  1:629151/1‑251 (MQ=255)
                                                                                              ccATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATCTACCATGCCgtc                                                                                                                                               >  2:144357/1‑249 (MQ=255)
                                                                                                 ccAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCg                                                                                                                                            <  1:591218/250‑1 (MQ=255)
                                                                                                     cGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCAc                                                                                                                                        >  1:544714/1‑251 (MQ=255)
                                                                                                          cAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGGCCGCCACGAAAc                                                                                                                                   >  1:21169/1‑251 (MQ=255)
                                                                                                           aaacccctgcatggtccagtcGACGCCATGAAGACCGGGAGTAAACGCCAGCAGAGCAAAGACATCGAGAATAATACCGGTCCACATCGGGGGAAAATTAAAAGAAAAACTCAGGTCTCGCGCAACGGCGTTCGCCTGAATCCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATTCCATCCCGGTCCGCCATGAAACt                                                                                                                                  <  2:676757/231‑1 (MQ=255)
                                                                                                              ccccTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGcc                                                                                                                               >  2:234640/2‑251 (MQ=255)
                                                                                                               cccTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGccc                                                                                                                              <  1:234640/251‑1 (MQ=255)
                                                                                                                          cctggCGACGCCATGAAGCCCGGCATAGATCGCTAGCAGAGCAAAGATAGCGTGAAGAATTCCTGTCACCAGCGGGAGAACATTAAAAGAAAAACGCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCGGAAAATTATGCCGTAGTCGTTGAGCAAAAAGACGGCGAACAGACCGCCCATCCACCGCATCCCCAGCCCGCGCGCCATATACCATGCTGTTCCGCCTCGAAACTGTCCAtgggcagctc                                                                                                                   <  2:14699/248‑10 (MQ=255)
                                                                                                                           cGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACCGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACg                                                                                                                  >  2:571211/1‑251 (MQ=255)
                                                                                                                                          aGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCACCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAgctg                                                                                                   <  1:148519/251‑4 (MQ=255)
                                                                                                                                                      atcgCCAGCAGAGCAAAGACAGCGAGAATAATACCTGGCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCAGCACTTTCTTTATAAAGCTGTGCGAGGGAACa                                                                                       <  1:76969/248‑1 (MQ=255)
                                                                                                                                                        cgCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACAt                                                                                      <  1:563648/249‑1 (MQ=255)
                                                                                                                                                         gCCAGCAGAGCAAAGACAGCGAGAATAATACCTATCACCAGCGGGGGAAAATCAACAGAAAAACTCAGGGCTCGCGCAACGGCGATCGCCTGAATTCCGCTGAAAATTATGACGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATACCCAGCCCGCGCGCCATATACCATGCCGTTCCGACACGCAACTGCCCATTGGCATACCGATTTTTATAAAGCTGTGCGAGGGAACATTc                                                                                    >  2:676572/1‑251 (MQ=255)
                                                                                                                                                           cagcagAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCgg                                                                                  >  1:274353/1‑251 (MQ=255)
                                                                                                                                                                 agcaaaGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGCCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCaaa                                                                              >  1:283527/4‑249 (MQ=255)
                                                                                                                                                                   caaaGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACg                                                                            >  1:444787/2‑249 (MQ=255)
                                                                                                                                                                    aaaGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTggcatc                                                                                                                     >  2:234920/1‑202 (MQ=255)
                                                                                                                                                                    aaaGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTggcatc                                                                                                                     <  1:234920/207‑6 (MQ=255)
                                                                                                                                                                    aaaGACAGCGAGAATAATAACTGTCAACAGCGGGGGAAAGTCAAAAGAAAAACTCAGGGATCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAAAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCAATTGGCATCACGTTATTTATAAAGCTGTGCGAGGGAACCTTCGGCAAACGAgg                                                                         >  1:588138/1‑251 (MQ=255)
                                                                                                                                                                          aGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGcc                                                                    >  1:675250/1‑250 (MQ=255)
                                                                                                                                                                              agaATAATACCTGTCAACAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGAGCAACGGCGGTCGCCTGAATTCCGCTGAAAATTATGCAGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCAGCCACGACACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGCACATTCGGAAAACGAGGTCGcactgcc                                                               >  2:468907/1‑245 (MQ=255)
                                                                                                                                                                                 ataataCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCACGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCCTATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGCGGTCGCCATTCCGa                                                             >  2:330703/1‑250 (MQ=255)
                                                                                                                                                                                   aataccTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATa                                                           >  1:617306/5‑250 (MQ=255)
                                                                                                                                                                                      accTGTCACCAGCGGGGGAAAATCAAAAGGAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAGTTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGATGAGCGAGGGAACATTCGGCAATCGAGGCCGCCATGCCGATAAAcg                                                       >  2:46497/2‑251 (MQ=255)
                                                                                                                                                                                         ggTCACCAGCGGGGGAAAAGCAAAAGAAAAACGCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCAGGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCgg                                                    <  1:86498/250‑1 (MQ=255)
                                                                                                                                                                                           tCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGc                                                   >  1:140285/1‑250 (MQ=255)
                                                                                                                                                                                           tCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGAAAACGAGGGCGCCCTGCCGATAAACGCGGc                                                   >  2:508498/1‑250 (MQ=255)
                                                                                                                                                                                                   gggggAACATCAAACGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGCAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCCTATACCATGCCGGTCCGCCCCGCAACTGCCCCTGGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGACAAACGAGGACGCCCTGCCGATAAACGCGGCCACCcacac                                          >  2:154775/1‑250 (MQ=255)
                                                                                                                                                                                                        aaactCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCCCGCAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGCAGCCCTGCAGATAAACGCGTAAACCCACATcc                                        >  2:195973/5‑248 (MQ=255)
                                                                                                                                                                                                             cAAAAGAAAAACTCAGGGCACGCGAAACGGCGTTCGCTTGAAGTCCGCTGAAAATTATGCCGCAGGCGATGAGCAAAAAGACGGCGAAAAGAACCCCCATCCAGCGCACCCCCAGCCCGCGCGCCTTATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGAc                                 <  2:58691/250‑1 (MQ=255)
                                                                                                                                                                                                              aaaaGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGCGCAAAAGGACGGCGAACAGAACGCCCCTCCAGCGCATCACCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGACACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGTGGTCGCCATGCCGCTAAACGCAGCAACCCACATCCAGAAGACgg                               >  1:320716/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cAACGGCGTTCGCCTGAATTCCGATGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAAATCGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGGGGACGCCATGCAGTTAAACGCGGAAACCCACATCCAGACGACTGCCACAGGTCCACccgcggcaaat        >  2:657629/1‑246 (MQ=255)
                                                                                                                                                                                                                                      aaCGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCACCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAttg       <  1:334784/251‑3 (MQ=255)
                                                                                                                                                                                                                                      aaCGGCGTTCGCCTGAATCCCGCGGAAAATTATGCCGTAGGCGATGACCAAAAAGACGGCGAACAGAACCCCCATCCAGCGCTCCCCCAGCCCGCGCGCCATATCCCAGCCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAttg       <  1:263093/251‑3 (MQ=255)
                                                                                                                                                                                                                                            gatCGCCTGAATTCCGCTGAAAATTATGCCGAAGGCGATGCGCAAAAAGACGGCGAACAGCACGCCCAACCGGCGCATCCCCAGCCCGCGCGCCATATACCATGACGGGCCGCCACGAAACTTCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTAGGCAACAGAAGTAGCCCTGCCGAACAACGCGACAACACACATACAGAAGCAGGCTCCACATCCACcgccgggaaaatacagc  >  2:376219/3‑239 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
CCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGC  >  NC_000913/7232‑7717

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