Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 624,319 A→C D64E (GAT→GAG fepB ← ferrienterobactin ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913624,3190AC75.8% 119.0 / 47.6 66D64E (GAT→GAGfepBferrienterobactin ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (7/9);  new base C (20/30);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GGGGCGATGGTGGAAAGCTGATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATT  >  NC_000913/624089‑624529
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                  
ggggCGATGGTGGAAAGCTGATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCg                                                                                                                                                                                                             >  2:143762/1‑238 (MQ=255)
ggggCAATGGTGGAAAGCTGAACACACGGGGCCGGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGc                                                                                                                                                                                                            <  2:1010147/239‑1 (MQ=255)
       tggtggAAAGCTGATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATTACCGGAGCCTCAATCGCCAGCAGt                                                                                                                                                                                                     >  2:36846/1‑239 (MQ=255)
               cgCTGATCATACAGTGCCAGCGCCGAACCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGt                                                                                                                                                                                             <  1:158532/238‑1 (MQ=255)
               aGCTGATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGt                                                                                                                                                                                             <  2:761560/239‑1 (MQ=255)
                           aGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGAc                                                                                                                                                                                    >  1:103427/1‑236 (MQ=255)
                             tGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCtgg                                                                                                                                                                               >  2:50878/1‑239 (MQ=255)
                              gCCAGCGCGCAATCCCCGCCGGTTGCGCTAATTAAAAGCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCCGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCtggt                                                                                                                                                                              <  1:960427/239‑1 (MQ=255)
                              gCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCtggt                                                                                                                                                                              >  1:820634/1‑239 (MQ=255)
                               ccAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGt                                                                                                                                                                                             >  2:934923/1‑223 (MQ=255)
                               ccAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGt                                                                                                                                                                                             <  1:934923/223‑1 (MQ=255)
                                  gcgcCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACg                                                                                                                                                                                   >  2:345156/1‑230 (MQ=255)
                                  gcgcCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGaa                                                                                                                                                                          >  2:469264/1‑239 (MQ=255)
                                  gccccGAACCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACg                                                                                                                                                                                   <  1:345156/227‑1 (MQ=255)
                                   cgcCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGaa                                                                                                                                                                          >  2:346783/1‑238 (MQ=255)
                                   cgcCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGaa                                                                                                                                                                          >  2:193574/1‑238 (MQ=255)
                                    cccgacaccccGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCGGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCtggtgg                                                                                                                                                                            <  1:329919/228‑1 (MQ=255)
                                        aaTCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCt                                                                                                                                                                                 <  2:1128201/226‑1 (MQ=255)
                                        aaTCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCt                                                                                                                                                                                 >  1:1128201/1‑226 (MQ=255)
                                                  ggTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCg                                                                                                                                                             >  2:725884/1‑236 (MQ=255)
                                                     tGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCggctgg                                                                                                                                                        <  1:143762/238‑1 (MQ=255)
                                                        gcTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCtt                                                                                                                                                     <  1:469264/238‑1 (MQ=255)
                                                          tAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTc                                                                                                                                                   <  2:103427/238‑1 (MQ=255)
                                                          tAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTc                                                                                                                                                   <  1:342853/238‑1 (MQ=255)
                                                          tAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTc                                                                                                                                                   <  1:197891/238‑1 (MQ=255)
                                                           aaTTAAAATCAGATCTGGCATTTGCGCGGCAACGGCTTCCGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCa                                                                                                                                                 <  2:226796/239‑1 (MQ=255)
                                                           aaTTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCa                                                                                                                                                 <  1:36846/239‑1 (MQ=255)
                                                            aTTAAACTCGGATCAGACATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAg                                                                                                                                                <  2:133267/239‑1 (MQ=255)
                                                                 aaTGAGTTCTGGTATTTGGGCGCCAACGGGTACGGCGTTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCCCCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGGGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTAt                                                                                                                                           <  2:1015272/239‑1 (MQ=255)
                                                                                 tGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACg                                                                                                                                 >  1:532822/1‑233 (MQ=255)
                                                                                                tCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAg                                                                                                                                                                                                  <  1:825397/153‑1 (MQ=255)
                                                                                                tCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAg                                                                                                                                                                                                  >  2:825397/1‑153 (MQ=255)
                                                                                                        cGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCa                                                                                                     >  1:517346/1‑238 (MQ=255)
                                                                                                            tCGCCGATATAGGGCCGGTGCAGTTTGCGTTCTTTCGCCCCCTTGCTCCACTGGCTTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGGGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGc                                                                                                <  2:874009/239‑1 (MQ=255)
                                                                                                             cGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGcc                                                                                               <  2:852237/239‑1 (MQ=255)
                                                                                                                   tataGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGa                                                                                            >  1:928571/1‑236 (MQ=255)
                                                                                                                           cGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCgg                                                                                     <  1:973577/235‑1 (MQ=255)
                                                                                                                            gTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTc                                                                                                                                                   >  1:947368/1‑172 (MQ=255)
                                                                                                                            gTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTc                                                                                                                                                   <  2:947368/172‑1 (MQ=255)
                                                                                                                             ttGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGc                                                                                                                               >  2:116883/1‑191 (MQ=255)
                                                                                                                             ttGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGc                                                                                                                               <  1:116883/191‑1 (MQ=255)
                                                                                                                             ttGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTAtt                                                                                >  1:86733/1‑238 (MQ=255)
                                                                                                                             ttGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACa                                                                                                                                                                       <  2:129747/151‑1 (MQ=255)
                                                                                                                             ttGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACa                                                                                                                                                                       >  1:129747/1‑151 (MQ=255)
                                                                                                                             ttGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTAtt                                                                                >  2:640679/1‑238 (MQ=255)
                                                                                                                                          tCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCt                                                                                                                              <  2:94636/179‑1 (MQ=255)
                                                                                                                                          tCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCt                                                                                                                              >  1:94636/1‑179 (MQ=255)
                                                                                                                                            tttCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTc                                                                <  1:870015/239‑1 (MQ=255)
                                                                                                                                            ttgcGCCACCTTTCTCCACTGGCGGAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGGGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTc                                                                <  1:785785/236‑1 (MQ=255)
                                                                                                                                              tctcagccTTTCACCAGTGCCGTAGAAAGCACGCGTCATCCGCTGAGCGGATATTCGGCGGGGTCGCGCCGCTGGCGATCAACGAAGCCTCAAACGCCAGCAGTGAGCGGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCt                                                              <  2:1050171/233‑1 (MQ=255)
                                                                                                                                                        gCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACg                                                                                                                                                                                   <  1:710805/112‑1 (MQ=255)
                                                                                                                                                        gCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACg                                                                                                                                                                                   >  2:710805/1‑112 (MQ=255)
                                                                                                                                                           ccACTGGCGGAAAAAGCCCTGGTCATCCGCGACGCGGTTCTTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAggg                                                 <  1:198543/239‑1 (MQ=255)
                                                                                                                                                             aCTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCg                                               <  1:50878/239‑1 (MQ=255)
                                                                                                                                                                  cGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGt                                              <  2:86733/235‑1 (MQ=255)
                                                                                                                                                                   gTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCCATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCgg                                         >  2:846549/1‑239 (MQ=255)
                                                                                                                                                                            ccTGGTCCTCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTagag                                    <  1:300765/235‑1 (MQ=255)
                                                                                                                                                                                  cATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAg                           <  1:707698/238‑1 (MQ=255)
                                                                                                                                                                                           cgcgGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAAt                 <  1:640679/239‑1 (MQ=255)
                                                                                                                                                                                                ttATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGc             <  1:725884/238‑1 (MQ=255)
                                                                                                                                                                                                        cGTGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGtt     <  1:280202/238‑1 (MQ=255)
                                                                                                                                                                                                          tGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCtggt                                                                                                                                                                              >  2:829990/1‑67 (MQ=255)
                                                                                                                                                                                                          tGGTCGCGCCGCTGGCGATCACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCtggt                                                                                                                                                                              <  1:829990/67‑1 (MQ=255)
                                                                                                                                                                                                          tGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTCTTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGttatt  >  1:480592/1‑239 (MQ=255)
                                                                                                                                                                                                                             cACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCg                                               <  1:1068645/175‑1 (MQ=255)
                                                                                                                                                                                                                             cACCGGAGCCTCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCg                                               >  2:1068645/1‑175 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                  
GGGGCGATGGTGGAAAGCTGATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCCGATATAGAGCCGTTGCAGTTTGCGTTCTTTCGCCACCTTGCTCCACTGGCGTAAAAAGCCCTGGTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATT  >  NC_000913/624089‑624529

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