Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,514,139 A→C I57S (ATC→AGC)  fecB ← ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,514,1390AC96.8% 196.6 / ‑3.9 63I57S (ATC→AGC) fecBferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (28/33);  total (29/34)
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.87e-01

TTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGC  >  NC_000913/4513895‑4514383
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ttGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGt                                                                                                                                                                                                                                                 >  1:214520/1‑251 (MQ=255)
    aaaTTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTTAACCCCCGCATGGCGACTGCTGGCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCAc                                                                                                                                                                                                                                             <  1:439351/251‑1 (MQ=255)
    aaaTTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCAc                                                                                                                                                                                                                                             >  1:222468/1‑251 (MQ=255)
     aaTTTTCAGCGTAGGTTTCGTTGCGGGCCTTAAGCAGCATTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGTTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCAcg                                                                                                                                                                                                                                            <  1:19605/251‑1 (MQ=255)
               gTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAgcg                                                                                                                                                                                                                                  >  1:19416/1‑251 (MQ=255)
                  ggTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCACGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTACCCGGCGGcccgcccat                                                                                                                                                                                                                               >  2:118378/1‑244 (MQ=255)
                     ttCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCg                                                                                                                                                                                                                             <  1:472776/250‑1 (MQ=255)
                             gggACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACg                                                                                                                                                                                                                      >  1:211642/1‑249 (MQ=255)
                               gACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTAGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACgag                                                                                                                                                                                                                    >  1:226822/1‑249 (MQ=255)
                                  ttAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTc                                                                                                                                                                                                               >  1:277189/1‑251 (MQ=255)
                                         gcagttccgggtcgATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGcacc                                                                                                                                                                                                        <  2:403123/239‑1 (MQ=255)
                                                   cgcgATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGtt                                                                                                                                                                                              <  1:42915/251‑1 (MQ=255)
                                                    gcgATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTtg                                                                                                                                                                                             <  1:13879/251‑1 (MQ=255)
                                                      gATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTtgtg                                                                                                                                                                                           <  1:284954/251‑1 (MQ=255)
                                                      gATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTACAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTtgtg                                                                                                                                                                                           <  2:159760/251‑1 (MQ=255)
                                                          tgctgcAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGGTTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCg                                                                                                                                                                                        <  2:167762/250‑1 (MQ=255)
                                                           gctgcAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCTGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGtt                                                                                                                                                                                      <  2:38200/251‑1 (MQ=255)
                                                                   gCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCACGCGCATCGGCGGACGGGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCg                                                                                                                                                                                >  2:381937/1‑249 (MQ=255)
                                                                        gTAAACCCCCGCATGGCGACTGCTGTCGGCAATGTTCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCCCGGCGGCCAGCGCATCGGCTAACGCGAGATCCAGCAACACCATCCGTTGGGGCGTTTTTTCGAGTGTa                                                                                                                                                                          >  2:104738/1‑250 (MQ=255)
                                                                         aaaccccccGCTTGGCGACTGAGTGCGGCAATGTTCAGGGCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCAACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAAGCCGTTGTGGCGTTTTTTCGAGTGTaaa                                                                                                                                                                        <  2:230023/247‑1 (MQ=255)
                                                                                ccGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCg                                                                                                                                                                 <  1:374372/251‑1 (MQ=255)
                                                                                      ggCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCCATGGTTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTc                                                                                                                                                            <  2:269749/250‑1 (MQ=255)
                                                                                      ggCGACTGCGGTCGGCAATGGTCAGGTCTGGTTTCAGGGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGTACTGCCACGGTTTCAGGGGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGGTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAATCGTGCCGTGTTc                                                                                                                                                            <  2:336879/250‑1 (MQ=255)
                                                                                          cctgcggTCGGCAATGTGCAGGTCTGGTTGCAGAGCGTCAAGGGCTTCCATGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGGCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTGCGAGTGTAAACGTGCCGTGTTCGTCCt                                                                                                                                                       <  2:369812/250‑1 (MQ=255)
                                                                                          acttctgTCGGCAATGATCTGGTCGGGTTTCAGAGGTGCAATGGCTGCCAGGCTCGGTTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCCCTTCGGGCAGGATGCGTTTTGCAGCGTTATCGGCGGCAATACCGATCGGGCTGACGTCCACGTCTGCCAGCGCATCGGCTAACTAGAGTTCCAGCACCTCAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCt                                                                                                                                                       <  2:366894/247‑1 (MQ=255)
                                                                                               tgTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCCGCGCATCGGCGAACGAGAGGTCCAGAACCACAATCCGTTGATGCGTTTTTTCGAGAGAAAAACTGCCGGGTTCGACCTGAAcc                                                                                                                                                  >  2:262410/1‑251 (MQ=255)
                                                                                                     gCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTggc                                                                                                                                             >  1:333513/1‑251 (MQ=255)
                                                                                                          gATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCa                                                                                                                                         <  1:345264/250‑1 (MQ=255)
                                                                                                           aTCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCaa                                                                                                                                        >  2:25671/1‑250 (MQ=255)
                                                                                                             cAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAgg                                                                                                                                     >  1:312842/1‑251 (MQ=255)
                                                                                                                    ggTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTgg                                                                                                                                <  1:365888/249‑1 (MQ=255)
                                                                                                                         cAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATc                                                                                                                          >  1:328347/1‑250 (MQ=255)
                                                                                                                            ggCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACca                                                                                                                      <  1:317155/250‑1 (MQ=255)
                                                                                                                                  aaTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCATGGCTGATCACcagcagc                                                                                                                 >  1:24351/1‑250 (MQ=255)
                                                                                                                                   aTGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACcagcagcag                                                                                                               <  1:15725/251‑1 (MQ=255)
                                                                                                                                         tCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCt                                                                                                           <  1:349405/249‑1 (MQ=255)
                                                                                                                                              gCTCGGCTGCGCGCGCGTTCCGCCGGACTGCCCCTGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTTCATCGTTATCGTCGGCAATTCCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGGTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCaaaaa                                                                                                    <  2:316569/251‑1 (MQ=255)
                                                                                                                                               cTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCaaaaa                                                                                                    >  2:364948/1‑250 (MQ=255)
                                                                                                                                                 cGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCaaaaacaa                                                                                                 <  2:188682/251‑4 (MQ=255)
                                                                                                                                                   gCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACaaaa                                                                                               <  1:25671/251‑1 (MQ=255)
                                                                                                                                                      gcgcgcgcGTTCCGACGGACTACCACGGTTTCAGGTGCGCACGCCCTTCGGGCAGGATGCGGTTTGCATCGGTATCGGCGGCAATCCCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACgg                                                                                            <  2:404468/251‑1 (MQ=255)
                                                                                                                                                         ggcgcgTTCCGCCGGACTGCCACGGTTTCAGGTGCGCACGCCCTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATa                                                                                         <  2:107089/250‑1 (MQ=255)
                                                                                                                                                              gTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCACAGGCGTGGCTGATCACCAGCAGCACGCCTGCAAAAACAAAACGGATAAATGc                                                                                    >  2:326556/1‑251 (MQ=255)
                                                                                                                                                                  cGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGTATAAATGCCAAc                                                                                <  1:111608/251‑1 (MQ=255)
                                                                                                                                                                      ggACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCg                                                                                                                                                           >  1:246663/1‑172 (MQ=255)
                                                                                                                                                                      ggACTGCCACGGTTTCAGGTGCGCACGCACTGCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCTCGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTGTCGAGTGGAAACGTGCCGTGTTCg                                                                                                                                                           <  2:246663/172‑1 (MQ=255)
                                                                                                                                                                       gACTGCCACGGTTTCAGGTGCGCACGCACTTCTGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCCGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGGGTGTAAACGTGCCGTATTCGTCCTGAAACGTGCCGGCCAAGGCGTGGCTGTTCACCAGCAGCAGGCCTGCAAAAACACAACGGGTAAAAGCCAACATAAt                                                                          >  2:223017/1‑251 (MQ=255)
                                                                                                                                                                           gCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATcac                                                                       <  1:326556/250‑1 (MQ=255)
                                                                                                                                                                               cgggttCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGGGTTTTTTCGAGTGTAAACGTGCCGGGTTCGTCCTGAACCGGGGCGGCAAAGGCGAGGCTTATCACCAGCAGCACGCCTGCACAAACAAACAGGATAAATGCCaaaaaaaacccctacc                                                                   >  2:360286/5‑237 (MQ=255)
                                                                                                                                                                                  tttCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGc                                                               >  1:14090/1‑251 (MQ=255)
                                                                                                                                                                                        gTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTaaa                                                           >  1:373083/1‑249 (MQ=255)
                                                                                                                                                                                            gcacgcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCg                                                     <  2:421829/251‑1 (MQ=255)
                                                                                                                                                                                                  acTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGc                                               >  1:385302/1‑251 (MQ=255)
                                                                                                                                                                                                       gggCAGGATGCGTTTAGCATCGGTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGATCCAGCACCACAATCCGTTGTGGAGTTTTTTCGAGTGTAAACGGGCCGTGTTCGTCCTGAACCGGGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCAAGCTAAAAGCCCGGCAAGACGGGc                                          >  1:470310/1‑251 (MQ=255)
                                                                                                                                                                                                              aTGCGGTTTGCATCGTTCTCGTCGGCAATACCGATCGCGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACAGAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAaca                                   <  2:4252/251‑1 (MQ=255)
                                                                                                                                                                                                                 cGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCAAAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGACGGCCTGCAAAAACCAAACGGAAAACTGCCACACTAATCACATTCCAGCTAAAAACCCGGCAAGCCGGGCGTTAACACAt                                >  2:392923/1‑251 (MQ=255)
                                                                                                                                                                                                                  gTTTGGCATCGTGATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTGCAGCACCACAATCCGTTGTTGCGGTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGTCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGGCAACATATTCTCATTCCAGCTAAAAGCCCGGCAAGCCGTGCGTTAACACATc                               <  1:277983/251‑1 (MQ=255)
                                                                                                                                                                                                                    tttGCATCGTTATCGTCGGCAATACCGATCGGGCTGCCGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGGTCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACAAAATACCCTTCCAGCTAAAAGCCCGACAAGCCGGGCGTTAAAACATc                               >  2:170758/1‑249 (MQ=255)
                                                                                                                                                                                                                         aTCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAACGGCGTGGCTGATCACCAGCATCAGGCCTGCAAAAACAAAACGGATAAATGCCAAAATAATCACATTCCAACTAAAAGCCCGGCAAGCCGGGCGTTAACACATCCAAACTTc                       >  2:73961/1‑251 (MQ=255)
                                                                                                                                                                                                                            gTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCa                      >  2:39436/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aaTACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGcct         >  2:255349/1‑249 (MQ=255)
                                                                                                                                                                                                                                            ccGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATAc     <  1:255349/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aTCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGc  <  2:328347/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
TTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGC  >  NC_000913/4513895‑4514383

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