Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,479 C→T V161I (GTT→ATT)  yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,4790CT96.0% 223.8 / NA 75V161I (GTT→ATT) yaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (37/35);  minor base G (0/3);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAA  >  NC_000913/7233‑7721
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
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                                                                                                                                                      >  2:228132/1‑208 (MQ=255)
                                                                                                                                     aTGAAGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGcc                                                                                                                                                      <  1:228132/208‑1 (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)
                                                                                                                                                               gagcaaaGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAAc                                                                                 >  1:592607/5‑251 (MQ=255)
                                                                                                                                                                agcaaaGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGATCGCGAAACGGAGTTCGCCTGAATTCCGCTGATAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGACCGCGCGCCAAATACCAGGCCGGACCGCCACGAAACTGCCCATTGGCATCACGTTGTTTATAAAGCTGTGCGCGGGAACATTCAGCAAACg                                                                                >  1:277259/4‑251 (MQ=255)
                                                                                                                                                                agcaaaGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACGCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACg                                                                                <  1:105825/248‑1 (MQ=255)
                                                                                                                                                                     aGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTc                                                                           <  2:200336/251‑1 (MQ=255)
                                                                                                                                                                        cAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGcc                                                                        >  1:137501/1‑251 (MQ=255)
                                                                                                                                                                                 aaataccTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATa                                                               <  1:370898/246‑1 (MQ=255)
                                                                                                                                                                                  aataccTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCAGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGCCCGCCACGAAACTGCCCCTTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATa                                                               >  2:555721/5‑250 (MQ=255)
                                                                                                                                                                                      ccTGTCACCAGCGGGGGAAAATCAAAAGAAAAACGCAGGGCGCGCGCAACGGCGTGCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAAcgc                                                          <  1:319605/251‑1 (MQ=255)
                                                                                                                                                                                       cTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCACCGGCGTTAGCCTGAATTCCGCTGCAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGAATCCCCAGCCCGCGCGCCATATACCATGCCGGTACGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAAcgcg                                                         >  2:607514/1‑251 (MQ=255)
                                                                                                                                                                                        tGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGACGCCATGCCGATAAACGCgg                                                        >  1:451862/1‑251 (MQ=255)
                                                                                                                                                                                                      gAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGACCATCCAGCGAATCCCCAGCCCGCGCGCAATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATcc                                            >  1:583644/1‑249 (MQ=255)
                                                                                                                                                                                                       aaaaTAAAAAGAAAAACCCAGGGCTCGCGCAACGGCGTTCACCTGAATTCCGCGGAAAATTAGGCCGTAGGCGATGGGCAAAAAGACGGCGAACAGAACGCCCATCCAGCCCATCCCCAGCCCGCGCGCCTTATACCATGCCGGTCCGCCACGAAACGGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGCCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCaga                                         <  2:519143/251‑3 (MQ=255)
                                                                                                                                                                                                            cAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACg                                    <  1:232223/251‑1 (MQ=255)
                                                                                                                                                                                                             aaaaGTAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACgg                                   <  1:620114/251‑1 (MQ=255)
                                                                                                                                                                                                                 aaaaaaCCCAGGGCTCGGCCAACGGCTTTCGCCTGAATGCCGCTGAAAATGATGCCGGAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCTCGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGTTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAATACGGCTcc                               <  2:288932/250‑1 (MQ=255)
                                                                                                                                                                                                                    aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAgg                            >  1:483979/1‑251 (MQ=255)
                                                                                                                                                                                                                    aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCAAgg                            >  2:260682/1‑251 (MQ=255)
                                                                                                                                                                                                                     aaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGt                           >  2:548345/1‑251 (MQ=255)
                                                                                                                                                                                                                      aCTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTc                          >  2:86900/1‑251 (MQ=255)
                                                                                                                                                                                                                         cAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTcc                         <  1:18816/249‑1 (MQ=255)
                                                                                                                                                                                                                                ggcgcagcGGCGTTCGCCTGAATCCCGCTAAAAATTAGGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCATCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAATCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACcggcg                  <  2:195630/242‑1 (MQ=255)
                                                                                                                                                                                                                                     aaCGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGcgag                                                                                                 >  2:127960/1‑162 (MQ=255)
                                                                                                                                                                                                                                     aaCGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGcgag                                                                                                 <  1:127960/165‑4 (MQ=255)
                                                                                                                                                                                                                                              cgcctGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAATGCCAGCGCac  >  2:279800/5‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
CAAAAGCACTCTCAAAAATAGACCAAATGACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAA  >  NC_000913/7233‑7721

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