Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,395 T→C T189A (ACT→GCT)  yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,3950TC98.6% 226.0 / ‑4.8 70T189A (ACT→GCT) yaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (35/33);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

CATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAA  >  NC_000913/7156‑7641
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
cATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCga                                                                                                                                                                                                                                             <  1:221938/251‑1 (MQ=255)
   aCTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCgaga                                                                                                                                                                                                                                           <  1:541696/250‑1 (MQ=255)
                    tAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGc                                                                                                                                                                                                                         >  1:429498/1‑251 (MQ=255)
                     aGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCg                                                                                                                                                                                                                        >  1:342256/1‑251 (MQ=255)
                             ccTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGaaaa                                                                                                                                                                                                                <  1:334190/251‑1 (MQ=255)
                                   ttAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGCAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGCGTGATCGCCAGCAGAGCAAAGACAGCGCGAATACTAACTGTCACCCGCGGGGGAAAATCaaa                                                                                                                                                                                                           >  2:307333/1‑250 (MQ=255)
                                     aaGGTATATCCCGCCGCGCCGCCGGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTTCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGa                                                                                                                                                                                                        <  2:472022/251‑1 (MQ=255)
                                        gTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGaaa                                                                                                                                                                                                      >  1:548305/1‑250 (MQ=255)
                                          atatCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGACGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTCCCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGCCGCCATGCAGACCGCGAGTGATCGACAGCAGAGCAAAGACAGCGCGATTAATACCTGTCACCAGCGGGGGAAAATCAAACGAAAAAc                                                                                                                                                                                                   >  1:627593/1‑251 (MQ=255)
                                            atCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTc                                                                                                                                                                                                 >  1:159119/1‑251 (MQ=255)
                                              cccGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTc                                                                                                                                                                                                 >  2:370898/1‑249 (MQ=255)
                                              cccGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCa                                                                                                                                                                                                >  2:495155/1‑250 (MQ=255)
                                               ccgccgCGCCGCCTTACGATTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCCAATGACGTGGGGCAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGGGCATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAgg                                                                                                                                                                                              >  2:68935/1‑251 (MQ=255)
                                                     cgccgcCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGgctcgc                                                                                                                                                                                        >  1:297141/1‑247 (MQ=255)
                                                        cgccTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGgctcgcgca                                                                                                                                                                                     <  2:157195/251‑8 (MQ=255)
                                                         gccTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCATAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGTGGGAAAATCAAAAGAAAAACTCAGGGTGCGGGCaa                                                                                                                                                                                    <  2:318265/251‑1 (MQ=255)
                                                         gccTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAAGAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAATCCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGAGCGCCAGCAGAGCAAAGACAGCTAGAATAAGAGCTGTCACCATCGGGGGAAAATCTAAAGAAAAACTCAGGGCTCGCGCaa                                                                                                                                                                                    <  2:591513/251‑1 (MQ=255)
                                                            tGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACgg                                                                                                                                                                                 >  1:579691/1‑251 (MQ=255)
                                                             gCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACg                                                                                                                                                                                  <  1:168602/249‑1 (MQ=255)
                                                                             caAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCCCCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTGCGCCTGAATTCCg                                                                                                                                                                <  2:114360/251‑1 (MQ=255)
                                                                                       ctcAAAAAGAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAAt                                                                                                                                                        <  2:429498/249‑1 (MQ=255)
                                                                                         cAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTa                                                                                                                                                      <  2:43628/249‑1 (MQ=255)
                                                                                           aaaaTAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGcc                                                                                                                                                  >  1:173366/1‑251 (MQ=255)
                                                                                            aaaTAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGccg                                                                                                                                                 >  2:383276/1‑250 (MQ=255)
                                                                                                         gACGTGGGGAAGTTGCCCGATATTCATTACGCAAAACACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTGCGCCTGAAGTCCGCTGAAAATTATGCCGTAGGCGATGAGCa                                                                                                                                    <  2:83348/251‑1 (MQ=255)
                                                                                                            gTGGGGAAGTTGCCCTATATTCATTACGCAAATCAACAGGCGGGACATGACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCAGGAAGACCGCGAGTGAGCGCCAGCAGAGAAAAGACAGCGAGAATAATTCCTGTCACCTGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATGCCGCTGTAAATTATGCCGTAGGTGCTGAGCaaaa                                                                                                                                 <  2:77371/251‑1 (MQ=255)
                                                                                                                   agtagCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGc                                                                                                                          >  2:53236/5‑251 (MQ=255)
                                                                                                                   aGTTGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGc                                                                                                                          >  1:467639/1‑251 (MQ=255)
                                                                                                                   aGTTGCCCGATATTCATTACGCAAAGCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGCGCAACGGCGTGCGCCTGAAGTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGc                                                                                                                          <  2:579691/251‑1 (MQ=255)
                                                                                                                      tGCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGa                                                                                                                        >  1:116635/1‑250 (MQ=255)
                                                                                                                       gCCCGATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAAGTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGaa                                                                                                                       <  2:467639/250‑1 (MQ=255)
                                                                                                                           gATATTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAACGACGGCGAACAGa                                                                                                                   >  2:324198/1‑250 (MQ=255)
                                                                                                                              atTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGa                                                                                                                                       >  1:102560/1‑227 (MQ=255)
                                                                                                                              atTCATTACGCAAATCACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGa                                                                                                                                       <  2:102560/227‑1 (MQ=255)
                                                                                                                                               ccAGGCTGGACAGTACCCAGCTTATCGCCATCAACGGGACAAAACCCTGCATGAGACGGGCGACGCCCTGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATAACTGTCACCAGCGGGGGAAAATAAAAAGAAAACCTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCCAAAAGCCGGCGAACAGAACGCCCATCCAGCGCATcccc                                                                                              >  2:340817/1‑251 (MQ=255)
                                                                                                                                                    cTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGccc                                                                                         >  1:596751/1‑251 (MQ=255)
                                                                                                                                                       gTCAGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCcgcg                                                                                      <  1:544040/251‑1 (MQ=255)
                                                                                                                                                          aGTACCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGCAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCgcgcgc                                                                                   >  2:440603/1‑251 (MQ=255)
                                                                                                                                                             aCCCAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATAACTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCat                                                                                >  2:18816/1‑251 (MQ=255)
                                                                                                                                                                cAGATTATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCatat                                                                              >  2:68324/1‑250 (MQ=255)
                                                                                                                                                                    ttATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCAt                                                                         >  1:135027/1‑251 (MQ=255)
                                                                                                                                                                    ttATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCAAGCGCATCCCCAGCCCGCGCGCCATATACCAt                                                                         >  2:531352/1‑251 (MQ=255)
                                                                                                                                                                     tATCGCCATCAACGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCAt                                                                         <  2:135027/250‑1 (MQ=255)
                                                                                                                                                                        cGCCATCAACGGGACAAACCCCTTCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGcc                                                                      <  1:431859/250‑1 (MQ=255)
                                                                                                                                                                           cATCAACGGGACAAACCCCGGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCAGCCCCAGCCCGCGCGCCATATACCATGCCGGTc                                                                  <  1:307333/251‑1 (MQ=255)
                                                                                                                                                                            aTCAACGGGACCAACCCCTGCATGCGCAGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGCGCAAAGACAGCGAGAATAATACATGTCACCAGCGGGGGCAACTCAAAAGACAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAACAAGACGGCGAACAGAACGCCCATCCAGCGCCTCCCCAGCCCGCGCGCCATATACCATGCCGGTcc                                                                 >  2:623283/1‑251 (MQ=255)
                                                                                                                                                                                 cGGGACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCAc                                                            <  1:144251/251‑1 (MQ=255)
                                                                                                                                                                                  gggACAAACCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACg                                                           >  1:273379/1‑251 (MQ=255)
                                                                                                                                                                                         accccTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGc                                                    <  2:182476/249‑1 (MQ=255)
                                                                                                                                                                                            ccTGCATGAGCCGGGAGACGCCAAGAAGACCGCGAGTGATCGCCCGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGCATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGCCGGCGACCAGAACGCCCATCCAGCGCATCCCCAGCCCGGGCGCCATATACCATGCCGGTCCGCCCCGAAACTGcccc                                                 >  2:362857/1‑250 (MQ=255)
                                                                                                                                                                                               gCATGAGCCGGGCGACGCCAGGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTGCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCAtt                                               <  1:444803/250‑1 (MQ=255)
                                                                                                                                                                                                cATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTgg                                             <  1:312202/251‑2 (MQ=255)
                                                                                                                                                                                                 aTGAGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTggc                                            <  1:73588/251‑3 (MQ=255)
                                                                                                                                                                                                   gAGCCGGGCGACGCCAAGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATTCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGGTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTggcat                                          <  2:169411/251‑5 (MQ=255)
                                                                                                                                                                                                    aGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCCAAAGAACAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGcc                                                   >  2:302972/1‑241 (MQ=255)
                                                                                                                                                                                                    aGCCGGGCGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGcc                                                   <  1:302972/241‑1 (MQ=255)
                                                                                                                                                                                                     gCCGGGCGACGCCATGAAGACCGCTAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCGGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTggcatca                                        <  1:410081/251‑7 (MQ=255)
                                                                                                                                                                                                      ccGGGCGACGCAATGAAGACCGCGAGTGATCGCCAGCCGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTTAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCCGCGCATCCCCAGCCAGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTggcatcac                                       >  1:278318/1‑244 (MQ=255)
                                                                                                                                                                                                       cGGGCGACGCCATGAAGCCCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGTGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACg                                      <  1:463049/251‑1 (MQ=255)
                                                                                                                                                                                                           cGACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCt                                  <  1:419416/251‑1 (MQ=255)
                                                                                                                                                                                                            gACGCCATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCtt                                 >  2:182650/1‑251 (MQ=255)
                                                                                                                                                                                                                 cATGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGCAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCCTCCCGTTCTTtata                             >  2:538903/1‑250 (MQ=255)
                                                                                                                                                                                                                  aTGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGcc                                                                      <  1:228132/208‑1 (MQ=255)
                                                                                                                                                                                                                  aTGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGcc                                                                      >  2:228132/1‑208 (MQ=255)
                                                                                                                                                                                                                       gACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCtgt                      <  1:169284/251‑1 (MQ=255)
                                                                                                                                                                                                                          cgcgAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCTGTCCGCCACGAAACTGCCCATTGGCATCACGtt                                    <  1:142612/234‑1 (MQ=255)
                                                                                                                                                                                                                          cgcgAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGtt                                    >  2:142612/1‑234 (MQ=255)
                                                                                                                                                                                                                             gAGTGATCGCCATCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCAGCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTGCTTTATAACGCTGTGcgagg                <  2:56730/251‑5 (MQ=255)
                                                                                                                                                                                                                                  atcgCCAGCAGAGCCAAGACAGCGAGAATAATACCTGTCAACAGCGGGGGAAAATCAAAAGAAAAACTCAGAGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGCTGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCCTCCCCAGCCCGCGCGCCATATACCATGCAGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTCTAAAGCTGTGCGAGGGAAc            >  2:249570/4‑250 (MQ=255)
                                                                                                                                                                                                                                           agagCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCCAACAGAACGCCCATCCAGCGCATCCCCATCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGTCATCACGTCCTTTATAAAGCTGTGCTAGGGATCATTCGGCaaa  <  1:484529/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
CATACTGCGCTGAAAACCGTTAGTAATCGCCTGGCTTAAGGTATATCCCGCCGCGCCGCCTGCCGCTTCCTGCCAGCCAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAA  >  NC_000913/7156‑7641

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