Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,517,3550CT56.0% 138.2 / inf 414Q121Q (CAG→CAA) 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 C (85/97);  new base T (124/108);  total (209/205)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 421 total aligned reads displayed.

CAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCA  >  NC_000913/4517205‑4517503
                                                                                                                                                      |                                                                                                                                                    
cAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGc                                                                                                                                                      >  8:1412676/1‑151 (MQ=255)
 aGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGGAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCt                                                                                                                                                     >  3:448192/1‑151 (MQ=255)
  ggACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGttgttg                                                                                                                                                    >  8:1268822/1‑151 (MQ=255)
    acacGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGAc                                                                                                                                                  >  1:1317790/1‑151 (MQ=255)
     cacGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACg                                                                                                                                                 <  5:1542589/151‑1 (MQ=255)
        ggAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCt                                                                                                                                              <  7:707965/151‑1 (MQ=255)
        ggAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCt                                                                                                                                              <  7:658682/151‑1 (MQ=255)
        ggAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGCACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGCAGCCATCTTCCAGTTGTTGACGGCt                                                                                                                                              >  6:1102974/1‑151 (MQ=255)
         gAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTg                                                                                                                                             <  2:1415713/151‑1 (MQ=255)
          aaGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTg                                                                                                                                             >  6:405756/1‑150 (MQ=255)
            ggggCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACt                                                                                                                                          <  1:791229/151‑1 (MQ=255)
                cGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGc                                                                                                                                      >  5:466597/1‑151 (MQ=255)
                   tGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCtt                                                                                                                                   >  3:194149/1‑151 (MQ=255)
                    gCAGGGCATCTTTCGCGGTTGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCttt                                                                                                                                  <  4:1667825/151‑1 (MQ=255)
                     cAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTg                                                                                                                                 <  8:1585718/151‑1 (MQ=255)
                        ggCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcg                                                                                                                              <  3:116427/151‑1 (MQ=255)
                        ggCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGgcg                                                                                                                              >  3:1526909/1‑151 (MQ=255)
                         gCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcgg                                                                                                                             <  4:1025582/151‑1 (MQ=255)
                          cATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggt                                                                                                                            >  3:1614898/1‑151 (MQ=255)
                            tCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGACCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGgcggtgc                                                                                                                          >  3:889230/1‑151 (MQ=255)
                              tttCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtgcgg                                                                                                                        >  8:1504832/1‑151 (MQ=255)
                               ttCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtgc                                                                                                                          <  4:402182/148‑1 (MQ=255)
                                 cgcgGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTa                                                                                                                      <  3:1378968/150‑1 (MQ=255)
                                 cgcgGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTa                                                                                                                      <  3:1028853/150‑1 (MQ=255)
                                       agtaATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGGATTCAGCGTGAGCAGGGAGCCATCTTCGAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGc                                                                                                                <  1:985141/150‑1 (MQ=255)
                                         taATTGCGATTTCACCGTACCAGAGCCGTACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGAAGGATGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCg                                                                                                             >  4:964148/1‑151 (MQ=255)
                                          aaTGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCg                                                                                                             <  1:1094743/150‑1 (MQ=255)
                                          aaTGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCg                                                                                                             <  6:227456/150‑1 (MQ=255)
                                          aaTGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGc                                                                                                            >  4:522299/1‑151 (MQ=255)
                                                gATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCt                                                                                                      <  5:410427/151‑1 (MQ=255)
                                                   ttCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCg                                                                                                   >  7:328745/1‑151 (MQ=255)
                                                   ttCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTc                                                                                                    <  5:1344395/150‑1 (MQ=255)
                                                   ttCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTc                                                                                                    <  2:1652950/150‑1 (MQ=255)
                                                    tCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCg                                                                                                   <  5:1119890/150‑1 (MQ=255)
                                                    tCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCg                                                                                                   <  5:1659803/150‑1 (MQ=255)
                                                      aCCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCg                                                                                                >  5:56641/1‑151 (MQ=255)
                                                      aCCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCg                                                                                                >  8:221407/1‑151 (MQ=255)
                                                      aCCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTAAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCg                                                                                                <  3:644504/151‑1 (MQ=255)
                                                       ccGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCgg                                                                                               >  1:1024199/1‑151 (MQ=255)
                                                        cGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGt                                                                                              >  8:680415/1‑151 (MQ=255)
                                                         gTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGtt                                                                                             >  3:1202527/1‑151 (MQ=255)
                                                           aCCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGttt                                                                                            <  5:763210/150‑1 (MQ=255)
                                                                  ccGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGc                                                                                    >  2:1310661/1‑151 (MQ=255)
                                                                   cGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGcc                                                                                   >  1:533307/1‑151 (MQ=255)
                                                                    ggACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGTGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCa                                                                                  <  8:1362786/151‑1 (MQ=255)
                                                                      aCGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCaga                                                                                <  7:1184210/151‑1 (MQ=255)
                                                                            cgcTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTCTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTgac                                                                          >  6:776752/1‑149 (MQ=255)
                                                                               tGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGc                                                                       >  1:416481/1‑151 (MQ=255)
                                                                               tGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGc                                                                       >  1:173460/1‑151 (MQ=255)
                                                                                gATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGcc                                                                      >  6:97419/1‑151 (MQ=255)
                                                                                gATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGcc                                                                      <  6:481964/151‑1 (MQ=255)
                                                                                   gcgcATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTcc                                                                   >  7:835439/1‑151 (MQ=255)
                                                                                     gcATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGc                                                                 >  5:1439392/1‑151 (MQ=255)
                                                                                       aTCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGcc                                                                >  6:1167190/1‑150 (MQ=255)
                                                                                       aTCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGcc                                                                >  4:1584800/1‑150 (MQ=255)
                                                                                         cataaCGCACTTCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGGCTTCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCGCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTGCGCCAgc                                                             <  2:1102723/148‑1 (MQ=255)
                                                                                         cAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAgc                                                             <  5:195946/151‑1 (MQ=255)
                                                                                         cAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAgc                                                             >  7:1680388/1‑151 (MQ=255)
                                                                                          aaaaCGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAgcg                                                            >  5:334110/1‑151 (MQ=255)
                                                                                          aaaaCGCACATCCGCCGCGCTTTGGGTAGTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAgcg                                                            <  2:551241/151‑1 (MQ=255)
                                                                                           aaaCGCACATCCGCCGCGGTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGTTAATCTGCCCGCAGACCTTCGCGGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAgcgc                                                           <  4:1491254/151‑1 (MQ=255)
                                                                                             atgCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCg                                                         >  7:132209/3‑151 (MQ=255)
                                                                                              cGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGa                                                        <  5:1285017/151‑1 (MQ=255)
                                                                                                cacaTCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGc                                                      <  8:1194129/151‑1 (MQ=255)
                                                                                                   aTCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAc                                                   <  4:194193/151‑1 (MQ=255)
                                                                                                       gccgcGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCCGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAcagc                                                >  7:389509/1‑150 (MQ=255)
                                                                                                        ccgcGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAAcagcag                                              <  8:859910/151‑1 (MQ=255)
                                                                                                         cgcgCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGt                                             >  4:1095074/1‑151 (MQ=255)
                                                                                                         cgcgCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGt                                             <  7:271254/151‑1 (MQ=255)
                                                                                                           cgcTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTa                                            >  1:145672/1‑150 (MQ=255)
                                                                                                            gcTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACTGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAAt                                          <  1:1228897/151‑1 (MQ=255)
                                                                                                               ttGGGTGTTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCt                                       >  4:810319/1‑151 (MQ=255)
                                                                                                               ttGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATc                                         <  8:1036770/149‑1 (MQ=255)
                                                                                                                tGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGACCGCAGACCTTCTCCGGTTTCCGACTGCCAGCGCTGCCATCCTCCGCAAGCGCCGAGCAACAGCAGTAATCCtt                                      >  7:1633657/1‑151 (MQ=255)
                                                                                                                 gggTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGATTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCttt                                     <  1:321782/151‑1 (MQ=255)
                                                                                                                  ggTATGCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCAGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTtc                                    <  4:714147/151‑1 (MQ=255)
                                                                                                                   gTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCCGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTtca                                   <  4:1275890/151‑1 (MQ=255)
                                                                                                                    tATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTtcat                                  >  2:613296/1‑151 (MQ=255)
                                                                                                                     aTTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTtcat                                  <  8:1652021/150‑1 (MQ=255)
                                                                                                                         aGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGTGCCGAGCAACAGCAGTAATCCTTTCATCACGt                             >  6:572408/1‑151 (MQ=255)
                                                                                                                         aGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGt                             <  8:167108/151‑1 (MQ=255)
                                                                                                                          gCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTg                            <  5:930055/151‑1 (MQ=255)
                                                                                                                                  aGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCggg                    >  6:1041376/1‑151 (MQ=255)
                                                                                                                                  aGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGCCCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCggg                    >  4:303915/1‑151 (MQ=255)
                                                                                                                                  aGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCggg                    >  8:327107/1‑151 (MQ=255)
                                                                                                                                     gAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGa                 <  6:1389098/151‑1 (MQ=255)
                                                                                                                                       gCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAgg               >  2:449737/1‑151 (MQ=255)
                                                                                                                                       gCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAgg               >  3:741801/1‑151 (MQ=255)
                                                                                                                                        ccATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGc              <  2:216895/151‑1 (MQ=255)
                                                                                                                                         cATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCg             <  1:3921/151‑1 (MQ=255)
                                                                                                                                         cATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCg             >  3:30234/1‑151 (MQ=255)
                                                                                                                                          agcTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCgg            <  5:974994/149‑1 (MQ=255)
                                                                                                                                           tCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCggg           >  5:851579/1‑151 (MQ=255)
                                                                                                                                             ttCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGt          >  6:1536356/1‑150 (MQ=255)
                                                                                                                                             ttCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGt          >  8:474148/1‑150 (MQ=255)
                                                                                                                                             ttCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGt          >  8:635012/1‑150 (MQ=255)
                                                                                                                                               ccaggcgctgacggctgacgGTCCCGTTGGCGGCGCGGGAAGCTCCCCGCAGACCTTCGCCGGGTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATc       <  7:1472103/145‑1 (MQ=255)
                                                                                                                                                cAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCg      <  7:506215/151‑1 (MQ=255)
                                                                                                                                                  gTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTg    <  6:316620/151‑1 (MQ=255)
                                                                                                                                                    tGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGc   <  2:75534/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCA  >  NC_000913/4517205‑4517503

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: ‑