Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,518,627 C→A P34H (CCT→CAT)  insA7 → IS1 family repressor protein InsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,518,6270CA93.3% 207.5 / 3.8 75P34H (CCT→CAT) insA7IS1 family repressor protein InsA
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (52/18);  total (55/20)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

GCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGG  >  NC_000913/4518388‑4518872
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
gcgtcttgCATGTGTCGGTCAGAATGACTCATGATGTTGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTTCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTaa                                                                                                                                                                                                                                              <  1:628474/242‑1 (MQ=255)
       gCATGTGCCGGTCAGAATGACTCATGATGTGGACTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGAGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGACGGACATCAGCGCTATCTCTGCTATCATTGCCGTAAAACATGGc                                                                                                                                                                                                                                      >  2:76367/1‑250 (MQ=255)
             gtCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATGTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGAAAAAGCACTGCCGGACATCAGCGCTATCGCTGCTCTCATTGCCGTAAAACATGGCAACTACa                                                                                                                                                                                                                               <  1:174681/251‑1 (MQ=255)
                ggTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGtt                                                                                                                                                                                                                            >  1:355505/1‑251 (MQ=255)
                   cAGAATGACTCATGAAGTGGTCTGCTATTATTGACATCCTCACGGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCa                                                                                                                                                                                                                          <  1:162536/250‑1 (MQ=255)
                      aaTGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACtt                                                                                                                                                                                                                       >  1:244330/1‑250 (MQ=255)
                      aaTGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACtt                                                                                                                                                                                                                       >  2:209079/1‑250 (MQ=255)
                      aaTGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTa                                                                                                                                                                                                                      <  1:210797/251‑1 (MQ=255)
                       aTGACTCATGATGTGGTCTGCTATTATTGACATCCTCACGGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTac                                                                                                                                                                                                                     <  2:594671/251‑1 (MQ=255)
                        tGACTCATGATGAGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGAGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGTAGGCGTGGTGCGTAACGGCAAAAGAACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCACCTACAGTTCACTTac                                                                                                                                                                                                                     >  2:499060/1‑250 (MQ=255)
                            tCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGc                                                                                                                                                                                                                >  1:416131/1‑251 (MQ=255)
                                gatgTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCctct                                                                                                                                                                                                            >  1:124853/1‑251 (MQ=255)
                                gatgTGGACTGCTATTATTGACATCCTCACTGCCCTCAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCATTCCTGCTCCGCTACTGAAGTCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTCCAGTTCACTTACACCGCctct                                                                                                                                                                                                            >  2:598222/1‑251 (MQ=255)
                                     ggTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGACACTCCAGTTCACTTACACCCCCTCTCAGcc                                                                                                                                                                                                       >  2:432513/1‑251 (MQ=255)
                                      gttttcTATTTTTAACATCCTCTCGGCCCTAAGGGATGGGGATTTCGGTAATGCTCCCAACTTATTGTTTTAGTGTATGATGGTGATTTTAAGGTGCTGGCGTGGTTTCCATTCCCACCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGAAAAAGCACTGCCGGACATCAGCTCTACCTCTGCTCTCATGGCCTTAAAACATGGAAACTACATTTCCCTTACACCGCCTCTCAGccc                                                                                                                                                                                                      <  2:375032/246‑1 (MQ=255)
                                      gTCTGCAATTCTTGCCATCCTAACTGCCCTAACGGACGGGGATTTCGGTACTGCTGTCAACTTACTGATTTAGAGTATGATGGTGATTTTACGGTGCTTGCGTGGTTTACATTTCCCTCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACAGCAATCGCCCTGCCGGCCATCAGCGCTATCTCTGCTCTCATTGCAGTAAAACATGGCAACTACAGTTCACTTCCACCGCCCCTAAGccc                                                                                                                                                                                                      >  2:324735/1‑251 (MQ=255)
                                          gCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTAAGTGTATAATGGAGATTTTAAGGTGCTTGCGGGGATTCCATTTCAATCAGATGTCCTTCCTGCACCGCTACTGAAGGCGGGGTGCGTAACGGCAAAAGCACTGCCGGACATACGCGCTATCTATCCTCTCATTGCCAGAAAACAGGGCAGCTACAGGTCACTTACACGGCCTCTCAGCCCGGt                                                                                                                                                                                                   >  2:159394/1‑250 (MQ=25)
                                              ttattGACATCCTCACTGCCCTAACGGCTGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGCTTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTacgc                                                                                                                                                                                               >  1:145022/1‑247 (MQ=255)
                                                  tGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCAcc                                                                                                                                                                                            >  1:384643/1‑249 (MQ=255)
                                                       tCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAAt                                                                                                                                                                                     >  1:250939/1‑251 (MQ=255)
                                                           cccTGCCCTAAAGGAGGGGGATTTCGGAAATGCTGCCAACTGACTGATTTAGTGTATGATGGTGATTTAAAGGGTCGGGCGTGGCTTCCATTTCCATCAGATGCCCTTCCGGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACAACAGTTCGCTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCAt                                                                                                                                                                                  <  2:338291/248‑1 (MQ=21)
                                                           caccgccttaagggATGGGGTTTTCGTTAATGCGGACATTTTCCTTATTTAGTGGATGTTGGTGATTTTAAGGTGCTGGTGGGGCTTCCTTTTCCATCAGATTTCTTTCCTGCCTCGCTACTGATGGCGTGGTGCGTAACGGCAAAAGCACGGCCGGACACCAGCGCTATCTCTGCTCCCCTTGCAGTAAAACATGCCAACTACATTTTCCTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCAtt                                                                                                                                                                                 <  2:124853/239‑1 (MQ=255)
                                                             cTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTg                                                                                                                                                                                >  1:148038/1‑250 (MQ=255)
                                                                cccTAAAGGGTGGGGATTTCGGTAAAGCTGACAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGACAATCATTGatat                                                                                                                                                                            >  1:453415/1‑251 (MQ=255)
                                                                   tAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGc                                                                                                                                                                         >  1:104046/1‑251 (MQ=255)
                                                                     aaGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTa                                                                                                                                                                                                                                               <  2:544725/179‑1 (MQ=255)
                                                                     aaGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTa                                                                                                                                                                                                                                               >  1:544725/1‑179 (MQ=255)
                                                                      aGGATGGGGATTTCGGTAATGCTGCCAAATTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCAt                                                                                                                                                                      >  1:578288/1‑251 (MQ=255)
                                                                                    ggTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGAt                                                                                                                                                        >  1:578228/1‑251 (MQ=255)
                                                                                      tAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGt                                                                                                                                                      >  1:67154/1‑251 (MQ=255)
                                                                                          gctgcCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTTTCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTCCACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTcgcg                                                                                                                                                  <  1:355246/251‑1 (MQ=255)
                                                                                           ctgcCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCgcgc                                                                                                                                                 >  1:75880/1‑251 (MQ=255)
                                                                                              ccAACTTACTGCTTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGATTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTACGCCCGGTACGCACCAGAACATCATTGATATGGACATGAATGGCCTCGGATGTCGAGcc                                                                                                                                                >  2:330850/1‑249 (MQ=18)
                                                                                                   ttACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCAc                                                                                                                                         >  1:123662/1‑251 (MQ=255)
                                                                                                        gATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCAtt                                                                                                                                    >  1:261068/1‑251 (MQ=255)
                                                                                                          tttatgtttggttggggATTTTAGGGTGCTTGTGGGGTTTCAATTCCAATCAGATGTCCTCCCTGCCCCCCTACGGAAGGCGTGGGTCGTAACGCCAAAAGCACTGCCGGACAGCAGCGCTACCTCTGCTCTCATGGCCGAAAAACATGGCAACTACAGTTCATTTCCACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTa                                                                                                                                   <  2:578228/234‑1 (MQ=12)
                                                                                                             aGTGTATGATGGTGCTTTTAAGGAGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTTTggg                                                                                                                               >  2:9514/1‑251 (MQ=255)
                                                                                                             aGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTCCTGAAGGCGTGGTGCGTAACGGCAAAACCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCCCCAGAAAATCATTGATATGGCAATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATggg                                                                                                                               >  1:436884/1‑251 (MQ=255)
                                                                                                                gtATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGt                                                                                                                            >  1:169527/1‑251 (MQ=255)
                                                                                                                 tATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTAGGATGTCGCGCCAGTGCACGCATTATGGGCGtt                                                                                                                           >  1:183760/1‑251 (MQ=255)
                                                                                                                  atgatgGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTg                                                                                                                          <  1:209079/251‑1 (MQ=255)
                                                                                                                    tttgGTTATTTTAGGGTGCTTGCTTGTCTTCCATTTCCATCCGTTTTCCTTCCTGCTCCGCTCCTGAAGGCGTGGTGCTTACCGCCAAAAGCCCTGCCGGACATCAGCGCTATCTCTCCTCTCATTGCCGTAAAACTTGGCAACTCCATTTCCCTTACTCCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTgg                                                                                                                         <  1:159394/248‑1 (MQ=18)
                                                                                                                       tgtgtttttAAGTTGCTGCCGTGGCTTCCATTTCCACCAGATGTCCTCCCTGCTCCGCGACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCACCTCTATCTCTTCTCTCATCGCCGTAAAACATGGCAACTCCAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCt                                                                                                                      <  1:494643/245‑1 (MQ=18)
                                                                                                                        gTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCCCCGGTGCACGCCTTATGGGCGTTGGCCTc                                                                                                                     >  2:524962/1‑250 (MQ=255)
                                                                                                                                    tGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTAc                                                                                                        >  1:395023/1‑251 (MQ=255)
                                                                                                                                        tGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGt                                                                                                                                             >  1:52072/1‑210 (MQ=255)
                                                                                                                                         tcGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCCCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCa                                                                                                    <  1:76367/249‑1 (MQ=255)
                                                                                                                                          cGTGGCTTCCATTTGCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAACGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCCCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGACAGGCGTCGGATGTCGCTCCAGTGCACGCATTTTGGGCGTTGGCCTCAACACGGTTTTACGTCACt                                                                                                  >  2:608144/1‑251 (MQ=21)
                                                                                                                                          cGTGGCTTCCATTTACATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACt                                                                                                  >  1:454920/1‑251 (MQ=255)
                                                                                                                                          cGTGGCATCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGGAACGGCAAAAGCACTGCCGGCCATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGATCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCAATGACTGGCGTCGGCTGTCGCCCCAGTGCCCGCATTATGGGCGTTGGCCTCCACCAGGTTTAACGTCACt                                                                                                  >  2:361168/1‑251 (MQ=18)
                                                                                                                                            tGGCTTCCATTACCATCAGATATCCTTCATGCTCTGCTACTGAAGGCGTGGTGCTAAACGGCAAAAGCACTGCCGGCCATAAGCGCTATCTATGCTCTCATTGACGCAAACCATGGCAACTACAATTCACTTACACCGCCTCTCAACCCGGCACGCACCAGAAAATCATTGATATGGCCATGAATGACCTCGGATGTCGCGCCAGTGCACGCATTCTGGGACTTGGCCTCACCACGAGTTTACGTCACtt                                                                                                 >  2:410558/1‑250 (MQ=12)
                                                                                                                                                    atttccatCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTc                                                                                        >  1:393899/1‑251 (MQ=255)
                                                                                                                                                            cAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCATCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGATTTACGTCACTTAAAAAACTAAGGcaga                                                                                 >  2:234401/1‑247 (MQ=255)
                                                                                                                                                            aagagtctccgcctGCTCCGCTACTAAAGGAGTGGGTCGTAACGGAAAAAGCACTCCCGGACACCATCGCTATTCCTGTCCTCATTCCCGTAAAACATGACATCTAAATGGCCCTTACCCCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAATCTCAGTCCgta                                                                                <  2:649806/240‑3 (MQ=12)
                                                                                                                                                               aTGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTc                                                                             >  2:257713/1‑251 (MQ=255)
                                                                                                                                                                      tcctGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCACAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAACTCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTATCGTCCCTTAAAAAACTCAGGACGCAGTCGGTAAcc                                                                      >  2:22265/1‑251 (MQ=18)
                                                                                                                                                                             ccgtaacgggtgggggggggCGTAACGGCAAAAGCTCCGCCGTACATCAGCGCTATCGCTGTTCTCACTGCTGTAAAACATGGCAACCGCAGTTCACTTACACCGGGTTTCAACCCGGGACGCACCAGAAAATCATTGATAGGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGTCTCAACACGGTTTGACGGCACTTAAAAAACGCAGGCCGTTGGCGGTATCCTcgcga                                                                <  2:402373/231‑2 (MQ=255)
                                                                                                                                                                                cTACTGCAGGCGTGGGGCGCAACGACAAAAGCAATGCCGGACATCAGCGCTATCTCTGCGCTCATTGCCGTAAAACATGGCAACTAAAGTTCACTTACACCGCCTCTCAGCCAGGTACGCACCCGCAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGAACGCATTATGGGCGTTGACCTCAACACGGTTTAACGGCACTTCAAAAACTCCGACCGAAGTAGGCAACATCGCGCATAc                                                            >  2:318694/1‑251 (MQ=255)
                                                                                                                                                                                    tGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTCCGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGAATACAAcc                                                        >  2:416177/1‑251 (MQ=255)
                                                                                                                                                                                               tgCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGACTGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAgtgatg                                             >  1:53771/1‑251 (MQ=255)
                                                                                                                                                                                                   tAACGGCAAAAGCACTGCAGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGAGTCGGATGTCGCGCCAGTGAACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCCATTACAAAACTCAGGCCGAATGCGGTAACCACACGCCTCCAACCGGGCAgagatctgc                                          >  2:527949/1‑246 (MQ=17)
                                                                                                                                                                                                                 cTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGaa                            >  1:20402/1‑250 (MQ=255)
                                                                                                                                                                                                                   gCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATg                         >  1:289915/1‑251 (MQ=255)
                                                                                                                                                                                                                       gACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGTTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACg                     >  1:361940/1‑251 (MQ=255)
                                                                                                                                                                                                                        aCATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTc                                                                                                                     >  1:330183/1‑154 (MQ=255)
                                                                                                                                                                                                                        aCATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTaa                                                                                               >  1:673272/1‑176 (MQ=255)
                                                                                                                                                                                                                        aCATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGGTTTACGTCCCTTAAAAAACTCAGGCCGTGGTCGATGACCTCGCGCATACAACCGGGCAGAGATGGGATTGACTGCCATGCAATGGACg                     >  2:494490/1‑250 (MQ=17)
                                                                                                                                                                                                                            cAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGAGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGCAATGGACGAAc                  >  2:13274/1‑249 (MQ=255)
                                                                                                                                                                                                                              gcgcTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTg              >  1:507770/1‑251 (MQ=255)
                                                                                                                                                                                                                               cgcttgctctGCCCTCATTGCCGTAAAACAGGGCAACTACAGTTACCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTgg             <  1:324983/245‑1 (MQ=18)
                                                                                                                                                                                                                                     ccctcctctCATTGCCTTAAAACATGCCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTAc       <  2:261068/246‑1 (MQ=255)
                                                                                                                                                                                                                                       ctGCTCTCATTGCCGTAAAACATGGCCAATACAGTTAACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATAGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAACACCTCAGGCCTTCGTCGGAAACCTCGCGCATACAGCCGGACGGTGATGTGATTGTCTGCGATGAAACGGCCGACCGGTGGGGCTACgg     >  2:219263/1‑250 (MQ=17)
                                                                                                                                                                                                                                        tGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGAAGTCGGTAACCTCCCGCATACAACCGGGCAGTGTTGTGCTTGTCTGCGCTGAAATGGAAGAACGATGGGGCTActcc    >  2:329160/1‑248 (MQ=18)
                                                                                                                                                                                                                                         gCTCTCATCGCCGAAAAACTGGGCAACTACAGTTCCATTCCACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGACATGGCCATGAATGGCGTCGGATTTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACgcc    <  1:417646/250‑3 (MQ=255)
                                                                                                                                                                                                                                          ctctCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCgg  >  1:596941/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
GCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGG  >  NC_000913/4518388‑4518872

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