Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,517,556 T→G 6.6% E54D (GAA→GAC) ‡ fecR ← anti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,517,5560TG 6.6% 362.7 / 13.2 121E54D (GAA→GAC) ‡fecRanti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator
Reads supporting (aligned to +/- strand):  ref base T (62/51);  new base G (3/5);  total (65/56)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

Only 100 of 123 total aligned reads displayed.

CCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGG  >  NC_000913/4517407‑4517701
                                                                                                                                                     |                                                                                                                                                 
ccGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGttt                                                                                                                                                   <  1:682800/150‑1 (MQ=255)
 cGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGtttt                                                                                                                                                  <  2:464417/150‑1 (MQ=255)
  ggTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTc                                                                                                                                                 <  2:735470/150‑1 (MQ=255)
    tttCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCaa                                                                                                                                               <  2:1259450/150‑1 (MQ=255)
    tttCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCaa                                                                                                                                               >  1:1217854/1‑150 (MQ=255)
      tCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAAGTTTTCAAcc                                                                                                                                             <  2:1279341/150‑1 (MQ=255)
       ccGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCCAccg                                                                                                                                            >  1:756048/1‑149 (MQ=255)
        cGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACctg                                                                                                                                           >  2:1608055/1‑150 (MQ=255)
          ccTGCCAGAGCTGCCAGCCCCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGCCGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACctgct                                                                                                                                         <  1:1301168/149‑1 (MQ=255)
           cagCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGATGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTAGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgcgc                                                                                                                                        <  1:817875/148‑3 (MQ=255)
           cTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACctgctg                                                                                                                                        <  2:611579/150‑1 (MQ=255)
           cTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgcgc                                                                                                                                        <  2:1368410/150‑3 (MQ=255)
            tGCCAGAGCTGCCAGCCTCCGCCCGCGCCGCGCAACAGCAGAAATCCATTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgcgca                                                                                                                                       <  2:1051245/150‑4 (MQ=255)
            tGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgcgca                                                                                                                                       <  1:1153067/150‑4 (MQ=255)
            tGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgcgca                                                                                                                                       <  1:1147229/150‑4 (MQ=255)
             gCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTTcc                                                                                                                                      >  2:1533972/1‑150 (MQ=255)
              ccAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCa                                                                                                                                     >  1:526563/1‑150 (MQ=255)
                  cgCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGcc                                                                                                                                 <  1:1202390/149‑1 (MQ=255)
                      gCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACt                                                                                                                             >  1:462009/1‑150 (MQ=255)
                        cAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTgg                                                                                                                           >  1:1594939/1‑150 (MQ=255)
                          gCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGtt                                                                                                                         <  2:646804/150‑1 (MQ=255)
                           ccTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGttt                                                                                                                        <  1:602750/150‑2 (MQ=255)
                            cTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGc                                                                                                                                  <  2:681586/139‑1 (MQ=255)
                              ccGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcc                                                                                                                     <  2:565936/150‑1 (MQ=255)
                              ccGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcc                                                                                                                     <  1:1629392/150‑1 (MQ=255)
                                  cAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACCCCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCTCACTGGTTatcctgat                                                                                                                 >  2:318491/1‑150 (MQ=255)
                                  cAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctgat                                                                                                                 <  1:550582/150‑1 (MQ=255)
                                   aGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctgatc                                                                                                                <  2:1130267/150‑1 (MQ=255)
                                   aGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctgatc                                                                                                                >  2:1225965/1‑150 (MQ=255)
                                    gcgcCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctgatcc                                                                                                               <  1:1677868/150‑1 (MQ=255)
                                       ccGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGGTATCCTGATCCTGt                                                                                                            >  1:501225/1‑150 (MQ=255)
                                         gAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTc                                                                                                          >  2:542229/1‑150 (MQ=255)
                                              aCAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATAcc                                                                                                     <  2:231266/150‑1 (MQ=255)
                                              aCAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATAcc                                                                                                     <  2:944006/150‑1 (MQ=255)
                                                agcagTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCAc                                                                                                   >  2:1439141/1‑150 (MQ=255)
                                                 gcagTATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGCTCCTGTTCATACCACTg                                                                                                 >  2:930391/1‑150 (MQ=255)
                                                  cagTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTg                                                                                                 <  2:1611824/150‑1 (MQ=255)
                                                      aaTCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCCCTGGGTATCCTGATCCTGTTCATACCACTGTTGc                                                                                             >  1:113443/1‑150 (MQ=255)
                                                       aTCCTTTCATCACGTGACGGCTGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCTGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGcc                                                                                            <  1:1265733/150‑1 (MQ=255)
                                                          cTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAgc                                                                                         <  2:432618/150‑1 (MQ=255)
                                                          cTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAgc                                                                                         <  1:1225965/150‑1 (MQ=255)
                                                              catcaCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAgcgcgc                                                                                     >  1:1119268/1‑150 (MQ=255)
                                                                tcaCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCtt                                                                                   >  2:608560/1‑150 (MQ=255)
                                                                 caCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGCTCCTGTTCATACCACTGTTGCCAGCGCGCTTc                                                                                  >  1:104521/1‑150 (MQ=255)
                                                                  aCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTcc                                                                                 <  1:552812/150‑1 (MQ=255)
                                                                   cGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCt                                                                                >  2:469466/1‑150 (MQ=255)
                                                                            gggTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGc                                                                       <  2:1319419/150‑1 (MQ=255)
                                                                              gTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTg                                                                     <  2:1110252/150‑1 (MQ=255)
                                                                                 aGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGAcg                                                                  >  1:949507/1‑150 (MQ=255)
                                                                                  ggCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCACCCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCCGCGCGCTTCCTGTTGTGGGCTGAcgc                                                                 >  1:1485520/1‑150 (MQ=255)
                                                                                  ggCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGAcgc                                                                 >  1:1231045/1‑150 (MQ=255)
                                                                                    cGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACgcgc                                                               <  1:1521499/150‑1 (MQ=255)
                                                                                     gggTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCGGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCGGACGCGCt                                                              <  1:1316244/150‑1 (MQ=255)
                                                                                        tATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGAGCTGACGCGCTcgc                                                           <  2:840206/150‑1 (MQ=255)
                                                                                           cGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGGTCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTcgccgc                                                        >  1:458614/1‑150 (MQ=255)
                                                                                               cAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAg                                                    >  2:1271587/1‑150 (MQ=255)
                                                                                                    ccccGCTGGCAACGTCGCCAGGCACACCGCCAAGCGGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCAGATCCTGTTCATACCACTGTTGCCCGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACgg                                               <  2:171113/150‑1 (MQ=255)
                                                                                                    cccGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACgg                                               >  2:946818/1‑150 (MQ=255)
                                                                                                     ccGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGATGTGGGCTGACGCGCTCGCCGCTTAGCACGGc                                              <  2:174022/150‑1 (MQ=255)
                                                                                                       ggCTGGCAACGGCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTGCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATCCCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCGTAGCACGGCAt                                            <  2:79812/150‑1 (MQ=255)
                                                                                                       gccTGGCGACGTTGCCAGGCACACCTCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCAt                                            >  2:292915/3‑150 (MQ=255)
                                                                                                        gCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTcc                                                                                 <  1:213944/112‑1 (MQ=255)
                                                                                                        gCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCACTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATa                                           >  1:330965/1‑150 (MQ=255)
                                                                                                         cTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCAGCTTAGCACGGCATAc                                          >  2:1218483/1‑150 (MQ=255)
                                                                                                          tGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATAcc                                         <  2:526563/150‑1 (MQ=255)
                                                                                                           ggCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGTGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCa                                        >  1:588848/1‑150 (MQ=255)
                                                                                                           ggCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCa                                        >  1:1143098/1‑150 (MQ=255)
                                                                                                             cAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGt                                      >  1:1497803/1‑150 (MQ=255)
                                                                                                             cAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGt                                      >  2:430400/1‑150 (MQ=255)
                                                                                                               aCGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTgg                                    <  1:927839/150‑1 (MQ=255)
                                                                                                               aCGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTgg                                    <  1:608560/150‑1 (MQ=255)
                                                                                                                 gTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGa                                  >  2:1708552/1‑150 (MQ=255)
                                                                                                                    gCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGCAAGc                               >  2:322411/1‑150 (MQ=255)
                                                                                                                     ccAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCt                              >  1:1402651/1‑150 (MQ=255)
                                                                                                                      cAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTg                             <  1:57515/150‑1 (MQ=255)
                                                                                                                        ggCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGaa                           >  1:408577/1‑150 (MQ=255)
                                                                                                                        ggCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGaa                           >  1:1124560/1‑150 (MQ=255)
                                                                                                                        ggCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGaa                           >  2:525854/1‑150 (MQ=255)
                                                                                                                         gCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAAc                          <  1:1533972/150‑1 (MQ=255)
                                                                                                                            cacCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCa                       <  1:430400/150‑1 (MQ=255)
                                                                                                                             acCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAg                      <  1:522621/150‑1 (MQ=255)
                                                                                                                             acCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAg                      >  1:1416817/1‑150 (MQ=255)
                                                                                                                             acCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCAAGTCCCCCTGGTTATCCTGATCCTGTTCATACCACTTTTGCCAGCGCGCTTCCGGTTCTGGGCTGACGCGCTATCCGCTTAGCACGGCATACCAGTTGGAAGCTGAACGCAg                      >  1:1034064/1‑150 (MQ=255)
                                                                                                                                 ccAAGCTGGTTGCGCATGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCAAGGCATACCAGTGGGAAGCTGAACGCAGCGcc                  <  2:1213963/150‑1 (MQ=255)
                                                                                                                                     gCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAc              <  1:1271587/150‑1 (MQ=255)
                                                                                                                                      cTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGCcg             >  1:1547252/1‑150 (MQ=255)
                                                                                                                                       tGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcgg            >  1:474715/1‑150 (MQ=255)
                                                                                                                                        ggTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcggc           >  1:140545/1‑150 (MQ=255)
                                                                                                                                        ggTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcggc           >  2:761732/1‑150 (MQ=255)
                                                                                                                                        ctggcagcAGGTTGAAAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcggc           >  1:1324107/7‑150 (MQ=255)
                                                                                                                                        ctggcagcAGGTTGAAAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcggc           >  1:258584/7‑150 (MQ=255)
                                                                                                                                           tGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCggg        >  1:975653/1‑150 (MQ=255)
                                                                                                                                           tGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCggg        <  2:761946/150‑1 (MQ=255)
                                                                                                                                            gcgcAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGa       >  2:1562940/1‑150 (MQ=255)
                                                                                                                                            cagcAGGTTGAAAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGa       >  2:686490/3‑150 (MQ=255)
                                                                                                                                             cgcAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGaa      >  1:4915/1‑150 (MQ=255)
                                                                                                                                                 ggTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCgg  <  2:1594939/150‑1 (MQ=255)
                                                                                                                                                 ggTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCgg  <  2:38410/150‑1 (MQ=255)
                                                                                                                                                 ggTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCgg  >  1:1145168/1‑150 (MQ=255)
                                                                                                                                                 ggTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCgg  <  1:1608055/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGG  >  NC_000913/4517407‑4517701

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