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→GGGpotIputrescine 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

GCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGC  >  NC_000913/897004‑897465
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
gCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGt                                                                                                                                                                                                                                      <  1:858643/239‑1 (MQ=255)
     catcatGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGAGGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGGGCCGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGGGtt                                                                                                                                                                                                                                 >  1:768295/1‑239 (MQ=255)
     catcatGTTGCTGCTGCTAATTGTGCCGATAATATGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTTTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGGAGATGGTGCCGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAggt                                                                                                                                                                                                                                 >  1:244210/1‑237 (MQ=255)
              gctgctgctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCTTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAAGGCTGATGCTGGTTTTCTATTCGTTTAACAGCTCGAAGCTGGTGCAGGTGTGGGCCGGCTGGGCACGGCGCTGGGATGGGGGGTTATGGggcgg                                                                                                                                                                                                                        >  2:626070/1‑238 (MQ=255)
              gctgctgctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGCTGCTGGTTATCTATTGCTTTAACAGCCCGAAGCGGGTGCCGGGGTGGGCCGGCTGGTCACCGCGCCGGCAGGGGGAGTTTTTGcgcgg                                                                                                                                                                                                                        >  2:865146/1‑238 (MQ=255)
                tgtttcaaaTTGGTCCGATAATGTGGTTCCAAAAACCCCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGcgcgatg                                                                                                                                                                                                                      <  1:764223/232‑1 (MQ=255)
                  ctgctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGcgcgatga                                                                                                                                                                                                                     >  1:877248/1‑238 (MQ=255)
                    gctAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGCATGGTGAGTTATTGcgcgatgacgc                                                                                                                                                                                                                  >  1:939039/1‑239 (MQ=255)
                              ccGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGGGAGTTATTGCGCGATGACGGGGTGATGAGt                                                                                                                                                                                                        >  2:1061114/1‑239 (MQ=255)
                               cGATAATGTTGTTCCACAAAGCCCAGGAAAAAAGCGGGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTTGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGt                                                                                                                                                                                                        <  1:711774/238‑1 (MQ=255)
                                           tCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTTTGGTGAGTTATTGCGCGGTGACGGCTTGATGATGGGGGTTGGTTTaa                                                                                                                                                                                           >  1:845246/1‑239 (MQ=255)
                                            ccACAAACCCCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTCAGCGCCAATGCTGATGCGGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGCCGGCGGGGGCCGGCGGTCCAACGCCGGGGTATGGGGAGTGATTGCGCGAGGGGGGGATGAGGGGGGGGGTTGGTTTCAg                                                                                                                                                                                          >  2:650815/1‑239 (MQ=255)
                                            ccACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGCGTGCGGTGGGTTTAAg                                                                                                                                                                                          >  1:1024613/1‑239 (MQ=255)
                                            ccACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGGGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTCTTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCGGGCTGGTCAACGCGCTGGAATGGTGAGTGATTGGGCGATACCACGAGGAGGCGGGCGGTGGGTTTAAg                                                                                                                                                                                          >  2:105166/1‑239 (MQ=255)
                                               caAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAg                                                                                                                                                                                          >  1:1063893/1‑236 (MQ=255)
                                                 aaCCCCAGCAAAAAAGGGGTGGAAAACCGGGATGAATAATTTACCGGTGGTTCGTTCGCCCTGGGTGATTGTGATTTTGCGGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGa                                                                                                                                                                                     <  2:1024613/239‑1 (MQ=255)
                                                 aaCACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGa                                                                                                                                                                                     <  1:14895/239‑1 (MQ=255)
                                                  acacCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGCTCTGGTGAGTTTTTGCGCTATGACGCGGTGTGGGGTGCGGTTGGTTTatggttgac                                                                                                                                                                                    >  2:955964/1‑231 (MQ=255)
                                                       aGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTTTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGACGTTGAGTGCGGTTGGTTTAGGTCTGACAATTg                                                                                                                                                                               >  1:3026/1‑239 (MQ=255)
                                                       aGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTg                                                                                                                                                                               >  1:997750/1‑239 (MQ=255)
                                                           taagagtcagggAGAACACGGATGAATAATTTCCCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGc                                                                                                                                                                           <  1:165566/230‑1 (MQ=255)
                                                           aaaaaGGGGGGGAGAATCGGGGTGAATAATTTCCCGGTAGTGCGTTCGCCCTGGGGGATTTGGATTTTTCTGCTGGCCTTTACCTTTCACTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGc                                                                                                                                                                           <  2:718329/239‑1 (MQ=255)
                                                                   tGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGGTGATGAGTGCGGTTGGTTTAAGTCTGCCAATTGCGGCCTGTgcg                                                                                                                                                                    >  1:703333/1‑238 (MQ=255)
                                                                    gggaaaaaccGGATTAATTATTTACCGGTAGTTCGTTCGCCTTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTgcggc                                                                                                                                                                  <  2:759978/230‑1 (MQ=255)
                                                                               gatttataatTTCCCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGTGGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGAcggcggcgg                                                                                                                                                       <  2:939039/234‑1 (MQ=255)
                                                                                aTGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGCAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGGATGGTGGGTTATTGCGCGCTGACGCGGTGATGAGTGCGGTTGGGTTAAGTCTTACCATTGCGGCCTGTGCGGCGAcggtcgcttt                                                                                                                                                      >  2:165566/1‑232 (MQ=255)
                                                                                aTGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACggcggcggc                                                                                                                                                      >  1:759978/1‑239 (MQ=255)
                                                                                 ttaatttttttCCGGTTGTTCGTTCGCCCTGGGGGATTGGGATTTTGCTTCTGGGCTTTACCTTGCTCTCCGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGgcggcggcg                                                                                                                                                     <  2:3026/232‑1 (MQ=255)
                                                                                        tttACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGGGGCGAGCCTCg                                                                                                                                              >  1:687491/1‑239 (MQ=255)
                                                                                          tACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGGTCTTCGGa                                                                                                                                            >  1:525876/1‑239 (MQ=255)
                                                                                          tACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGGGAGTTATTGCTGGGTGACGCGTGGTTGAGGGCGGTTGGTTTAAGGTTGACAATTGCGGCCTGGGCGGGGAcggcgggggggggccccgga                                                                                                                                            >  2:562127/1‑228 (MQ=255)
                                                                                          tACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGGGAGTTATTGCGCGATGACGCGATGATGAGGGCGGTTGGTTTAAGTTTGACAATTGGGGGCTTGgcggggggggggggggggggccgggga                                                                                                                                            >  2:599993/1‑216 (MQ=255)
                                                                                             cGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCACCGCGCTGGTATGGTGAGTTATTGCGCGATGCCGCGCTGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGGGGGGCGACGGCGGCGGCGGTTCTCGGaaat                                                                                                                                         >  2:87835/1‑237 (MQ=255)
                                                                                                 tgttcggttGCCCGGGGGGTTTGTGATTTGGCTGTTGGGGTTTACCTTTCTCAACGCGCCAATTCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGACGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTg                                                                                                                                     <  2:768295/238‑1 (MQ=255)
                                                                                                     cgttcgGCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCt                                                                                                                                                <  1:744059/224‑1 (MQ=255)
                                                                                                     cgttcgCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGGGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTgcggc                                                                                                                                 >  1:67718/1‑239 (MQ=255)
                                                                                                       ttcgCCGTGGCGGATTGTGATTTTGCTGCTGGGCTTTTCCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGcggcgg                                                                                                                               <  1:501261/239‑1 (MQ=255)
                                                                                                             cTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGGGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGGGACGGCGGCGGCGATCCTCGGAACTATTGCGGCggtggggg                                                                                                                         >  1:616793/1‑236 (MQ=255)
                                                                                                               ggCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGGGTGCGGTTGGTTTAAGTCTGGCAATTGCGGCCTGGGCGGGGGCGGCGGGGGCGGGCCCCGGAAcatttctggttcttttaaac                                                                                                                       >  2:744059/1‑220 (MQ=255)
                                                                                                                   gATTGGGATTTTGCGTCTGGGTTTTCCCTTTCTCTACGCGGCAATGCTGATGCTGGTTATCTATTCGTTTAACCGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGtggtgctggtg                                                                                                                    <  1:423677/238‑1 (MQ=255)
                                                                                                                      ttgtttttttttgtggggTTTTACCTTTTCCTAGGCGCCAATGCTGATGGTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATTCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCt                                                                                                                <  2:877248/224‑1 (MQ=255)
                                                                                                                              tgctgcTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGggtggctgggggccacttgcggg                                                                                                        >  1:822868/1‑225 (MQ=255)
                                                                                                                                   tttggtttttCCTTTTTCTACGCGCCAATGGTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGAGGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTcg                                                                                                   <  2:997750/233‑1 (MQ=255)
                                                                                                                                   ctgGGCTTTCCCTTTCCTAGCGCGCCAATGTTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTcgc                                                                                                  <  1:1068914/239‑1 (MQ=255)
                                                                                                                                     tgggtttACCTTTCTCTACGGGCAAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGGGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTcgcg                                                                                                <  2:1059736/235‑1 (MQ=255)
                                                                                                                                            aCCTTTCCCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTGGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCaaa                                                                                        <  1:28584/239‑1 (MQ=255)
                                                                                                                                                  ctctACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGCGGGGCTGGGGCGCTTTGTCTGGTTTTGCGGGCCCAAAGGCtt                                                                                  >  1:128563/1‑239 (MQ=255)
                                                                                                                                                    ctACGCGCCAATGCTGATGCTGGTTATCTATTCGGTGAACAGCGTGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTg                                                                                <  1:543226/239‑1 (MQ=255)
                                                                                                                                                       tcgcgCCAATGTTGATGTTGGTTTTTTTTTCTGTTAACGGCTCGAGGCTGGTGACGGGGTGGGCCGGCTGGTCAACGCGCTGGTTTGGTGAGTTATTGCGGGATGACGGGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGcc                                                                              <  2:251276/238‑1 (MQ=255)
                                                                                                                                                        cgcgCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGGTGGGTGGCGTTGGGAGGATTTGCGGGTTCAAAGGGGTTTg                                                                                >  1:468204/1‑237 (MQ=255)
                                                                                                                                                        cgcgCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGCTGGGGGGCTTTGGCAGGTGTCTCgggttcaataggctgtg                                                                                >  1:233515/1‑222 (MQ=255)
                                                                                                                                                        cgcgCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGGGGGGCTGGTGGGCCTTTGCGGGGTTTGCGGGTCCAAAGGGttttc                                                                               >  1:251276/1‑235 (MQ=255)
                                                                                                                                                         gcgcCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGTGGCTGGGGGGCTTTGGGAGGGTTTCCCGATCAACAGGCttttgcct                                                                            >  1:412448/1‑234 (MQ=255)
                                                                                                                                                                                        ttaacgggctgagggggggtgcGGGGTGGGCCGGGTGGTCATGGCGTTGGTTTGGTGAGTTATTGCGGGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGAt                                          <  2:1063893/218‑1 (MQ=255)
                                                                                                                                                                                               cTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGGGGGCTTTTTCAGGTTTTCCGGGTCAAAATGCTTTTCCTTTTTGTTTCACCCGCCCCCGGGGATTCCAGATGttctc                                    >  1:804184/1‑234 (MQ=255)
                                                                                                                                                                                                       tGGTGACGGTGGGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCt                             <  1:562127/237‑1 (MQ=255)
                                                                                                                                                                                                        gggtgcGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGTTTGGTTTAAGTCTGACAATGGCGGCCTGTGCGGGGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGt                          <  1:87835/234‑1 (MQ=255)
                                                                                                                                                                                                         gTGACGGTGTGGGCCGGTTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGt                          <  1:786893/238‑1 (MQ=255)
                                                                                                                                                                                                             aggggtgGGCCGGCTGGTCAACGCGGTTGTCTGGTGCGTTATTGCGCGATGACGGGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCtg                     <  1:187780/235‑1 (MQ=255)
                                                                                                                                                                                                               gggtgGGGCGGCGGGGCGCCGCGGTGGTTGGGGGAGTTGTTGCGGGGTGACGGGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGGGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCtgt                    <  2:128563/236‑1 (MQ=255)
                                                                                                                                                                                                                        gggtGGTAAACGGGGTGGTTTGGGGGGTTAGTGCGGGATGACGGGATGGTGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTattatt          <  2:935548/236‑1 (MQ=255)
                                                                                                                                                                                                                          cTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGCGGGTGGGCTTTGGGCGGGTTCGCGGGTCAAAATGGCTTGGCTTTTTTGTTCCCCCGCCGCAGGGAAGGCCCAGATGTCTTCCCGGGCTTTTGGCGTTGtttttttttt        >  2:1035099/1‑237 (MQ=255)
                                                                                                                                                                                                                          cTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGGGTTGGGGGTGGGCGTTGGGTGGGTTTGCTGATCCAACAGGTTTCGCTTTTTATTACCCCCCCCCCGCGTAATGGCAGATGTCACCTCggggtgttttctgttggttttttttc        >  2:604018/1‑216 (MQ=255)
                                                                                                                                                                                                                           tGGTCACGGGGCTGGTTTGGTGGGTTATTGCGCGATGAGGCGATGAGGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATtcgt       <  2:616793/239‑1 (MQ=255)
                                                                                                                                                                                                                            gggcAACGCGTTGGTATGTGGAGTTAGTGCGCGATGACGCGATTGTGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATtcgt       <  1:955964/235‑1 (MQ=255)
                                                                                                                                                                                                                                aaCGCGCTGGTAGGGGGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCgcgc  <  1:230772/239‑1 (MQ=255)
                                                                                                                                                                                                                                     gcTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAg                                            <  1:390864/192‑1 (MQ=255)
                                                                                                                                                                                                                                       tGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCt                             <  1:212876/205‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
GCGATCATCATGTTGCTGCTGCTAATTGTGCCGATAATGTGGTTCCACAAACACCAGCAAAAAAGCGTGGGAGAACACGGATGAATAATTTACCGGTAGTTCGTTCGCCCTGGCGGATTGTGATTTTGCTGCTGGGCTTTACCTTTCTCTACGCGCCAATGCTGATGCTGGTTATCTATTCGTTTAACAGCTCGAAGCTGGTGACGGTGTGGGCCGGCTGGTCAACGCGCTGGTATGGTGAGTTATTGCGCGATGACGCGATGATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGC  >  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: