Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,517,355 C→T Q121Q (CAG→CAA 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,3550CT97.4% 1436.7 / 15.2 428Q121Q (CAG→CAAfecRanti‑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 C (8/3);  new base T (206/209);  total (215/213)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 434 total aligned reads displayed.

GGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGC  >  NC_000913/4517207‑4517502
                                                                                                                                                    |                                                                                                                                                   
ggACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGttgttg                                                                                                                                                   <  6:495423/151‑1 (MQ=255)
  acacGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGAc                                                                                                                                                 >  1:141053/1‑151 (MQ=255)
    acGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACgg                                                                                                                                               >  6:1156379/1‑151 (MQ=255)
     cGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGc                                                                                                                                              >  7:767039/1‑151 (MQ=255)
      ggAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGcg                                                                                                                                             >  4:251605/1‑150 (MQ=255)
      ggAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCt                                                                                                                                             >  3:1038636/1‑151 (MQ=255)
              cGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTGCAGTTGTTGACGGCTGACTGTGc                                                                                                                                     >  2:1142089/1‑151 (MQ=255)
               gTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAACCGCACATCCGCCGCGCTTTGGGTATTCCGCGTGCGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGcc                                                                                                                                    >  6:666095/1‑151 (MQ=255)
                ttGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCt                                                                                                                                   >  2:986577/1‑151 (MQ=255)
                  gCATGGCATCGTTCGCGGTAGTAATGGCGAATTCAACGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCttt                                                                                                                                 <  5:90008/151‑1 (MQ=255)
                  gCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCttt                                                                                                                                 >  5:1275401/1‑151 (MQ=255)
                  gCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCttt                                                                                                                                 >  1:920279/1‑151 (MQ=255)
                   cAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTg                                                                                                                                >  2:316448/1‑151 (MQ=255)
                    aGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTgg                                                                                                                               >  1:449994/1‑151 (MQ=255)
                        cATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggt                                                                                                                           >  7:1283975/1‑151 (MQ=255)
                        cATCTTTCGCGGTAGTAATCGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGGGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggt                                                                                                                           <  1:1446742/151‑1 (MQ=255)
                             ttCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtgcggt                                                                                                                      >  6:1064003/1‑151 (MQ=255)
                              tCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtgcggt                                                                                                                      <  3:820015/150‑1 (MQ=255)
                               cgcgGTAGTAATGGCGGTTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTaa                                                                                                                    <  5:549894/151‑1 (MQ=255)
                               cgcgGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTa                                                                                                                     <  1:245375/150‑1 (MQ=255)
                                gcgGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAAt                                                                                                                   >  4:573903/1‑151 (MQ=255)
                                 cgGTAGTAATGGCGATTTCACCGTACCAGAGCCGGCCGGTGCGCTGATGCGCATCAAAACGCCCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCAACTTCCAGTTCTTGACGGCTGACTGTGCCATTGGCGGTGCGGTAATc                                                                                                                  >  4:1138954/1‑151 (MQ=255)
                                  ggTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCt                                                                                                                 <  4:915847/151‑1 (MQ=255)
                                     agtaATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGc                                                                                                               <  5:122159/150‑1 (MQ=255)
                                      gtaATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGccc                                                                                                             >  2:685433/1‑151 (MQ=255)
                                        aaTGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGc                                                                                                           >  1:556802/1‑151 (MQ=255)
                                           ggCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAg                                                                                                         <  5:449580/150‑1 (MQ=255)
                                            gCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGAc                                                                                                       >  8:1280552/1‑151 (MQ=255)
                                              gATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCt                                                                                                     >  7:160593/1‑151 (MQ=255)
                                              gATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCt                                                                                                     <  7:73600/151‑1 (MQ=255)
                                                   cACCGGACCAGAGCCGGACGGTGCGCTGATGCGCATCAAACCGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGcc                                                                                                <  7:905873/151‑1 (MQ=255)
                                                         aCCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGttt                                                                                           <  3:1624459/150‑1 (MQ=255)
                                                            tgagCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCg                                                                                        >  3:1505078/2‑150 (MQ=255)
                                                            agagCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCgc                                                                                       >  3:1303261/1‑150 (MQ=255)
                                                                ccGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTGAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGc                                                                                   <  3:1669399/151‑1 (MQ=255)
                                                                       gTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCt                                                                            >  1:1302657/1‑151 (MQ=255)
                                                                        tGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTg                                                                           >  3:56525/1‑151 (MQ=255)
                                                                         gcgcTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGc                                                                          >  1:1196808/1‑151 (MQ=255)
                                                                          cgcTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGcc                                                                         >  2:1222732/1‑151 (MQ=255)
                                                                          cgcTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGcc                                                                         >  4:719262/1‑151 (MQ=255)
                                                                           gcTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGcc                                                                         >  6:986568/1‑150 (MQ=255)
                                                                           gcTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCa                                                                        <  8:1486099/151‑1 (MQ=255)
                                                                              gATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGcc                                                                     <  1:1003330/150‑1 (MQ=255)
                                                                               aTGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCt                                                                    <  4:1329890/151‑1 (MQ=255)
                                                                                 gcgcATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTcc                                                                  <  5:1098294/151‑1 (MQ=255)
                                                                                   gcATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGc                                                                >  6:1524761/1‑151 (MQ=255)
                                                                                    cATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGcc                                                               >  8:1113479/1‑151 (MQ=255)
                                                                                     aTCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGcc                                                               >  2:1477889/1‑150 (MQ=255)
                                                                                          aaCGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGcc                                                         >  4:1548139/1‑151 (MQ=255)
                                                                                          aaCGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGcc                                                         >  4:486060/1‑151 (MQ=255)
                                                                                           aCGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCg                                                        >  8:1339791/1‑151 (MQ=255)
                                                                                            cGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGa                                                       <  7:1388066/151‑1 (MQ=255)
                                                                                              cacaTCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGc                                                     >  5:237420/1‑151 (MQ=255)
                                                                                                 aTCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAc                                                  <  5:1263287/151‑1 (MQ=255)
                                                                                                 aTCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAAAGCTGCCAGCCTCCGCCAGCGCCGAGCAAc                                                  >  7:1187779/1‑151 (MQ=255)
                                                                                                   ccgccgCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAca                                                 <  3:1171981/150‑1 (MQ=255)
                                                                                                     gccgcGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAcagca                                              <  5:755500/151‑1 (MQ=255)
                                                                                                      ccgcGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAcagcag                                             <  5:358657/151‑1 (MQ=255)
                                                                                                      ccgcGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCACCAGCGCCGAGCAAcagcag                                             <  1:1094854/151‑1 (MQ=255)
                                                                                                        gcgcTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTACCAGTTCTTGACGCCTGACTCTGCCTTTGGCGGTGCGGTAATCTGCCAGCAGACCTTCGCCGCTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGt                                            >  6:1628421/1‑150 (MQ=255)
                                                                                                            tttGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATcc                                       >  8:354187/1‑151 (MQ=255)
                                                                                                            tttGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATcc                                       >  6:745460/1‑151 (MQ=255)
                                                                                                             ttGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATcc                                       <  1:1686605/150‑1 (MQ=255)
                                                                                                              tGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCtt                                     >  2:626169/1‑151 (MQ=255)
                                                                                                              tGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCt                                      <  7:237840/150‑1 (MQ=255)
                                                                                                               gggTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCttt                                    <  6:567396/151‑1 (MQ=255)
                                                                                                               gggTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCtt                                     <  5:1503606/150‑1 (MQ=255)
                                                                                                                ggTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTtc                                   <  6:387033/151‑1 (MQ=255)
                                                                                                                   aTTCAGCGTGAGCAGGGAGCGATCTTCCTGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCGGCCAGCGCCGAGCAACAGCTGTAATCCTTtcat                                 <  8:1703690/150‑1 (MQ=255)
                                                                                                                   aTTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTtcatc                                <  3:406563/151‑1 (MQ=255)
                                                                                                                     tCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCAc                              >  3:114241/1‑151 (MQ=255)
                                                                                                                     tCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCACAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCAc                              >  3:121819/1‑151 (MQ=255)
                                                                                                                     tCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGCGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGGGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTtcatca                               <  4:964195/150‑1 (MQ=255)
                                                                                                                      cAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACg                             <  7:990707/151‑1 (MQ=255)
                                                                                                                       aGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGt                            >  7:65618/1‑151 (MQ=255)
                                                                                                                       aGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGt                            >  3:693619/1‑151 (MQ=255)
                                                                                                                        gCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTg                           <  2:850769/151‑1 (MQ=255)
                                                                                                                        gCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTg                           <  2:576517/151‑1 (MQ=255)
                                                                                                                         cGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTg                           >  3:1454936/1‑150 (MQ=255)
                                                                                                                           tGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGAcg                        >  2:1546139/1‑151 (MQ=255)
                                                                                                                           tGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGAcg                        >  8:499977/1‑151 (MQ=255)
                                                                                                                           ggAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGAcg                        <  1:1594799/150‑1 (MQ=255)
                                                                                                                            gAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGAcgg                       >  8:719994/1‑151 (MQ=255)
                                                                                                                             aGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGAcggc                      <  3:308431/151‑1 (MQ=255)
                                                                                                                                 gggAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGt                  >  6:277127/1‑151 (MQ=255)
                                                                                                                                  ggAGCCATCTTCCAGTTGTGGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTg                 <  5:1626508/151‑1 (MQ=255)
                                                                                                                                    aGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTg                 >  1:76310/1‑149 (MQ=255)
                                                                                                                                       cATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCg            <  6:1308853/151‑1 (MQ=255)
                                                                                                                                        aTCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCgg           >  8:656388/1‑151 (MQ=255)
                                                                                                                                        aTCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCgg           >  3:914826/1‑151 (MQ=255)
                                                                                                                                          cTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGCCCTTCGCCGGAGTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCGTTTCATCACGTGACGGCGGGTGAGGCGGGt         >  4:1670230/1‑151 (MQ=255)
                                                                                                                                           ttCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGt         >  8:1230014/1‑150 (MQ=255)
                                                                                                                                            tCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCgg           <  8:836474/147‑1 (MQ=255)
                                                                                                                                            tCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCgg           <  8:1619911/147‑1 (MQ=255)
                                                                                                                                             gcAGTTGTTGTCGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGGAGACCTTCGGCTGTTTCCGACTGGCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAGAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTTTc      <  1:342446/150‑1 (MQ=255)
                                                                                                                                              cAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCg     <  4:820008/151‑1 (MQ=255)
                                                                                                                                              cAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCg     >  7:510153/1‑151 (MQ=255)
                                                                                                                                              cAGTGGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCg     <  4:323945/151‑1 (MQ=255)
                                                                                                                                                gttgttGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTg   >  2:459950/1‑151 (MQ=255)
                                                                                                                                                gttgttGACGGCTGACTGTGCCTTTGGCGGGGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTg   <  6:1246296/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGC  >  NC_000913/4517207‑4517502

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