Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,516,378 G→A N100N (AAC→AAT fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,516,3780GA96.9% 176.5 / NA 64N100N (AAC→AATfecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (30/32);  minor base C (2/0);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCT  >  NC_000913/4516137‑4516620
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
cAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGGACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTccc                                                                                                                                                                                                                                           <  2:126135/251‑1 (MQ=255)
    gTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGGTGTCCGGCGCGCTGACGCCAGGGATGCGGTCAAGTACCTCACGCAGGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGGTCAAATACGTCGGTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCGTTTGGTGCGGGCGCGGGCTCCAGCTTCCAGCTATTATTTCCCAGCg                                                                                                                                                                                                                                       <  2:81875/251‑1 (MQ=255)
       cATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCATTCGCCTACCACGTTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGt                                                                                                                                                                                                                                     <  2:182077/250‑1 (MQ=255)
           gCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTAc                                                                                                                                                                                                                                >  1:135848/1‑251 (MQ=255)
            ccAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTAcc                                                                                                                                                                                                                               >  1:7281/1‑251 (MQ=255)
             cAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCATCGGTTTTACCt                                                                                                                                                                                                                              >  1:97646/1‑251 (MQ=255)
             cAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCt                                                                                                                                                                                                                              >  1:50879/1‑251 (MQ=255)
               ggTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCCCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGGCCAGCTATTATTTCCCAGCGGTTTTACCTGc                                                                                                                                                                                                                            >  2:196688/1‑251 (MQ=255)
                gTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCCGCGGTTTTACCTGc                                                                                                                                                                                                                            >  1:237399/1‑250 (MQ=255)
                         gCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCCCGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCt                                                                                                                                                                                                                    <  2:59491/249‑1 (MQ=255)
                         gCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCCGCATGTTCCAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCCGGCTCCAGCGTCCAGCTCTTATTTCCCCGCGGTTTTACCTGCAGGCCGCTGc                                                                                                                                                                                                                  >  2:218622/1‑251 (MQ=255)
                                ccgctgttTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTTCCTCACGCATGTGTGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCCCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGTTCAGGCCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTCTTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCcag                                                                                                                                                                                                           <  2:52172/247‑1 (MQ=255)
                                  gttgttTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCcagc                                                                                                                                                                                                          >  1:67238/1‑250 (MQ=255)
                                    tgttTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATAAGTCGTTTTCACGCGCATCACCCCCCCAGCCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTACGGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCCGCAgt                                                                                                                                                                                                       >  1:45812/1‑251 (MQ=255)
                                       tttCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAgttgt                                                                                                                                                                                                    >  1:196731/1‑251 (MQ=255)
                                        ttCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCCGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAgttgtt                                                                                                                                                                                                   >  1:175091/1‑251 (MQ=255)
                                         tCCGGCGCGCTGACGCCAGGGGTGAGATTAAGTACCTCCCGAATGGTGGTTGCGCCGGCTTTGGCGAACTCCTCCCGGCGGCTCCCGTCACGCGCGACAGGATGTTCAAATACGTCGTTTTCACGCGCATCACCCCACCAGCCGCCCACCACGGTCAGGGCATCTTTTTTTGGTGCAGGCGCCGGCACCCGCGTCCAGCTATTATTTCCCAGTGGTTTTACCTGCAGTCCGCTGCCGTCCCGCAGttgttg                                                                                                                                                                                                  >  2:249917/1‑251 (MQ=255)
                                              ggcgcTGACGATAGGGATGCGGTTAAGTCCCTCACGCATGGTGGTTGCGCCGGTGTTGGCGAAATCCTCACGGCGGATCTCGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTGCACGCGCATCACCCAGCCAGTCGCCGACCTCGGTCAGGGCATCGTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGc                                                                                                                                                                                             <  1:267874/250‑1 (MQ=255)
                                                       gCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCgc                                                                                                                                                                                    >  2:121257/1‑250 (MQ=255)
                                                          aGGGATGAGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTTACGCGCGCCAGCCTGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTCTTATTTCCCAGCGGTTTTTCCTGCAGTCCGATGCCGTCCAGCTGTTGTTGCAGGCCGCTCTCGAcgc                                                                                                                                                                                 >  2:36619/1‑250 (MQ=255)
                                                           gggATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGAcgtc                                                                                                                                                                                >  1:112456/1‑251 (MQ=255)
                                                           gggATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGGGCGGGCGCGGGCTCCAGCGTCCAGCTATTAGTTCCCAGCGGTTTTACCTGCAGTCCGCTTCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGAcgtc                                                                                                                                                                                <  2:230315/251‑1 (MQ=255)
                                                             gATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGCTCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCCTCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCACCCGTTGTTGCAGGCCGCTATCGAcgtcgc                                                                                                                                                                              >  2:224069/1‑250 (MQ=255)
                                                                gCGGTTAAGTCCCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAAt                                                                                                                                                                           <  2:119413/251‑1 (MQ=255)
                                                                gCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCGCGTCACGCGCGCCAGCATGTTCAAATACGGCGGTTTCACGCGCATCACCCAGCCAGGCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCTGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCTGTTTTACCTGCAGTCCTCTGCCGTCCATCAGGTTTTGCAGTCCGCTCTCGACGTCGTAAt                                                                                                                                                                           <  2:269202/251‑1 (MQ=255)
                                                                      aaggcccTCACGCATGGTGGTCGCGCCGGTTTTGGCGAAATCCGCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATTCGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGt                                                                                                                                                                     <  2:150218/246‑1 (MQ=255)
                                                                          aCCTCACGCTGGGTGGTTGGGCCGGGTGTGGGGAGATGCGTACGGCGGATCTCGTCACGCGCGCCAGCATGGTAAAATGGGTCGGTTTCCCGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCGTCTTTTGGTGCGGGGGCGGTCTCCAGCGTCCAGCTATTATTTCCCAGCGGGTTTACCTGCAGTCCGCTGCCGTCGAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAg                                                                                                                                                                 <  2:19317/251‑1 (MQ=255)
                                                                               ccGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGt                                                                                                                                                            <  2:225998/250‑1 (MQ=255)
                                                                                 gCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTAGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCAAGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGACGTGCAGGCCGTTg                                                                                                                                                          >  1:155715/1‑251 (MQ=255)
                                                                                  cATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCTCCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGc                                                                                                                                                         <  1:128303/251‑1 (MQ=255)
                                                                                            tGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTg                                                                                                                                                 <  1:26555/250‑1 (MQ=255)
                                                                                                                 cgcACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCGCTGCTTGCCGCGCGTCAGGCTGGCGTCaa                                                                                                                           <  2:228038/249‑1 (MQ=255)
                                                                                                                   cgcggcggATCACGTCACGCGCGCAAGCATGTGCAAATACGTCGTTTTCACGCGCCTCCCCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCGGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCGCTCGACGTCGGAATCGCCGTGCAGGCCGTTGCGCTGCTTGCCGCGCGTCAGGCTGGCGTCaa                                                                                                                           <  2:175091/247‑1 (MQ=255)
                                                                                                                   cACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGGTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCCGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAAcc                                                                                                                         >  2:225392/1‑251 (MQ=255)
                                                                                                                             cacgtcacgCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTa                                                                                                                >  2:163532/1‑250 (MQ=255)
                                                                                                                               cgtcacgCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCCGGCTGGCGTCAACCGCGAGGGTAAAt                                                                                                             >  1:122454/1‑251 (MQ=255)
                                                                                                                                gtcacgCGCGCCAGCATGTTCAAATGCGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCTACCACGGTCAGGGCATCTTCTTTTTGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCTGTTGTTGCGGGCCGCTCTCGACTTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATc                                                                                                            <  1:127192/251‑1 (MQ=255)
                                                                                                                                        cgcCAGCATGTTCAAATACGTCGTTTTCACGAGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCCGAGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCCGGCTGGCGTCACCCGCGAGGGTAAATCCGCtgtgt                                                                                                    >  2:239175/1‑251 (MQ=255)
                                                                                                                                                tGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTAGCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTACTCGCCGCGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCAt                                                                                             >  1:107237/1‑250 (MQ=255)
                                                                                                                                                  ttCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGAGCGGCATACt                                                                                          >  1:59181/1‑251 (MQ=255)
                                                                                                                                                  ttCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACAACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCCGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCAGCTCTAGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGCCTGGCGGCACACGAGAGGGTAAATCTGTTGTGTGCGGCATACt                                                                                          >  1:156807/1‑251 (MQ=255)
                                                                                                                                                   tCAAATACGTCGTTTTCCCGCGCCTCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGGGCGGGCGCGGGCTTCAGCGTCCAGCTAGTATTTCCCAGCGGTGTTACCTGCAGTCCGCGGCCGGCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAAGCGCCGTGCAGGCCGTTGCTCTGCTTTCCGCGCGTCAGGCTGGCGTTTACCGAGAGGGTAAATCCGCTGTTTGCGGCATACt                                                                                          <  2:175828/250‑1 (MQ=255)
                                                                                                                                                     aaaTACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGAt                                                                                       <  1:225926/251‑1 (MQ=255)
                                                                                                                                                                 ttCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTAt                                                                           >  2:200877/1‑251 (MQ=255)
                                                                                                                                                                  tCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATc                                                                          >  1:33141/1‑251 (MQ=255)
                                                                                                                                                                   cccgcgCAGCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCTTTGCTCTGCTTGCCGCGCGGCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCGGTGTGCGGCAGACTGATTGAGCGCTTTATc                                                                          <  2:29168/248‑1 (MQ=255)
                                                                                                                                                                     cgaccATCACCCAGCCAGGCGCCTACCACGGTCAGGGCATCGTCTGTTGGTGCGGGCGCGGGCTCCAGGGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTGGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAg                                                                       <  1:12499/247‑1 (MQ=255)
                                                                                                                                                                      ccgcATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGc                                                                      <  2:8485/250‑1 (MQ=255)
                                                                                                                                                                          aTCACCCAGCCAGTCGCAGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTTCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATTCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATc                                                                  <  1:84358/251‑1 (MQ=255)
                                                                                                                                                                               ccagccagTCGCCGACCACGGTCAGGGCATCTTATTTTGGTGCGGGAGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTCCCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGCGCGATCCCGGt                                                             >  1:255139/1‑251 (MQ=255)
                                                                                                                                                                                 agccagTCGCCGACCACGGGCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGGCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGc                                                           <  1:248898/251‑1 (MQ=255)
                                                                                                                                                                                  gccagTCGCCGTCCACGGTCAGGGCATCTTCTTTTGGTGCGGGAGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCACGGTGCg                                                          <  1:148858/251‑1 (MQ=255)
                                                                                                                                                                                        cGCCGACCACGGTCATGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCGGGCGTCATCCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTa                                                    <  1:238624/251‑1 (MQ=255)
                                                                                                                                                                                               cACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCGTTATCGAGCGATCCCGGTGCGATATTAACCtgtg                                             <  1:83798/251‑1 (MQ=255)
                                                                                                                                                                                                   gTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGGAAATCCTCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGGGCGATATTAACCTGTgcagc                                         >  2:221259/1‑251 (MQ=255)
                                                                                                                                                                                                                cttTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCTCGCGATCCCGGTGAGAGATTAACCTGTGCAGCAAAAGCGGaaaa                             >  2:78578/1‑250 (MQ=255)
                                                                                                                                                                                                                  tttGGTGCGGGCGCGGGCTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCAGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGCGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACg                           >  2:46894/1‑250 (MQ=255)
                                                                                                                                                                                                                        gcgggcgcgggcCCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGAcc                    <  1:147226/251‑1 (MQ=255)
                                                                                                                                                                                                                          gggcgcgggcTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGCCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCgg                  <  2:173825/251‑1 (MQ=255)
                                                                                                                                                                                                                             cgcgggcTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCa               <  2:196731/251‑1 (MQ=255)
                                                                                                                                                                                                                             cgcgggaTCCAGAGACCGGCTATTCTTTCCCAGCGGTTTTCCCTGCCGTCAGCAGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCACCGTGCCGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACAGCGAGGGTAAATCCGCTGTGGTCGGAGTACTGATTGCGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGGGAAGCAAAAGCGGAAAAAGAGAGACCGGccc               >  2:236885/1‑250 (MQ=255)
                                                                                                                                                                                                                                 ggcTCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCgg           <  1:95311/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tCCAGCGTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGgc          >  1:5289/1‑249 (MQ=255)
                                                                                                                                                                                                                                          gTCCAGCTATTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCt  <  1:189398/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
CAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCT  >  NC_000913/4516137‑4516620

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