Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 896,963 A→G 14.0% E277G (GAG→GGG)  potH → putrescine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913896,9630AG14.0% 172.0 / 2.1 72E277G (GAG→GGG) potHputrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base G (7/4);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTATGGTACTGCCGATTTATACCGCGTTGATTCGTATTGATTATTCGCTGGTGGAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAA  >  NC_000913/896725‑897199
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
tttATGGTACTGCCGATTTATACCGCGTTGATTCGTATTGATTATTCGCTGGTGGAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCCCGGTGGCCCGGACCGCATCATGATCGGGCGCGTGCAATGGCAaga                                                                                                                                                                                                                                                  >  1:896440/1‑239 (MQ=255)
  tATGGTACTGCCGATTTATACCGCGTTGATTCGTATTGATTATTCGCTGGTGGAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACGCCATCATGATCGGGCGCGGGCTATGGCAAGAGt                                                                                                                                                                                                                                                >  1:158378/1‑239 (MQ=255)
          tGCCGATTTATACCGCGTTGATTCGTATTGATTATTCGCTGGTGGAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGCTTATTGCCGGATCGATGCTGGGGTTTATACCGGCTGTGGGCGAGTTTGGGATCCCGGAACTGGCCGGTGGCCCGGACCGCATCCGAATCGGGGGAGGGCTTTGGGAAGGGATTTTCaaa                                                                                                                                                                                                                                        >  2:717479/1‑238 (MQ=255)
          tGCCGATTTATACCGCGTTGATTCGTATTGATTATTCGCTGGTGGAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACGGCATTATGATCGGGCGCGTGCTATGGCAAGAGGTTTTTAAt                                                                                                                                                                                                                                        >  1:559260/1‑239 (MQ=255)
          tGCCGATTTATACCGCGTTGATTCGTATTGATTATTCGCTGGTGGAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCaa                                                                                                                                                                                                                                         >  1:832084/1‑238 (MQ=255)
                         cgTTGATTCGTATTGATTATTCGCTGGTGGAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGGCCGGACACCATCATGATCGGGCGGGAGCTATGGCAGAAGTTTTTCAATAAcccggcccgtccg                                                                                                                                                                                                                         >  1:716373/1‑228 (MQ=255)
                         cgTTGATTCGTATTGATTATTCGCTGGTGGAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGAGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGGGTTTGTGATCCCGGAACCGCCCGGTGGCCCGGGCAGCATCATGATCGGGGCCGTGCTATGGCAAAGGTTTTTAAAAAACcgcgcaaggcgg                                                                                                                                                                                                                         >  2:119056/1‑231 (MQ=255)
                          gTTGATTCGTATTGATTATTCGCTGGTGGAAGCAGAGCTGGATCTCGGTGCACGCTCGTTGAAAACGTTTTTTACTGTGATCGTGCTGCTGACTAAAGGTGGGATTATTGCCGGCTCGGTGATGGGGTTTATACCGGCTGTGGGCGAGGTGGTGCTCCCGGAACTGGTCGGTGGTCCGGGCAGCATCCTTATCGGGCGGGGGCTATGGCAATGGGTTTTTaaacgctcttaatacaca                                                                                                                                                                                                                          >  2:285357/1‑223 (MQ=255)
                                tcgtttttagttttCTCTGGTGGAAACAGGGGTGGTTTTCGGTGCACGTACGCGGAAAACGGTCTTTGCTGTGGTGGTGCCGCTGATTAAAGGTGGGATTTTTGCCGGATCGATGCTGGTGGTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACGGCTCGGTGGCCCGGACAGCATCATGATCGGGCTCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCt                                                                                                                                                                                                                  <  2:98246/227‑1 (MQ=255)
                                        attattCGCTGGTGGAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCt                                                                                                                                                                                                                  <  1:510774/231‑1 (MQ=255)
                                        attattCGCTGGTGGAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTa                                                                                                                                                                                                          >  1:635287/1‑239 (MQ=255)
                                             tCGCTGGTGGAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATAATGATCGGGCGGGTGATATGGCAAGGGTTTTTCAATAACCGAGACTGGCCGGTGGCCTcggggggcaggac                                                                                                                                                                                                     >  1:419869/1‑232 (MQ=255)
                                                  ggtggAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTTTGGCAAGAGTTTTTCAATAACCCCGCCTGGCCGGGGGCCCcggcggggagggtgtt                                                                                                                                                                                                  >  2:510774/1‑227 (MQ=255)
                                                             cgcTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGGGTCCCGGACAGCATCATACTCGGGCGCGTGCTCTGGCAAGAGTTTTTCAATAACCACGACTAGACGGTGGCCTcgggggggacggacgtaacgttgtacggg                                                                                                                                                                                     >  2:585357/1‑216 (MQ=255)
                                                                  gATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCGATGGCAAGGGTTTTTCAATAACCCCGCTTGGCCGGGGGCCCCGGCGGTAGGGGACCTCATGttcgttcgcttact                                                                                                                                                                                >  2:647933/1‑227 (MQ=255)
                                                                     ccggggtgccgcgctgtGAAACTGTTTTTTATTGGGTTCGTGCCGTTGACTAAAGGGGGGTTTATTGCCGGTTCGTGGTTGGTTTTTCTCCCGGCTGTGGGGTAGTTTGGGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGGGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATtgt                                                                                                                                                                             <  2:449261/223‑1 (MQ=255)
                                                                      cgctttgctcgccctttgaaaCTTTTTTTTCTGTTTTTGTGCTTCTGAATGAAGGGGGGATTATTCCCGGATCGATGCTGGTGTTTACCCCGGCTGGTGGCGAGTTTTTGATCCCGGAGCTGCTCGGTTGCCCGGACAGCAGCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATtgtg                                                                                                                                                                            <  2:360555/221‑1 (MQ=255)
                                                                       gggTGCACGACCCGGGAAATCTTTTTTTTCTGTGATTGTGCCGCTGCCAAAAGGTGGGTTTTTTGCCGGATCGATGCTGGTGTTTACCCCGGTTGTGGGCGAGTTTGTGATCCCGGAACTGCGCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGTAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGc                                                                                                                                                                           <  2:642157/238‑1 (MQ=255)
                                                                       cGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGGAGCGATCATCATGTTGCTGCTGCTAATtgtg                                                                                                                                                                            >  1:528933/1‑238 (MQ=255)
                                                                        ggtgcagtgccGTTGAAATGGTTTTTTACGGTGTTCGTGGCGCTGGCTAAGGGTGGGTTTATTGCCGGATGGATTCTGGTGGTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTGGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGcc                                                                                                                                                                          <  2:37000/230‑1 (MQ=255)
                                                                           gcccggccGCTGAAAGCGTTCTTTACGGGGATCGGGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGa                                                                                                                                                                        <  1:119056/232‑1 (MQ=255)
                                                                            cACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGGGTTTTTCAATAACCGCGACTGGCCGGTGACCTCGGCGGTAGGGGTCATCATGTTGGTGCTGCTAATTGTGCCGaga                                                                                                                                                                      >  1:433442/1‑237 (MQ=255)
                                                                             aCGCCCGTTGAAATCGTTTTTTTCTGGGTTCGTGTCGCTGACTAAAGTTGGGATTATTGCCGGATCGATGCTGGGGTTTATCCCGGCTGTGGGCGAGTTTGTGAGCCCGGAACTGCTCGGTGGCCCGTACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCATCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATaa                                                                                                                                                                     <  2:201865/239‑1 (MQ=255)
                                                                                            gTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTACCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGACCTGCTCGGTGGCCCGGCCAGCCTCGTGATCGGGCGCGTGCTATGGCAAGGGTTTTTCAATACCCGCTACTGGCCGGTGGCCTCGGCGGACGCGGACATCCAGTTTTTGCTCCTAATTTTGcccctatgtgggggttccaaaa                                                                                                                                                      >  2:149562/1‑219 (MQ=255)
                                                                                                       tgATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGAGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCCGGGGCCTCGGGGGTAGCGATCATCATGTTGCCGCTGCGAAATGTGGCGGTAATGTGGTTCGCCAAACACACGCacaaa                                                                                                                                          >  2:257208/1‑239 (MQ=255)
                                                                                                           tgttccGGTGGCTAAAGGGGGGTTTTTTTCCGGGTCGTGGCTGGTGTTTTGCCCGGCTGGGGGCGGGTTTGTGATCCCGGAACTGCCCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCg                                                                                                                                     <  2:744905/235‑1 (MQ=255)
                                                                                                                     gcttagggTGGGATTTTTGCCGGATGTATGCTGGTGTTTATCCCGGCTGTGGGCGGGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGGGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAaca                                                                                                                          <  2:728716/233‑1 (MQ=255)
                                                                                                                     attaaggggggtttttttGCCGGATCGTGTTTGGTTTTGACCCCGGTTTTGGGCGTGTTTGTGATCCGGGAATGGCTGGGTGGCCCGGACAGCATCATGTTCGGGCTGGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAtca                                                                                                                          <  2:433442/223‑1 (MQ=255)
                                                                                                                      ctaagtgggggATTTTTGCCGGAGCGGTGCTGGTGTTTACCCCGGGTGTGGGCGATTTTGTGATCCCGGAACTGCGCGGGGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATATCCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAaca                                                                                                                          <  2:559260/230‑1 (MQ=255)
                                                                                                                      cTAAAGGTGGGATTATTGCCGGAGCGATGCTGGTGTTTATCCCGGCTGGGGGCTAGTTTGGGATCCCGGAACGGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAaca                                                                                                                          <  1:647933/238‑1 (MQ=255)
                                                                                                                               ggATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGGTCATCAAGTTGCTGCTGCTAATTGTGCCTATAATGTGGTTCCACAAAAACCCGCAAAAAAGCGTGGGAGaaaccgggggagta                                                                                                                >  1:418326/1‑227 (MQ=255)
                                                                                                                               ggATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGCTCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGGTCCACAAACACCAGCAAAAAAGCGTGGGGGAAACCGGTTGaaaa                                                                                                                >  1:860038/1‑237 (MQ=255)
                                                                                                                                   tatgGCCGGATGGGTTCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAAGTGCTCGGTGGCCCGGACAGCATCAGGATCGGGCGGGTGCTATGGAAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAAttt                                                                                                            <  1:310393/239‑1 (MQ=255)
                                                                                                                                    attGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGGTAATTTGGTTCCACAAAAACCAGCAAAAAAGCCTGGGAGAACACGGATGAAAAAttt                                                                                                            >  1:855309/1‑238 (MQ=255)
                                                                                                                                          ggttcgtggttGGTGTTTATCCCGGATGTGGGCGGGTTTGTGATCCCGGAACTGCCCGGTGGCCCGGACAGCATCATGTTCGGGCGCGTTCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTa                                                                                                     <  2:896440/229‑1 (MQ=255)
                                                                                                                                                 tGCTGGTGTTTTCCCCGGCTGTGGGCGATGTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAgttcgtt                                                                                              <  2:253179/239‑1 (MQ=255)
                                                                                                                                                 tGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGACCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAgttcgt                                                                                               <  1:585357/238‑1 (MQ=255)
                                                                                                                                                       tgtTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGAGATCATCATGTTGCGGCTGCTAATTTTGCCGGTAAGTTGGTTCCCCAAAACCCAGCAAAAAAGCGAGGGAGGAACAGGAGGAATAAttttacggaggttggttgacccc                                                                                        >  2:538831/1‑219 (MQ=255)
                                                                                                                                                         tttATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGCGAACACGGATGAAAAATTTACCGGTAGTTCGTTCGCCCTgg                                                                                      >  1:253179/1‑239 (MQ=255)
                                                                                                                                                           tATCCCGGCTGTGGGCGAGTTTGTGATCCCGGACCTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGTGTTTTTCAATAAACGCGACTGGCCGGCGGCCTCGGCGGTAGTGATCCTCCTGTATCTGCTGCTGCTTGTGGCGATAACTGTGTTGAGACACCCCCAACAACAACGCGCGGGGGAAACCGGATGAATAATTTACCGGTTgttagttcccacatgcg                                                                                    >  2:398697/1‑230 (MQ=255)
                                                                                                                                                              cccGGCTGTGGGCGATTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGAt                                                                                 <  1:562850/239‑1 (MQ=255)
                                                                                                                                                                  gCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCCTCATGTTGCTGCTGCTAATTGGGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAAAACGGATGAATAATTTAACGGTAGTTCGTTCGCCCTGGGGGAGtggg                                                                             >  2:809430/1‑237 (MQ=255)
                                                                                                                                                                      tgGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTTTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGGTGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGGGAttt                                                                         >  1:864067/1‑239 (MQ=255)
                                                                                                                                                                        tgCGTGTTTGTGATCCCGGGACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTGCAATAACCGCGACTTGCCGTGGGCCGCGGCGGCAGAGGTCCTCAAGTGCCTGCTGCTACTTGTGCCGGTAATGTGGGTCCCCAAACACCAGCCAAAAAGAGCGGGAGAAAACCAGTTGAGAATTAACCCGGCTTTCGTTCCCCCGGGCCGAGTGGGATTTtg                                                                       >  2:517366/2‑239 (MQ=255)
                                                                                                                                                                              tttggtctCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGAGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTgctgct                                                                  <  1:789347/231‑1 (MQ=255)
                                                                                                                                                                              tttGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTgctgct                                                                  <  1:257208/238‑1 (MQ=255)
                                                                                                                                                                                tgtgATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGAGGAATAATTTACCGGAAGTTCGTTCGCCCTGGCGGGTTGGGATTTTGCTGCTggg                                                               >  1:519367/1‑239 (MQ=255)
                                                                                                                                                                                      cccGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGGAGGGATCATCCTGTTGGTGCCGCTAATTGTGCCGATAAAGGGGGTCCACAAAAACCGGCAAAAAAGCGTGGGGGAAACCGGATGAATAAATAAACGGTAATTGGTTCGCCCTGGGGGGTTGGAGTTTTTGTGCTGGGGttt                                                           >  2:493809/1‑237 (MQ=255)
                                                                                                                                                                                       ccGGACCTGTCCGGGGGCCCGACCAGCATGATGGTCGGGCGCGTTCTATGGCAAGGGTTTTTCAAAACCCGCGCCGGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTAcc                                                        <  2:519367/239‑1 (MQ=255)
                                                                                                                                                                                          gactgtctcgtgtgCCCGGGCAGCTTCATGGTCGGGGGGGTGGTATGGGAAGGTTTTTTCAATAACCGCGGCGGGGCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCttt                                                     <  2:158378/226‑1 (MQ=255)
                                                                                                                                                                                           aaCTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCgg                                                                                                                                                                                                                        >  1:346390/1‑78 (MQ=255)
                                                                                                                                                                                           aaCTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCgg                                                                                                                                                                                                                        <  2:346390/78‑1 (MQ=255)
                                                                                                                                                                                               gcccGGTGGCCCGGACAGCATCATGATCGGGCGGGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTa                                                <  1:225549/236‑1 (MQ=255)
                                                                                                                                                                                                cTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGCTAATGTGGTTCCACAAACACCAGCAAAAAAGCGGGGGGGAACACGGATGAATAATTTACCGGGAGTTCGTTCGCCCGGGGGGGTTGTGGTTTTTTTGGTGGGGTTATCCttttctcac                                               >  2:225549/1‑233 (MQ=255)
                                                                                                                                                                                                cTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGTGTCCACAAACACCAGCAAAAAAGCGAGGGGGAAACAGGATGGATAATTTACCCGTAGGTCGTTCCCCCCGGGGGGTTGTGATTTTTGTTCTGGGCCTTACCTTTctctcc                                               >  2:637339/1‑237 (MQ=255)
                                                                                                                                                                                                cTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGATTAATTTTACGGGAATTCGTTCGCCCTGGGGGATTGGGATTTTTCTGCTGGGGTTTACCTTTATCTAc                                               >  2:129548/1‑239 (MQ=255)
                                                                                                                                                                                                cTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGTTGGGCTTTACCTTTTTCTa                                                >  1:154661/1‑238 (MQ=255)
                                                                                                                                                                                                cTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGGGGGTCCCCAACCACCAGCAAAAAAACGGTGGGGAAACAGGAGGGATAATTTACCCGAAGATCGTTCGCCCCGGCGGATTGGGGATTTGCTTTGGGGCTTTACATTTttctcg                                               >  2:296960/1‑237 (MQ=255)
                                                                                                                                                                                                     tGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCCGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACCCCAGCAAAAAAGCGGTGGAGAACACGGATGAAAAATTTACCGGACGTTCGTTTGCCCTGGCGGGTTGGGATTTTGTTGCTGGGCTTTACCTTTTTCTACGCGcc                                          >  2:738787/1‑239 (MQ=255)
                                                                                                                                                                                                      ggCCCGGACGGCATCATGACGGGGCGGGTGTTAGGGCAAGAGTTTTGCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGcc                                          <  2:864067/238‑1 (MQ=255)
                                                                                                                                                                                                               cgCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAatgctgatg                               <  1:129548/238‑1 (MQ=255)
                                                                                                                                                                                                               aGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAatgctgatg                               <  1:809430/239‑1 (MQ=255)
                                                                                                                                                                                                                  tttgtgggcggggggtggttATGGCAAGAGTTTTTCAATAACCGGGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAatgctga                                 <  2:635287/215‑1 (MQ=255)
                                                                                                                                                                                                                   ttttGATCGGGCCCGTGCTATGACACGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGGAACGGTCATCATGATGTTGCTGCAAATTGGGCCGATAAGGTGGTTGCACCAACCCCCAGATAACACCGTAGAAGAACAAGGGTTAAGAATGTACCGGTAGCTCGCTCCCCCTGGGGGATTTGGATTTTGCGTCTGGGGTTTCCCTTTTCCCCGGCGCCGaggctgatggggt                           >  2:501658/4‑235 (MQ=255)
                                                                                                                                                                                                                         tCGGGCGTGTTTTAGGGCAAGAGTTTTTCAATGCCCGCCGCTGGGCGGTGGCCTCGGCGGTAGCCATCATCATGTTGCTGCTGCTAATTTGGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCt                     <  2:731007/239‑1 (MQ=255)
                                                                                                                                                                                                                           ggggccgTGCTTTGGTATGAGTTTTTCAATAACCGGGGCTGGCCGGTGGCCTCGGGGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTAt                   <  2:441527/234‑1 (MQ=255)
                                                                                                                                                                                                                           gggctgggtcTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACCCCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTAt                   <  1:493809/230‑1 (MQ=255)
                                                                                                                                                                                                                           gggCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTTCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGTTGCTGCTGGTTATCTAt                   >  1:731007/1‑239 (MQ=255)
                                                                                                                                                                                                                            ggCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTTTCttt                   >  1:774377/1‑236 (MQ=255)
                                                                                                                                                                                                                               gcgTGCTATGGCAAGGGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCg                <  1:333192/238‑1 (MQ=255)
                                                                                                                                                                                                                                  tGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGCTCCACAAACACCCGCCAAAAAGCGTGGGAGAACACGGATGGATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTTGTTGTGGGCCTTACCTTTTCCTACCGTCCAatgtggttgctgttgttttctttgtttt            >  2:500210/1‑223 (MQ=255)
                                                                                                                                                                                                                                        ggaaaGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGGGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCt      <  1:717479/236‑1 (MQ=255)
                                                                                                                                                                                                                                           aaGGGTTTTTCAATACCCGGGGCGGGCCGGGTGCCTCGGGTGTAGCGATCAGCATGTTGCCGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACCGCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCg    <  2:528933/238‑1 (MQ=255)
                                                                                                                                                                                                                                            agagTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTTTGCCGATAATGTGGTTCCACAAAAACCAGCAAAAAATCGGGGGAGAACACCGGTGAATAATTTCCCGGTAGTTCGTTCGCCCTGGTGGATTGTGATTTTCCTCCTCGGCTTTACCTTTCTCCCCGCGCCAATGCGTGAGCTTGTTATTTTTTCTTTTAACAGTCCGaa  >  2:593637/1‑239 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
TTTATGGTACTGCCGATTTATACCGCGTTGATTCGTATTGATTATTCGCTGGTGGAAGCAGCGCTGGATCTCGGTGCACGACCGCTGAAAACGTTCTTTACTGTGATCGTGCCGCTGACTAAAGGTGGGATTATTGCCGGATCGATGCTGGTGTTTATCCCGGCTGTGGGCGAGTTTGTGATCCCGGAACTGCTCGGTGGCCCGGACAGCATCATGATCGGGCGCGTGCTATGGCAAGAGTTTTTCAATAACCGCGACTGGCCGGTGGCCTCGGCGGTAGCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAA  >  NC_000913/896725‑897199

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