Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 897,242 T→G 10.4% G53G (GGT→GGG)  potI → putrescine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913897,2420TG10.4% 145.5 / 2.0 67G53G (GGT→GGG) potIputrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base G (7/4);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGC  >  NC_000913/897004‑897465
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
gCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGt                                                                                                                                                                                                                                      <  1:858643/239‑1 (MQ=255)
     catcatGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGAGGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGGGCCGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGGGtt                                                                                                                                                                                                                                 >  1:768295/1‑239 (MQ=255)
     catcatGTTGCTGCTGCTAATTGTGCCGATAATATGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTTTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGGAGATGGTGCCGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAggt                                                                                                                                                                                                                                 >  1:244210/1‑237 (MQ=255)
              gctgctgctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCTTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAAGGCTGATGCTGGTTTTCTATTCGTTTAACAGCTCGAAGCTGGTGCAGGTGTGGGCCGGCTGGGCACGGCGCTGGGATGGGGGGTTATGGggcgg                                                                                                                                                                                                                        >  2:626070/1‑238 (MQ=255)
              gctgctgctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGCTGCTGGTTATCTATTGCTTTAACAGCCCGAAGCGGGTGCCGGGGTGGGCCGGCTGGTCACCGCGCCGGCAGGGGGAGTTTTTGcgcgg                                                                                                                                                                                                                        >  2:865146/1‑238 (MQ=255)
                tgtttcaaaTTGGTCCGATAATGTGGTTCCAAAAACCCCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGcgcgatg                                                                                                                                                                                                                      <  1:764223/232‑1 (MQ=255)
                  ctgctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGcgcgatga                                                                                                                                                                                                                     >  1:877248/1‑238 (MQ=255)
                    gctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGCATGGTGAGTTATTGcgcgatgacgc                                                                                                                                                                                                                  >  1:939039/1‑239 (MQ=255)
                              ccGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGGGAGTTATTGCGCGATGACGGGGTGATGAGt                                                                                                                                                                                                        >  2:1061114/1‑239 (MQ=255)
                               cGATAATGTTGTTCCACAAAGCCCAGGAAAAAAGCGGGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTTGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGt                                                                                                                                                                                                        <  1:711774/238‑1 (MQ=255)
                                           tCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTTTGGTGAGTTATTGCGCGGTGACGGCTTGATGATGGGGGTTGGTTTaa                                                                                                                                                                                           >  1:845246/1‑239 (MQ=255)
                                            ccACAAACCCCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTC‑AGCGCCAATGCTGATGCGGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGCCGGCGGGGGCCGGCGGTCCAACGCCGGGGTATGGGGAGTGATTGCGCGAGGGGGGGATGAGGGGGGGGGTTGGTTTCAg                                                                                                                                                                                          >  2:650815/1‑239 (MQ=255)
                                            ccACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGCGTGCGGTGGGTTTAAg                                                                                                                                                                                          >  1:1024613/1‑239 (MQ=255)
                                            ccACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGGGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTCTTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCGGGCTGGTCAACGCGCTGGAATGGTGAGTGATTGGGCGATACCACGAGGAGGCGGGCGGTGGGTTTAAg                                                                                                                                                                                          >  2:105166/1‑239 (MQ=255)
                                               caAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAg                                                                                                                                                                                          >  1:1063893/1‑236 (MQ=255)
                                                 aaCCCCAGCAAAAAAGGGGTGGAAAACCGGGATGAATAATTTACCGGTGGTTCGTTCGCCCTGGGTGATTGTGATTTTGCGGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGa                                                                                                                                                                                     <  2:1024613/239‑1 (MQ=255)
                                                 aaCACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGa                                                                                                                                                                                     <  1:14895/239‑1 (MQ=255)
                                                  acacCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGCTCTGGTGAGTTTTTGCGCTATGACGCGGTGTGGGGTGCGGTTGGTTTatggttgac                                                                                                                                                                                    >  2:955964/1‑231 (MQ=255)
                                                       aGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTTTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGACGTTGAGTGCGGTTGGTTTAGGTCTGACAATTg                                                                                                                                                                               >  1:3026/1‑239 (MQ=255)
                                                       aGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTg                                                                                                                                                                               >  1:997750/1‑239 (MQ=255)
                                                           taagagtcagggAGAACACGGATGAATAATTTCCCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGc                                                                                                                                                                           <  1:165566/230‑1 (MQ=255)
                                                           aaaaaGGGGGGGAGAATCGGGGTGAATAATTTCCCGGTAGTGCGTTCGCCCTGGGGGATTTGGATTTTTCTGCTGGCCTTTACCTTTCACTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGc                                                                                                                                                                           <  2:718329/239‑1 (MQ=255)
                                                                   tGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGGTGATGAGTGCGGTTGGTTTAAGTCTGCCAATTGCGGCCTGTgcg                                                                                                                                                                    >  1:703333/1‑238 (MQ=255)
                                                                    gggaaaaaccGGATTAATTATTTACCGGTAGTTCGTTCGCCTTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTgcggc                                                                                                                                                                  <  2:759978/230‑1 (MQ=255)
                                                                               gatttataatTTCCCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGTGGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGAcggcggcgg                                                                                                                                                       <  2:939039/234‑1 (MQ=255)
                                                                                aTGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGCAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGGATGGTGGGTTATTGCGCGCTGACGCGGTGATGAGTGCGGTTGGGTTAAGTCTTACCATTGCGGCCTGTGCGGCGAcggtcgcttt                                                                                                                                                      >  2:165566/1‑232 (MQ=255)
                                                                                aTGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACggcggcggc                                                                                                                                                      >  1:759978/1‑239 (MQ=255)
                                                                                 ttaatttttttCCGGTTGTTCGTTCGCCCTGGGGGATTGGGATTTTGCTTCTGGGCTTTACCTTGCTCTC‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGgcggcggcg                                                                                                                                                     <  2:3026/232‑1 (MQ=255)
                                                                                        tttACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGGGGCGAGCCTCg                                                                                                                                              >  1:687491/1‑239 (MQ=255)
                                                                                          tACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGGTCTTCGGa                                                                                                                                            >  1:525876/1‑239 (MQ=255)
                                                                                          tACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGGGAGTTATTGCTGGGTGACGCGTGGTTGAGGGCGGTTGGTTTAAGGTTGACAATTGCGGCCTGGGCGGGGAcggcgggggggggccccgga                                                                                                                                            >  2:562127/1‑228 (MQ=255)
                                                                                          tACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGGGAGTTATTGCGCGATGACGCGATGATGAGGGCGGTTGGTTTAAGTTTGACAATTGGGGGCTTGgcggggggggggggggggggccgggga                                                                                                                                            >  2:599993/1‑216 (MQ=255)
                                                                                             cGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCACCGCGCTGGTATGGTGAGTTATTGCGCGATGCCGCGCTGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGGGGGGCGACGGCGGCGGCGGTTCTCGGaaat                                                                                                                                         >  2:87835/1‑237 (MQ=255)
                                                                                                 tgttcggttGCCCGGGGGGTTTGTGATTTGGCTGTTGGGGTTTACCTTTCTCAA‑CGCGCCAATTCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGACGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTg                                                                                                                                     <  2:768295/238‑1 (MQ=255)
                                                                                                     cgttcgGCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCt                                                                                                                                                <  1:744059/224‑1 (MQ=255)
                                                                                                     cgttcgCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGGGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTgcggc                                                                                                                                 >  1:67718/1‑239 (MQ=255)
                                                                                                       ttcgCCGTGGCGGATTGTGATTTTGCTGCTGGGCTTTTCCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGcggcgg                                                                                                                               <  1:501261/239‑1 (MQ=255)
                                                                                                             cTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGGGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGGGACGGCGGCGGCGATCCTCGGAACTATTGCGGCggtggggg                                                                                                                         >  1:616793/1‑236 (MQ=255)
                                                                                                               ggCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGGGTGCGGTTGGTTTAAGTCTGGCAATTGCGGCCTGGGCGGGGGCGGCGGGGGCGGGCCCCGGAAcatttctggttcttttaaac                                                                                                                       >  2:744059/1‑220 (MQ=255)
                                                                                                                   gATTGGGATTTTGCGTCTGGGTTTTCCCTTTCTCTA‑CGCGGCAATGCTGATGCTGGTTATCTATTCGTTTAACCGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGtggtgctggtg                                                                                                                    <  1:423677/238‑1 (MQ=255)
                                                                                                                      ttgtttttttttgtggggTTTTACCTTTTCCTA‑GGCGCCAATGCTGATGGTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATTCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCt                                                                                                                <  2:877248/224‑1 (MQ=255)
                                                                                                                              tgctgcTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGggtggctgggggccacttgcggg                                                                                                        >  1:822868/1‑225 (MQ=255)
                                                                                                                                   tttggtttttCCTTTTTCTA‑CGCGCCAATGGTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGAGGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTcg                                                                                                   <  2:997750/233‑1 (MQ=255)
                                                                                                                                   ctgGGCTTTCCCTTTC‑CTAGCGCGCCAATGTTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGT‑GGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTcgc                                                                                                  <  1:1068914/239‑1 (MQ=255)
                                                                                                                                     tgggtttACCTTTCTCTA‑CGGGCAAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGGGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTcgc‑g                                                                                                <  2:1059736/235‑1 (MQ=255)
                                                                                                                                            aCCTTTCCCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTGGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑aaa                                                                                        <  1:28584/239‑1 (MQ=255)
                                                                                                                                                  ctctA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGCGGGGCTGGGGCGCTTTGTCTGGTTTTGC‑GGGCC‑CAAAGGCtt                                                                                  >  1:128563/1‑239 (MQ=255)
                                                                                                                                                    ctA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGGTGAACAGCGTGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTg                                                                                <  1:543226/239‑1 (MQ=255)
                                                                                                                                                       tcgcgCCAATGTTGATGTTGGTTTTTTTTTCTGTTAACGGCTCGAGGCTGGTGACGGGGTGGGCCGGCTGGTCAACGCGCTGGTTTGGTGAGTTATTGCGGGATGACGGGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGcc                                                                              <  2:251276/238‑1 (MQ=255)
                                                                                                                                                        cgcgCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGGTGGGTGGCGTTGGGAGGATTTGCGGGTTC‑AAAGGGGTTTg                                                                                >  1:468204/1‑237 (MQ=255)
                                                                                                                                                        cgcgCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGCTGGGGGGCTTTGGCAGGTGTCTC‑gggttcaataggctgtg                                                                                >  1:233515/1‑222 (MQ=255)
                                                                                                                                                        cgcgCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGGGGGGCTGGTGGGCCTTTGCGGGGTTTGC‑GGGTC‑CAAAGGGttttc                                                                               >  1:251276/1‑235 (MQ=255)
                                                                                                                                                         gcgcCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGTGGCTGGGGGGCTTTGGGAGGGTTTCC‑CGATC‑AACAGGCttttgcct                                                                            >  1:412448/1‑234 (MQ=255)
                                                                                                                                                                                        ttaacgggctgagggggggtgcGGGGTGGGCCGGGTGGTCATGGCGTTGGTTTGGTGAGTTATTGCGGGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGAt                                          <  2:1063893/218‑1 (MQ=255)
                                                                                                                                                                                               cTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGGGGGCTTTTTCAGGTTTTCC‑GGGTC‑AAAATGCTTTTCCTTTTTGTTTCACC‑CGCCCCCG‑GGGATT‑CCAGATGttctc                                    >  1:804184/1‑234 (MQ=255)
                                                                                                                                                                                                       tGGTGACGGTGGGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGATGTCATCACGGGCt                             <  1:562127/237‑1 (MQ=255)
                                                                                                                                                                                                        gggtgcGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGTTTGGTTTAAGTCTGACAATGGCGGCCTGTGCGGGGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGATGTCATCACGGGCTTGt                          <  1:87835/234‑1 (MQ=255)
                                                                                                                                                                                                         gTGACGGTGTGGGCCGGTTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGATGTCATCACGGGCTTGt                          <  1:786893/238‑1 (MQ=255)
                                                                                                                                                                                                             aggggtgGGCCGGCTGGTCAACGCGGTTGTCTGGTGCGTTATTGCGCGATGACGGGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGATGTCATCACGGGCTTGTCGCtg                     <  1:187780/235‑1 (MQ=255)
                                                                                                                                                                                                               gggtgGGGCGGCGGGGCGCCGCGGTGGTTGGGGGAGTTGTTGCGGGGTGACGGGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGGGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGATGTCATCACGGGCTTGTCGCtgt                    <  2:128563/236‑1 (MQ=255)
                                                                                                                                                                                                                        gggtGGTAAACGGGGTGGTTTGGGGGGTTAGTGCGGGATGACGGGATGGTGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGATGTCATCACGGGCTTGTCGCTGTTGTTattatt          <  2:935548/236‑1 (MQ=255)
                                                                                                                                                                                                                          cTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGCGGGTGGGCTTTGGGCGGGTTCGC‑GGGTCAAAATGGCTTGGCTTTTTTG‑TTC‑CCCCGCCGCAG‑GGAAGGCCCAGATGTCTTCCCGGGCTTTTGGC‑GTTGtttttttttt        >  2:1035099/1‑237 (MQ=255)
                                                                                                                                                                                                                          cTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGGGTTGGGGGTGGGCGTTGGGTGGGTTTGC‑TGATC‑CAACAGGTTTCGCTTT‑TT‑ATTACCCCCCCCCCGCGTAATG‑GCAGATGTCACCTCggggtgttttctgttggttttttttc        >  2:604018/1‑216 (MQ=255)
                                                                                                                                                                                                                           tGGTCACGGGGCTGGTTTGGTGGGTTATTGCGCGATGAGGCGATGAGGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATtcgt       <  2:616793/239‑1 (MQ=255)
                                                                                                                                                                                                                            gggcAACGCGTTGGTATGTGGAGTTAGTGCGCGATGACGCGATTGTGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATtcgt       <  1:955964/235‑1 (MQ=255)
                                                                                                                                                                                                                                aaCGCGCTGGTAGGGGGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCgcgc  <  1:230772/239‑1 (MQ=255)
                                                                                                                                                                                                                                     gcTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAg                                            <  1:390864/192‑1 (MQ=255)
                                                                                                                                                                                                                                       tGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGATGTCATCACGGGCt                             <  1:212876/205‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
GCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTA‑CGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGC‑GGATC‑AAATGGCTTTGCCTTTATG‑ATCACCGCGCCGCTG‑GTAATG‑CCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGC  >  NC_000913/897004‑897465

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