Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,513,044 T→C S122S (TCA→TCG fecC ← ferric citrate ABC transporter membrane subunit FecC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,513,0440TC95.2% 116.3 / ‑2.5 43S122S (TCA→TCGfecCferric citrate ABC transporter membrane subunit FecC
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (20/20);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

GTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCG  >  NC_000913/4512793‑4513287
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
gTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCGCATGCCGCGatgaat                                                                                                                                                                                                                                                     >  1:371980/1‑251 (MQ=255)
                   gCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCa                                                                                                                                                                                                                                   <  1:75291/251‑1 (MQ=255)
                             aaccaaccggCGGGCGTGGGAACCTCGTCCTGCCAGCCAGTAAAAGAGGCCGTAAGCAGGATCTGCGGCCAGCAGGAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACATCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCg                                                                                                                                                                                                                         <  2:212621/242‑1 (MQ=255)
                               ccTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCAAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCgg                                                                                                                                                                                                                        >  1:351996/1‑250 (MQ=255)
                                cTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGAc                                                                                                                                                                                                                      >  1:118696/1‑251 (MQ=255)
                                            gTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGACCGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTgg                                                                                                                                                                                                          >  1:279598/1‑251 (MQ=255)
                                                gACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACGCAGCGCGCTGGTAAg                                                                                                                                                                                                      <  1:296352/251‑1 (MQ=255)
                                                       cGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCGTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATa                                                                                                                                                                                               <  1:345738/251‑1 (MQ=255)
                                                                     tAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAgcgccgcgcc                                                                                                                                                                                 >  1:98441/1‑251 (MQ=255)
                                                                      aaaaGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGcgccgcgccg                                                                                                                                                                                >  1:228164/1‑251 (MQ=255)
                                                                          gATGCCGTAAGCATGACTTCCGGCCAGTAGCAGGGGGATGCGGGTAATGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTGCGGTCAGGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGt                                                                                                                                                                            <  2:98441/251‑1 (MQ=255)
                                                                               cGTAAGCATGATCTTCGGACAGCAGCAGGGTGATACGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGCGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATg                                                                                                                                                                       >  1:89411/1‑251 (MQ=255)
                                                                                  aaGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCg                                                                                                                                                                    >  1:219133/1‑251 (MQ=255)
                                                                                   aGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCg                                                                                                                                                                    >  1:275339/1‑250 (MQ=255)
                                                                                     cATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGGAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGc                                                                                                                                                                 >  1:165567/1‑251 (MQ=255)
                                                                                              cGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGTGCAGTGAAgg                                                                                                                                                         >  1:88139/1‑250 (MQ=255)
                                                                                                  cgagccgcaggcggATGGGGGTAAGGCCAATGCAAGAGGCCGGAAGCGCGATACGAGCGAGGGTTAGGTTGTGTCTGGCCGGTGGATGACGAAATCCGCCTCCGGCGGTCATGACCAGCATCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACGCAGCGCGCGTGTAAGCGCCATGGCCTGGGCGGCGCCGCTGTTAATGCCGGGTAGTGAAGGagag                                                                                                                                                     <  2:175040/238‑1 (MQ=255)
                                                                                                   cagcagcagGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAgg                                                                                                                                                    >  1:285929/1‑251 (MQ=255)
                                                                                                      cagcagGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGcc                                                                                                                                                  >  1:69710/1‑250 (MQ=255)
                                                                                                        gcagGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCAtt                                                                                                                                               <  1:181675/251‑1 (MQ=255)
                                                                                                        gcagGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCAtt                                                                                                                                               <  1:376200/251‑1 (MQ=255)
                                                                                                         ccgtgtaaTGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACTCGCGAGGCTCAGTTTGTTTCTGTCATGGGTATGACGAAGTCCGCCTCATGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGCTGAACGACAGGGAATAGCCTGAAATAGGCGTCGGACTACGCGTGCGGGTCAGCGACATGGCTAGCGCCGCGCGGCGGTTAATGCCTGGCAGTGACGGGGAgggaagtg                                                                                                                                              >  2:136111/8‑245 (MQ=255)
                                                                                                                  gCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTggg                                                                                                                                     <  1:143476/251‑1 (MQ=255)
                                                                                                                          ggCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGtt                                                                                                                             <  1:165376/251‑1 (MQ=255)
                                                                                                                            cccATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTg                                                                                                                           >  1:214674/1‑251 (MQ=255)
                                                                                                                              cATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCACTCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGACATTGGGTTGTGGGTCAGGGTTTgcc                                                                                                                         >  2:62079/1‑250 (MQ=255)
                                                                                                                                tACAAAAGGCAGACAGCGCGATACCCGCGAGTGTAAGTTTGATTCAGTCATGGGTCTGACGAAATCCGCCTCCGGGGGTCATGACCAGCAGACAGCTCACGCCACCCCCGCAGGCCGCGATGACAGACAGCGAATAGCCTGCAATCGGCGTCGGACTCAGCGAGCTGGGAAGCGACAGAGCCAGCGCCGCTCCGCAATTAATGCCACGCAGGGCAGGAGAGGCACTCGGGATGTGGGACAGGGATTgccgc                                                                                                                       >  2:42315/1‑248 (MQ=255)
                                                                                                                                                cgcgATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAAt                                                                                                                                                                                                                                 >  1:12481/1‑129 (MQ=255)
                                                                                                                                                cgcgATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAAt                                                                                                                                                                                                                                 <  2:12481/129‑1 (MQ=255)
                                                                                                                                                                  cAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATc                                                                                      >  1:8644/1‑250 (MQ=255)
                                                                                                                                                                                   ggggaTGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCtt                                                                    <  2:275339/247‑1 (MQ=255)
                                                                                                                                                                                   gggTATTACGAAAACCGCATCATGCGGTTATGACCATCCGCCTGCTCACGCCGCCCCCGCATGCCGCTATGAACGACAGAGAATAGCCTGCAAGCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATATCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCtt                                                                    <  2:214674/251‑1 (MQ=255)
                                                                                                                                                                                      gttagtgatatccccttccttgggTCATGGCCTGCAGCCAGTCCACGCCGACCTCGCATTCCGCGATGAACGACAGAGAATAGCCTGCACTCTGCGTCGGTCTCAGCGCGCGTGTAAGCGCCATAGCCTGTGCCTCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGTGTCAGGGTTTGCAGCAGCTTGCCCGCGAGCGCCTGGCTTGCGCCGATCAGTTCGTCGACCAGGCTTCg                                                                  <  2:119611/228‑1 (MQ=255)
                                                                                                                                                                                       aTGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGCAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTg                                                                >  1:357031/1‑251 (MQ=255)
                                                                                                                                                                                                              tCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAgc                                         >  1:83652/1‑251 (MQ=255)
                                                                                                                                                                                                               cctGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGCGCGGCGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGGGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAACCAAAGGGTTTGCACCAgcg                                        >  2:270255/3‑251 (MQ=255)
                                                                                                                                                                                                               cATGCCCAGCAGCCACCTCACGCCGCCCCCGCGTGCCGCGATGAACGACAGAGAAGAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGTGCCGCGCCGCTGTGAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTTGCAAACGAAGGTTTTGCACCAgc                                         <  2:226823/250‑1 (MQ=255)
                                                                                                                                                                                                                  gACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCt                                      <  1:17317/250‑1 (MQ=255)
                                                                                                                                                                                                                           ccagcTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCg                            <  1:118434/251‑1 (MQ=255)
                                                                                                                                                                                                                             agcTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGtt                          <  1:6704/251‑1 (MQ=255)
                                                                                                                                                                                                                              gcTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTg                         <  1:57423/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ccGCATGCCGCGATGATCGACAGAGAATAGCCGGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCGGCGTCGCTGATAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGTAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGTCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAgg            <  1:242850/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ccGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAgg            <  1:62079/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     cgATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGAAGGgcg  >  1:99956/1‑251 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
GTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCG  >  NC_000913/4512793‑4513287

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