Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,154,042 T→G 15.5% S160R (AGT→AGG yghB → general envelope maintenance protein; DedA family inner membrane protein; putative multidrug efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,154,0420TG15.5% 147.0 / 4.0 74S160R (AGT→AGGyghBgeneral envelope maintenance protein; DedA family inner membrane protein; putative multidrug efflux transporter
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base G (10/2);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 5.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCAT  >  NC_000913/3153809‑3154276
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
tgtgggttAATTTATATTCGGGGGGGTTGGTTAGGAAATACCAAAGCGGAGAAGGGCTGGCTGGCACAGTTTCCTGCTAAATATCACCAGCGCGCCCCCTGCATGTTTGCCGGCCCCGGTCTGCTGGGGTGGCTGGCTGGACGTTTTCTTGCATTTGTCCTTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGAt                                                                                                                                                                                                                                          <  1:867200/231‑1 (MQ=255)
tggtgGCTAATTTATATTCAGGGGGCGTGGTTAGGGAATCCAAAAACGGTGAAAGGCTGGCGGGCACGGTTTCCTGCAAAATATCACCAGCGCGCCTCCTGCATGTTTGACCGCCACGGTCTGGTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGAt                                                                                                                                                                                                                                          <  1:129347/235‑1 (MQ=255)
  ctgGCTAAGTTATATTCAGGGGCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCCGTTTTTCAACTGGTTAAGTGGAtt                                                                                                                                                                                                                                         >  1:855095/1‑237 (MQ=255)
          tttttttgggggggCGCTGTTGGGGAAAACCCAACCGGGGAAAGGGTGGGTGGTCACGGCTCCGTTCTAATTCCTCCCGGCGCGCCCCGTGTTGTTGTGCCCCCCGCGGTCTGCGGGCGCTGCTGGCGGGCCGTTTTCTGGACTTGGCCCGTACGCTGCTGCCACCCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTc                                                                                                                                                                                                                               <  2:748314/229‑1 (MQ=255)
               atTCAGGGGCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGCCTGCCAAACCGCCGCTTCCCGTTTTTCAACTGGTTAAGGGGATTGCTGTGGGGCAGCGt                                                                                                                                                                                                                          >  1:386451/1‑239 (MQ=255)
                     gggCGCTTGGTTGGGAATTCCCAAACGGGGGAAGGGCTGGGGGGCCCGGCTTCCGGCAAAATCCCCCAGGCGCGCCCCCTGCATGTTTACCCGCCACGGTCTGCGGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTaa                                                                                                                                                                                                                     <  2:931915/239‑1 (MQ=255)
                     gggCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGGATTTCCGGTCTGCCCAACCGCCGCTTCCAGTTTTTCAACTGGATAAGTGGATTGCTGTGGGTCAGCGGGGGaa                                                                                                                                                                                                                     >  1:1030129/1‑238 (MQ=255)
                     gggCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTTCCAAACCGCCCCCTCCCGTTTTTTCACGGGTTAAATGGATTGCTTGGGGTCCACGTGGTaaa                                                                                                                                                                                                                    >  2:294580/1‑238 (MQ=255)
                     gggCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCCGTTTTTCCACTTGTTAAGGGGATTTCTGTGGGTCAGCGGGGGaa                                                                                                                                                                                                                     >  1:1059795/1‑238 (MQ=255)
                     gggCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCCAGTTTTCAACTGGTTAAGTGGATTTCTGTTGGTCAGCGGGGAAAc                                                                                                                                                                                                                    >  1:263860/1‑239 (MQ=255)
                     gggCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCCCTTCCCGTTTTTTCACTGGTTAAGGGGATTGCTGTGGGTCAGCGGGGGAAc                                                                                                                                                                                                                    >  2:219913/1‑239 (MQ=255)
                     gggCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCCCTTCCAGTTTTTCAACTGGTTAAAGGGATTGCTTTGGGGCACCGGGGTaa                                                                                                                                                                                                                     >  1:337770/1‑238 (MQ=255)
                             gTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGGATTTCCGGTCTGCCACACCCCCGCCTCCCGGTTTTCAACTGGGTTAGGGGGTTGGTTGGGGTCGGCCGGGTAACCCGttttt                                                                                                                                                                                                            >  2:618442/1‑238 (MQ=255)
                                 gggAAACCCCAAACGGTGAAAGGCGGGGCGGCCCCTGTTCCTGGAAAAATTCCCCCGGCGCGCCCCCTCCTTTTTGGCCCGCCGGGGTCTGCTGGCGCTGTGGGCTGGACGTTTTCTTGCTTTTGCCCGTACGCTGCTGCCACCCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCt                                                                                                                                                                                                         <  2:1066907/239‑1 (MQ=255)
                                 ggaaaTCCCAAACCGGTGAAAGGTGGGTGGGCACGGCTTCCTGCAAAATTCCACAGGCGCGCCCCCGGCATGTTTGCCCGCCACGGTTTGCGGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTa                                                                                                                                                                                                        <  1:70823/236‑1 (MQ=255)
                                  ggAAAACCCAAACGGGTGAAGGGCGGGCTGCCACAGTTTCCTGCTAAATACCCCCAGCGCGCCCCCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTa                                                                                                                                                                                                        <  1:93713/239‑1 (MQ=255)
                                  gaaaTCCCAACACGGGAAAGGGTTGGTGGCCCAGCTTCCTGGTTAAATTCCCCCGGGGCGCTCCTGGTTTTTTGGCCCCCCACGGTCTGCGGGCGGTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTa                                                                                                                                                                                                        <  2:510718/236‑1 (MQ=255)
                                     aTACCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGGCTGCCAAACCGCCCCCTCCAGCTTTTCAACTGGTTAAGTGGAGTTGCGTGGGGCAGCGTGGTAAACAGGTTTTGCTAtt                                                                                                                                                                                                      >  2:156904/1‑236 (MQ=255)
                                       aCCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCCAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTTCTTTgg                                                                                                                                                                                                                                  >  1:366050/1‑207 (MQ=255)
                                       aCCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCCGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCCGTTTTGGCTATGCCtt                                                                                                                                                                                                 >  1:357590/1‑239 (MQ=255)
                                       aCCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCCCCGCTTCCAGTTTTTCAAATTGTTAAGTGGATTGCTCTGGGTACACGTGGTAACCCAGTTTGGCATTGCCtt                                                                                                                                                                                                 >  1:125198/1‑239 (MQ=255)
                                       aCCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGGTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCTTACGCTGCTGCCAACCATGGCGGGGATTTCCGGACTGCCAAAACGCCGCTTTCCGTCTTTTCACTGGTTAAATGGATTGCTGTGGGTCACCGGGGTAAACAGTTTTGGCAATGctta                                                                                                                                                                                                 >  2:137938/1‑236 (MQ=255)
                                           aaaCGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGGGGGTTGCTGTGGGTCAGCGTGGGAACCCGTTTTGGCTTGGCCTTaaag                                                                                                                                                                                             >  1:1065461/1‑237 (MQ=255)
                                                            gctggcACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCCAACCGCCGCTTCCAGTTTTTCAACTGGGTAAGGGGATTGCTTTGGGTCAGCGGGGTAACCAGTTTTGGCTAGGCCTTAAGTTTGGTTCCGTTCGTTa                                                                                                                                                                             >  1:434553/1‑239 (MQ=255)
                                                                gccacaGTTTCCGGTAAAATATCCCCGGCGCGCCCCCGGCTTGTTGGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTaaa                                                                                                                                                                           <  1:11034/236‑1 (MQ=255)
                                                                     gggttCCTGATAAATTCCACCAGGGCGCCACCTGCTTGTTTGCCCGCCACGGTCTGCTGGCGTTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTCCGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCAt                                                                                                                                                                     <  1:548005/236‑1 (MQ=255)
                                                                       cTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTaaagc                                                                                                                                                                         >  1:286346/1‑231 (MQ=255)
                                                                            tGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGGATTTCCGGTCTGCCACAACGCCGCTTCCAGCTTTTCAACGTGGTTATGGGGTTTCTGTGGGGCAGCCCTGTAAACCGTTTTTGCTATGCCCTTAAGAGTCTTTCCTTCCTTAAACACCCGGAAGaaa                                                                                                                                                              >  2:769039/1‑237 (MQ=255)
                                                                            tGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCACTTTTTTACATGGTTAAGGGGGTTGCTGTGGGTCCGCGGGGTAACCACATTTGGCTATGCCTTAAGGATGATTCCGTTCGTTAAACGCCCTGACGAta                                                                                                                                                              >  2:827492/1‑238 (MQ=255)
                                                                                 aaaTATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTTCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCCGCGTGGGAACCAGTTTTTGCTATGCCTTAAGTATGATTCCGTTCGTTAAAACCCATGAAAATCAggg                                                                                                                                                          >  2:845742/1‑238 (MQ=255)
                                                                                  aaTATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGCCTGCCAAACCGCCGCTTCCAGGTTTTCAACTGGTTAAGGGGATTGCTGTggg                                                                                                                                                                                                                                 >  1:657527/1‑168 (MQ=255)
                                                                                         ccAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGGGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGt                                                                                                                                                  >  1:804745/1‑239 (MQ=255)
                                                                                         ccAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTCTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGGTTGATTCCGTTCGTTAAAACCCATGAAATTCAGGTAATTACGt                                                                                                                                                  >  1:141169/1‑239 (MQ=255)
                                                                                         ccAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGAGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTTTGCCTTAAGTATGATTCCGTTCGTTAAACGCCCTGAAGATAAGGTAAGGCCGt                                                                                                                                                  >  1:19669/1‑239 (MQ=255)
                                                                                            gcgcgcCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCCAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAAACAGTTTTGGCCATGCCTTAAGTATGATTCCGTTCGTTAACAGCCATGAAGATCAGGGAATGGCGtt                                                                                                                                                 >  2:944636/1‑237 (MQ=255)
                                                                                            gcgcgcCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGGCTGGTAACCAGGTTTGGCTATGCCTTAAGTATGATTCCTGTCGTTAAACGCCATGCAGTTAAGGTAATgcgctt                                                                                                                                                 >  1:494525/1‑232 (MQ=255)
                                                                                             cgcgcCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGCTAAGTGGATTGCTTTGGGTAAGCGCGGGAACCAGTTTTTGTTATGACTTACGTATGATTCGCTTCCTTAAACTCCATGAAACTAAGGAAATGaactt                                                                                                                                                 >  2:444265/1‑232 (MQ=255)
                                                                                               cgcCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCtga                                                                                                                                            >  1:830412/1‑239 (MQ=255)
                                                                                                            tttGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATTAAGATCAGGGAATGACGTTCCTTAGGATCCTTCCAAt                                                                                                                               >  1:328113/1‑239 (MQ=255)
                                                                                                               tGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGGAATGACGTTCCTGATGATCCTGCCAAtt                                                                                                                              >  1:715153/1‑238 (MQ=255)
                                                                                                                 gACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAAAGCCATGAAGATCAGGGAATGACGTTCCTGATTATCCTGCCAATTg                                                                                                                             >  1:147049/1‑238 (MQ=255)
                                                                                                                       cACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGGATTTCCCGTCTGCCAAACCGCCGCCTCCAGTTTTTCAACTGGGTAAGTGGATTGCTGTGGGTCACCGTTGTAACCAGTTTTGGCTCTTCCTTAAGAATGTATCCCTTCGTTAAACGCCAGAAAGCTCCGGGAAGGAGGTTCCTCTCGATACTGCCCATGCCCTTGtt                                                                                                                      >  2:579808/1‑239 (MQ=255)
                                                                                                                             ctggtgGCGCTGCTGGCTGGCCGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGAAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCg                                                                                                                 <  1:219913/238‑1 (MQ=255)
                                                                                                                                      ttgtgggccggACGTTTGCTTCAATTGGTCCGTACGTTGTTGCCAACCATGGCGGGAATTGCCGGTCTGCAAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGt                                                                                                        <  1:174984/233‑1 (MQ=255)
                                                                                                                                      ctgctgGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGGGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCCGGTAATAACGTTCCTGATGATCCTGCCAATTGCCTTTTTAACCGCCGGCTtttt                                                                                                       >  1:351673/1‑239 (MQ=255)
                                                                                                                                          gggctggACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCGGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAgg                                                                                                    <  1:618442/237‑1 (MQ=255)
                                                                                                                                            gttggcttttttCTTTCTTTTTTCCCTTCGGTTCTCCAAACCTGGGGGGAATTTCCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCAc                                                                                                 <  2:916354/231‑1 (MQ=255)
                                                                                                                                             ctggACGTTTTCTTGCATTTGTCCGTACGCTGGTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCAc                                                                                                 <  1:827492/238‑1 (MQ=255)
                                                                                                                                              tggctgTTTTTTTGCATTTGTCCGTACGCTGTTGCCACCCATGGCGGGAATTTCGGGTCGGCCAATCCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACg                                                                                                <  1:1043603/233‑1 (MQ=255)
                                                                                                                                              tggACGTTTTCTGCCATTTGCCGTTACCTGGGTTCCAACCATGGCGGGAATTTCGGGTCGGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACg                                                                                                <  1:916512/238‑1 (MQ=255)
                                                                                                                                                                 gTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGTTTGCTGTGGGTTAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCTTTAAACGCCATGAAGATCAGGGAATGACGTTCCTGTTGATCCTGCCCATTGCCTTGTTAACCGCCGGCTTGTTTGGCACCGTGTTTGTGGGGTTTacaaa                                                                           >  1:633393/1‑239 (MQ=255)
                                                                                                                                                                 gTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTACCTTGTTAACCGCTGTCTTTTTTGGGCCGCGTTTTGGGGGTGTTTaaaa                                                                            >  2:248335/1‑239 (MQ=255)
                                                                                                                                                                  tCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCAAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTg                                                                                                              <  1:763328/204‑1 (MQ=255)
                                                                                                                                                                           tgctgcCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAAACGCCGGCTTGTTAGGCACGCTGTTTGTGGGGATTAAAAAAAAACACTGt                                                                 >  1:293719/1‑239 (MQ=255)
                                                                                                                                                                              ctccAACCATGGGGGAAATTTCAGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTTCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTaa                                                               <  1:610002/236‑1 (MQ=255)
                                                                                                                                                                                          cGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTTCCTTGTTAACAGCTGGCTTGTTAGGCACGCTGTTtg                                                                                        >  1:242897/1‑201 (MQ=255)
                                                                                                                                                                                                          tGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGGTTAAAAAAAAAAACTGTAACGCCTGACGATTTTCCCCGTTCCCGGtg                                   >  2:980070/1‑238 (MQ=255)
                                                                                                                                                                                                           tcccaccccccgcTTCCGGTTTTTCAATTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATGCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTg                                  <  2:804745/230‑1 (MQ=255)
                                                                                                                                                                                                             aaacccgccgGTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTTGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCt                                <  1:248335/235‑1 (MQ=255)
                                                                                                                                                                                                             aaacccgccgCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCt                                <  1:845742/235‑1 (MQ=255)
                                                                                                                                                                                                                    ccgcTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTGCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGGTTGTGGGGATTAAAAAAAAATACCGTAACGCTTGACGATTTTCCCCGGTCCCGGGGGCTGGACggg                         >  2:816885/1‑239 (MQ=255)
                                                                                                                                                                                                                    ccgcTTCCAGTTTTTCAACTGGTTAAGAGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGTTTTTCCCCGTTCCCGTTTTGTGTACCgg                         >  1:819662/1‑239 (MQ=255)
                                                                                                                                                                                                                      gtttCCATTTTTTCAATTTTTTAAGTGGATTTCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGa                       <  2:147049/237‑1 (MQ=255)
                                                                                                                                                                                                                      gtttCCAGTTTTTCAACTGGTTAAGTGGTTTGCTGTGGGTCAGCGTGGTAACCAGTTTGGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGa                       <  2:1059795/237‑1 (MQ=255)
                                                                                                                                                                                                                      gtttCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAATGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGa                       <  1:294580/237‑1 (MQ=255)
                                                                                                                                                                                                                      gcTTCCATTTTTTCAATTGGTTAAGGGGATGGGTGGGGGTCAGCGTGGTAACCAGTTTTGGCTAGGCCTTAAGTATTATTCCGTTCGTTAACCGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCCGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGa                       <  2:293719/239‑1 (MQ=255)
                                                                                                                                                                                                                      ctttcatttttttCAACTGTTTAATTGGATTGGTGTGGGGCAGCGTTGTAACCAGGTTTGGCTATGCCTTAGGTAGTATTCCGTTGGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGa                       <  2:125198/231‑1 (MQ=255)
                                                                                                                                                                                                                        ttCCATTTTTTCAACTGGTTAATGGGATTGGTGTGGGTCAGCGTGGTAACCAGTTTTGGCTTTGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGTCGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAAc                     <  1:816885/239‑1 (MQ=255)
                                                                                                                                                                                                                           cAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAACCGCCATTAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTa                  <  1:444265/239‑1 (MQ=255)
                                                                                                                                                                                                                                tttCAACTTGTTAAGTGGTTTGCTGTGGGGCAGGGTGGTAACCAGTTTTGGCTATGCCTTAAGTTTGATTCCGTTCGTTAAACGCCATGAAGAACAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTaa             <  2:855095/239‑1 (MQ=255)
                                                                                                                                                                                                                                     aCTGGGTAAGTGGGTTTGTGGGGGTCGGCGTGGTAACCAGTTTTGGGTATGCCTTAAGTATGATTCCGTGCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTccc        <  2:633393/239‑1 (MQ=255)
                                                                                                                                                                                                                                         ggTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCGTTACCGGGAACGAATTTAATTCCCCtc     >  2:155722/1‑238 (MQ=255)
                                                                                                                                                                                                                                           ttAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCa   >  1:780939/1‑239 (MQ=255)
                                                                                                                                                                                                                                            tAAGTGGATTGCTGGGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCAt  <  1:650761/239‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
TGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCACAGCTTCCTGCTAAATATCACCAGCGCGCCACCTGCATGTTTGACCGCCACGGTCTGCTGGCGCTGCTGGCTGGACGTTTTCTTGCATTTGTCCGTACGCTGCTGCCAACCATGGCGGGAATTTCCGGTCTGCCAAACCGCCGCTTCCAGTTTTTCAACTGGTTAAGTGGATTGCTGTGGGTCAGCGTGGTAACCAGTTTTGGCTATGCCTTAAGTATGATTCCGTTCGTTAAACGCCATGAAGATCAGGTAATGACGTTCCTGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCAT  >  NC_000913/3153809‑3154276

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