Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,619,269 A→G V26V (GTA→GTG rhsB → Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,619,2690AG97.0% 215.7 / NA 68V26V (GTA→GTGrhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (40/25);  minor base T (2/0);  total (43/25)
Fisher's exact test for biased strand distribution p-value = 5.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGG  >  NC_000913/3618999‑3619461
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                
tGGTTCCTTTGGTGAACGATCAACAAAGGGCCACTAGATCACCTGATTACTCTTTAAATACCCTATCCGAACCTTATCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTg                                                                                                                                                                                           <  4:634133/278‑1 (MQ=255)
tGGGTCCGTGTGTGAAACATAAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATCCCCTATCCGAAACTTCGCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGAATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTg                                                                                                                                                                                           <  4:410203/278‑1 (MQ=38)
   gTCCTTTGGTGAACGATCAACAAAGTGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAACCTTCACTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTg                                                                                                                                                                                           <  1:683600/275‑1 (MQ=255)
      cTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTAACCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTTTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGTGTCAGGGTGACATGACGCAGTATGGCGGAAGCATTGTTCAGGGTTCAGTCGGGGTGCGCATTGGTg                                                                                                                                                                                        >  1:459457/1‑275 (MQ=255)
        ttGGGGAACGAACAACAAAGGGCCAATAGCACAACTGATTCCTCTCTAAATACCCATTCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGcc                                                                                                                                                                                      <  3:332151/275‑1 (MQ=255)
            tGAACGATCAACAAAGGGCCAATAGCACAGCAGATTCCTCTCTCAATACCCTATCCGAACCGTCTCTTTTGGAACGTTCTAAAGATAATCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCCGAAAATCGGCAGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCa                                                                                                                                                                                  <  3:459456/275‑1 (MQ=38)
              aaCGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTTCGTATTGGGGCCTCCACCGGc                                                                                                                                                                             >  2:582155/1‑278 (MQ=255)
              aaCGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGAATTTTACGATGAGCGCAAAACCGGCGGCGCGCCAGGGTGACATGACGCATTTTGGCGTTAGACTTGTTAAGGTTTCACCCGGGTTGCGCATTGGTtcctcc                                                                                                                                                                                   >  1:426949/1‑272 (MQ=255)
                  ataaacaaAGGGCCACTTGAACACCTGATTACTCGCTAACTACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAAATCAATGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATGTTACGATGGGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGc                                                                                                                                                                             <  3:678741/271‑1 (MQ=37)
                  aTCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGc                                                                                                                                                                             <  3:546681/274‑1 (MQ=255)
                  aTCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGc                                                                                                                                                                             <  1:588548/274‑1 (MQ=255)
                                 cTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTTAGGGTTCCGCCGGGGTGCGCATTGGTGCCCCCACCGCCGTGGCCTGTTGTGTGTGcg                                                                                                                                                          >  4:224187/1‑277 (MQ=255)
                                 cTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCAGTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGCTCGGTGTGcc                                                                                                                                                          >  2:218604/1‑278 (MQ=255)
                                    gCACACCTGATTCCTTTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAACTCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCACTCGGGGTGCGCACTGGTGCCCCCCCCGGCGTGGCCTGTTCGGTGTGc                                                                                                                                                           >  1:235818/1‑274 (MQ=255)
                                      acacCTGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATGGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCGCCACCGCCGTGGCTGGTTCGGTGTGcc                                                                                                                                                          >  4:429570/1‑273 (MQ=255)
                                       cacCTGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcggc                                                                                                                                                     >  4:434333/1‑277 (MQ=255)
                                          cTGATTCCTCTCTAAATCCACTATCCGCACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGGATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGTAAAACCTGCGGCGCGTCAGGGTTACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAg                                                                                                                                                 <  2:137526/278‑1 (MQ=255)
                                           tGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGt                                                                                                                                                <  2:274931/278‑1 (MQ=255)
                                            gATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCCGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGTGTGCGCATTGGTGCCCCCCCCGGCGTGTCCTGGTCGGTGTGCCTCGGCGGAGt                                                                                                                                                >  4:502041/1‑277 (MQ=255)
                                               tCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGATCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCACCCACCGGCGTGGCCAGTTCGGTGTCCCCCTGCGGAGTGCCg                                                                                                                                            >  2:418178/1‑278 (MQ=255)
                                                  tctctAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCTTTGTTCAGGGTTCAGCCGGGGTGCGCATTGTTTCCCCCACCGGCTGGGCCTGTTCGGTGGCCCCcggcgggtggac                                                                                                                                             >  3:703514/1‑268 (MQ=255)
                                                  tctctAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGcc                                                                                                                                                                                      >  1:191563/1‑233 (MQ=255)
                                                    tctAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTc                                                                                                                                          >  1:165867/1‑275 (MQ=255)
                                                    tctAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTc                                                                                                                                          >  3:336989/1‑275 (MQ=255)
                                                      tAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCa                                                                                                                                                                                  <  3:191562/233‑1 (MQ=255)
                                                          tACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATGTATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCa                                                                                                                                    <  1:31936/275‑1 (MQ=255)
                                                          atcccTATCCGAACCTTCACGTTTTTAACGTTCTAGATATTTTCATAAACATCATCAATTCATTGTGACCACAAGTTTTTCTTGGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTATACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTANCATGGTTCAGGGTTCAACCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATcc                                                                                                                                 <  4:352936/276‑1 (MQ=21)
                                                               tctccgaggcatctcTATTGTAACGTTCTCAATCTAGTCCTGAAAAACTTCAATTCATTGTGACCACAAGTTTTTCTACGATTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGTTTTTACGATGAGCGGAAAACCGGAGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTNTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCTTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCa                                                                                                                            <  4:354814/267‑1 (MQ=21)
                                                                acccGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACCAGTTTTTCTTCGCTTTTTCGTCTGAAGATACTGTCATCAAAATAATAGAAAAGGATTTTACGATGAGCGGACAACCGGCGGCGCGTCCGGGTGACATGACGCAGTATGGCGGTAGCAGTGTTCAGGGTTCACGCGGGGTGCGAATTGGTGCCCGCACGGGGGTGGCCTGTTCGGGGTGACCCGGGGGTTTGTCGTCCGGCCCTcccgg                                                                                                                              >  3:453679/3‑272 (MQ=255)
                                                                aTCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGt                                                                                                                              >  1:727709/1‑275 (MQ=255)
                                                                aTCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGCATGGCGGTAGCATTGTTCAGGGTTCACGCGGGGGGCGCATTGGTTGCCCCAGCGGCGTGGCCGGTCCGGTGCTCCGCGGCGGGGTGGCGGCCGGCCATCAGGt                                                                                                                              >  3:111629/1‑275 (MQ=255)
                                                                       ccTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATccc                                                                                                                       >  1:505307/1‑275 (MQ=255)
                                                                           ctctTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATCCTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTTCGCATTGGTGCCCCCCCCGGCGTGGCCTGTTCGGTTTGCCCCGGGGGCGTGACGTCCGCCCCTCCGGTCACTCCCTTg                                                                                                                    >  3:331569/1‑274 (MQ=25)
                                                                            tctTTTGTAACGTTCTAAATATATTCCTTAACATATTCAATTCCTTTTGACACCAAGTGTTTCTTCGTTTTTTCGTATGAAGATACTTTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGTTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAACCCCCTGCt                                                                                                                  <  3:202263/275‑1 (MQ=21)
                                                                             ctTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTc                                                                                                                 <  1:336990/275‑1 (MQ=255)
                                                                              ttgtGTAACGTTCTAAATATTTTGCTAAAAATCTTCAATTCATGTGGACCTCAAGTTTTTCTTCGCTTGTTCGTATGAATATACTGTCATTAAAAGAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGGCAGGGTGACAGGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCg                                                                                                                <  1:111630/272‑1 (MQ=18)
                                                                                 tGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACCTCGGGCCACACGGTCAAACCCCTGCTCGGTg                                                                                                             >  1:528976/1‑275 (MQ=255)
                                                                                        ttCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAgt                                                                                                       >  3:368164/1‑273 (MQ=255)
                                                                                        ttCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGGGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAgg                                                                                                       >  1:584047/1‑274 (MQ=255)
                                                                                         tCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGGGCGGAAAGCCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATcc                                                                                                                                 <  2:111255/247‑1 (MQ=255)
                                                                                                  atTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGa                                                                                          >  2:408331/1‑277 (MQ=255)
                                                                                                  atTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGGCCTTCCCGGTGa                                                                                          >  2:631802/1‑277 (MQ=255)
                                                                                                  atTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCCCGGTGCAAAGGTCCTTCCCgg                                                                                             >  3:310209/1‑274 (MQ=255)
                                                                                                     ccTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGACAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATTACGCAGTATGGCGGTAGCATCGTTCAGGGTTCAGCCGGGGGCCGCACTGGTGCCCCTACCGGCGTGTCCTGGTGGGTGTTCCCCGCCGGTGTGCCGTCCGGCCATTCGGTCACCCCCCTGCCGGGTGCAACGGTCCTTCCCGGt                                                                                            >  3:56910/1‑272 (MQ=255)
                                                                                                         aaaaTCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCCCGGGGCAAAGGGCCTTCCCGGTGAAACCGAc                                                                                   >  4:423452/1‑277 (MQ=255)
                                                                                                                tCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAACCCGACATCGccc                                                                            >  2:177619/1‑277 (MQ=255)
                                                                                                                       aTTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCGGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTTTTCCCCGGCGGAGTGACGTCCGGCCACCCGGTCCATCCCCCGCCCGGTGCAAAGGTCCTTCCCGGCGAACCCGACATCGCCCTgcc                                                                        >  1:460762/1‑274 (MQ=18)
                                                                                                                            gACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAACCCGACATCGCCCTtccctgcc                                                                   >  3:57832/1‑274 (MQ=255)
                                                                                                                            gACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCATGGTGAAACCGACAACGCCCTgcccggc                                                                    >  1:118997/1‑273 (MQ=255)
                                                                                                                             aCCACAAGTTTTTGTTCTCTGTTTCGTCAGAGATGAATGGCATTAAAAAAATAGAAAAGGATTGCACGATGGGCGGAAAACCGGCGGCGCGTCAGGGTGACAGTACGCAGTGTGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATAGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTg                                                              <  4:408331/278‑1 (MQ=14)
                                                                                                                                 caAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTTCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCg                                                           >  2:194235/1‑277 (MQ=255)
                                                                                                                                  aaGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCATTCCCCGGCCCGGTTCAAAGGGCCTTCCCGGTCAAACCGACATCGCCCTGCCCGGCCCGCTGc                                                             >  3:435656/1‑274 (MQ=255)
                                                                                                                                       tttCTTCGCTTTTTCGTATGAAGATACTGTCATGAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATc                                                     <  2:429570/277‑1 (MQ=255)
                                                                                                                                            tcGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCAGTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCCTCCGGTCAATCTCCTGCTCGGTGCCAAGGTCCTTGCCGGTGAAACCGACATCGCCCTTCCCGGTCCCCTGCCGTTctcccc                                                   >  3:235411/1‑269 (MQ=17)
                                                                                                                                                 ttttCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGc                                            <  2:434333/276‑1 (MQ=255)
                                                                                                                                                 ttttCGTATGAAGAGACGGTCATTAAAATAATAGAAAAGGATTTTACGAGGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGc                                            <  2:502041/276‑1 (MQ=255)
                                                                                                                                                    tCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGc                                            <  3:727708/273‑1 (MQ=255)
                                                                                                                                                    tCGTAAGAAGATAATGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGGGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGc                                            <  1:301235/273‑1 (MQ=255)
                                                                                                                                                                  cTGTCATTAACATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGNTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGCCGTCCGGCCACCCGGTCAATCCACTGCTCGGTGCAAAGGGCCTTCCCGGTGAAACCGTCATCGCCCGGCCCGGCCCTCTGCCAGTCACCCTCGCCCGCACCTAccgcatgcaccg                           >  4:331233/1‑270 (MQ=14)
                                                                                                                                                                      cATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCTCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCACAGGTCCTTCCCGGTGAAACCGACATCGCCCTCCCCGGCCCGCTGCCGTTCACCCCCTCCCGCACCTACAGCATTTCCCGGACa                       >  2:188259/1‑277 (MQ=255)
                                                                                                                                                                      cATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGCCAATCCCCTGCTCGGTGCCAAGGCCCTTCCCGGTGAAACCGACATCGCCCTGCCCTGTCCGCTGCCGTTCATCCTCTCCCGCCTCTACACCACTTccctcacc                       >  4:64305/1‑272 (MQ=14)
                                                                                                                                                                       aTTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACCTCCCCCTTCCCGGCCCGCTGCCGTTCACCCCCTCCCGCACCTACAGCCTTTACCGGa                         >  3:298370/1‑274 (MQ=255)
                                                                                                                                                                                ataGGAAAGGGCTTTACGATGAGCGGAAAGCCAGCGGCGCGTCAGGGGGACATGGCGCAGTATGGCGGTAGCAGTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCCTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAAGCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGcc               <  3:234409/275‑1 (MQ=11)
                                                                                                                                                                                 taGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCcgc            >  2:180806/1‑277 (MQ=37)
                                                                                                                                                                                 taGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCAACCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGACACCGACCTCGCACTGCCCGGCCCGCTGCCGTTCATCCTCTCCCCCCCCTACCGCAGTTACCCGACAAAACCGcc               >  3:413052/1‑274 (MQ=11)
                                                                                                                                                                                            tttACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCCCTCCCGCACCTACAGCAGTTACCGGCCAAAAACGCCCGCGCCGGTggg   >  1:226918/1‑275 (MQ=17)
                                                                                                                                                                                              tACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCCCTGCCGTCCATCCTCTCCCGCCCCTTCAGCAGTTACCGGCCAAAAACGCCCGCGCCGGTgggg  >  1:116129/1‑274 (MQ=9)
                                                                                                                                                                                              tACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCACCGGCGGTGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCCTCCCCTCCCGCACCTACAGCAGTTACCGGACAAAACCCCCCGCGCCGGTgggg  >  1:709516/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                
TGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGG  >  NC_000913/3618999‑3619461

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