Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 248,203 G→A intergenic (+69/‑2695) rayT → / → dinB REP‑associated tyrosine transposase/DNA polymerase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913248,2030GA97.7% 148.0 / NA 45intergenic (+69/‑2695)rayT/dinBREP‑associated tyrosine transposase/DNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (25/18);  minor base C (1/0);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTG  >  NC_000913/247962‑248446
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
ataGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATa                                                                                                                                                                                                                                             >  1:21017/1‑250 (MQ=255)
ataGGCAACATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGTACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTAATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTATTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATa                                                                                                                                                                                                                                             >  2:46670/1‑250 (MQ=255)
  aGGCATCATGTTGATTATATTTATATAAATCCAGTAAGGCATGGTTGGGTAAAGCAAGTGAGTGATCGGCCATTCTCAACGTACCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGAAGCGACGCTAGAGCGTCTTCTCATGCCTACAAACTTGTGACGGATCGGTAAGCCGGATCAgg                                                                                                                                                                                                                                          >  1:104496/1‑251 (MQ=255)
   ggCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGCTTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGAGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATAGGTAAGCCGGATAAggg                                                                                                                                                                                                                                         >  2:241472/1‑250 (MQ=255)
          tGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATg                                                                                                                                                                                                                                  >  1:352651/1‑251 (MQ=255)
                    atttatatAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGCTCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCgg                                                                                                                                                                                                                        >  1:169067/1‑251 (MQ=255)
                    atttatatAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCg                                                                                                                                                                                                                         <  1:343278/250‑1 (MQ=255)
                       gatataAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCa                                                                                                                                                                                                                      <  1:291539/249‑1 (MQ=255)
                             aaTCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATgc                                                                                                                                                                                                               <  2:210335/251‑1 (MQ=255)
                              aTCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATgcg                                                                                                                                                                                                              >  1:234662/1‑251 (MQ=255)
                                  agcaaaGCATGTCTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCTATGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTGATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGa                                                                                                                                                                                                          <  2:139078/248‑1 (MQ=255)
                                    tAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATg                                                                                                                                                                                                        >  1:134263/1‑251 (MQ=255)
                                                gggTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGGAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGAc                                                                                                                                                                                             <  1:45469/250‑1 (MQ=255)
                                                     aaGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTTGGCGCAGATGCCTGATGCGACGCTGgc                                                                                                                                                                                       >  2:287265/1‑251 (MQ=255)
                                                            tgagtgaTTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCtt                                                                                                                                                                                >  2:97708/1‑251 (MQ=255)
                                                                 gaTTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGCTAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATcct                                                                                                                                                                           >  2:48157/1‑249 (MQ=255)
                                                                                   ttCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGAAGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAAACTTTTGCCt                                                                                                                                                         >  2:150252/1‑251 (MQ=255)
                                                                                    tCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTg                                                                                                                                                        >  1:306494/1‑251 (MQ=255)
                                                                                       aTCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCTTGTGCCTGAAc                                                                                                                                                     >  1:330729/1‑251 (MQ=255)
                                                                                             aTGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGCTTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATACCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAgg                                                                                                                                               >  1:167770/1‑251 (MQ=255)
                                                                                                     cgAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCTTGTGCCTGAACGGTAGGACGGATa                                                                                                                                        >  2:11792/1‑250 (MQ=255)
                                                                                                      gAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATa                                                                                                                                        >  1:64310/1‑249 (MQ=255)
                                                                                                             atatCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCa                                                                                                                               <  1:241472/251‑1 (MQ=255)
                                                                                                                       gATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATcc                                                                                                                     <  1:11792/251‑1 (MQ=255)
                                                                                                                        aTTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATg                                                                                                                                                                                                                                  >  1:377098/1‑141 (MQ=255)
                                                                                                                        aTTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATg                                                                                                                                                                                                                                  <  2:377098/141‑1 (MQ=255)
                                                                                                                              cGGGGGACGCTACGGATTTTAGTGCCGTGGAGCGTATAATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAATCTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGtt                                                                                                              <  2:64310/251‑1 (MQ=255)
                                                                                                                                  ggACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATg                                                                                                                                                                                                                <  2:314248/149‑1 (MQ=255)
                                                                                                                                  ggACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATg                                                                                                                                                                                                                >  1:314248/1‑149 (MQ=255)
                                                                                                                                         aCGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTACGTCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGATGGGGCGTCTTATCATGCCTACCAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCAACCGGCAGTTGTGCCCCaagg                                                                                                   >  2:83374/1‑249 (MQ=38)
                                                                                                                                            gATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTg                                                                                                <  1:150252/251‑1 (MQ=255)
                                                                                                                                               tttAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGa                                                                                              <  1:84987/250‑1 (MQ=255)
                                                                                                                                               tttAGTGCCGGGGAGCGTATCATTCCATAATTGTGCGCAGATGCCTGATGCGACTCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGAGCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGAGGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGGCACCAATGCTGGa                                                                                              <  1:277919/250‑1 (MQ=255)
                                                                                                                                                          ggAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCa                                                                                   >  1:186980/1‑250 (MQ=255)
                                                                                                                                                                tATCATTTCATAATTGTGCGCAGATGCCTGAGGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGt                                                                                                               <  1:270640/216‑1 (MQ=255)
                                                                                                                                                                        cATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGGACGGTAGGACGGATGAGGCGTTCACGCCGCATCCGGCAGTTGTGCCCCAATGCTGGATACTTACATCAACGCCCCCCTAATCt                                                                    >  2:29651/1‑251 (MQ=255)
                                                                                                                                                                           aaTTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTTTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTc                                                                  <  1:83374/250‑1 (MQ=255)
                                                                                                                                                                                tgCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTcc                                                            >  1:368983/1‑251 (MQ=255)
                                                                                                                                                                                                gCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACtt                                            <  2:134263/251‑1 (MQ=255)
                                                                                                                                                                                                 cGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGGTAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTc                                           >  1:41990/1‑251 (MQ=255)
                                                                                                                                                                                                        aGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATAc                                    >  1:176048/1‑251 (MQ=255)
                                                                                                                                                                                                          cgcgTCGTATCATGCCTACAAACTTGTGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGa                                  <  1:305549/251‑1 (MQ=255)
                                                                                                                                                                                                                   tCATGCCTACAAACTTGGGCCGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTCTCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCCTCCGGCAGTTGTGCACCCTTGCTGGATACATACATCCACGCCCCCCTCATCTTCCACTCCACCTCATCCGGCACCTCGTTATACGACAGccaatg                         >  2:217274/1‑246 (MQ=255)
                                                                                                                                                                                                                                      ccGGATCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAaacaac      <  1:29651/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aTCGGTAAGCCGGATAAGGCGTTTATGCCGCATCCGTCAGTTATGCGCAGATGCCGTATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATATGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTg  <  2:241796/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
ATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTG  >  NC_000913/247962‑248446

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