Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,517,146 T→G N191T (AAT→ACT)  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,1460TG98.6% 272.4 / NA 73N191T (AAT→ACT) 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 (0/1);  new base G (30/42);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGG  >  NC_000913/4516922‑4517377
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
tGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCgg                                                                                                                                                                                                                           <  1:995246/239‑1 (MQ=255)
ggCGGTAACGGGTTAGCGTGGCTAGCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCgg                                                                                                                                                                                                                           <  2:988597/238‑1 (MQ=255)
     tAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGt                                                                                                                                                                                                                      >  1:866083/1‑239 (MQ=255)
            ttAGCGTGCGTATCACCTCCCCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCt                                                                                                                                                                                                               <  2:415338/239‑1 (MQ=255)
             tAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCt                                                                                                                                                                                                               <  2:300855/238‑1 (MQ=255)
                     cTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTaaa                                                                                                                                                                                                      >  1:48665/1‑239 (MQ=255)
                        tccccgcaccCCGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCa                                                                                                                                                                                                   <  1:1151469/233‑1 (MQ=255)
                        tcacctcaccCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCa                                                                                                                                                                                                   <  2:966337/239‑1 (MQ=255)
                          acctcaccCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGt                                                                                                                                                                                                 >  2:455310/1‑239 (MQ=255)
                             tcaccCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAg                                                                                                                                                                                              <  1:295023/239‑1 (MQ=255)
                               accCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCATTGa                                                                                                                                                                                               >  1:296384/1‑236 (MQ=255)
                                            tCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCCCTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCCCTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               <  1:960179/239‑1 (MQ=255)
                                            tCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               <  2:48665/239‑1 (MQ=255)
                                            tCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               <  1:768075/239‑1 (MQ=255)
                                            tCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               <  1:222681/239‑1 (MQ=255)
                                             cGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               <  1:1128961/238‑1 (MQ=255)
                                                  aaGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCt                                                                                                                                                                                                               <  2:1007420/201‑1 (MQ=255)
                                                  aaGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCt                                                                                                                                                                                                               >  1:1007420/1‑201 (MQ=255)
                                                  aaGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGt                                                                                                                                                                                                             >  2:7392/1‑203 (MQ=255)
                                                  aaGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGt                                                                                                                                                                                                             <  1:7392/203‑1 (MQ=255)
                                                            tGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAgggg                                                                                                                                                               <  1:96227/239‑1 (MQ=255)
                                                              tCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               <  2:708335/221‑1 (MQ=255)
                                                              tCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               >  1:708335/1‑221 (MQ=255)
                                                               ccTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGGGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGt                                                                                                                                                            >  1:708595/1‑239 (MQ=255)
                                                                             gTACTCTCGTCATCCAGCGGTTTCGCTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCttt                                                                                                                                              <  2:301837/239‑1 (MQ=255)
                                                                                   tcGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGa                                                                                                                                                                          >  2:251188/1‑205 (MQ=255)
                                                                                   tcGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGa                                                                                                                                                                          <  1:251188/205‑1 (MQ=255)
                                                                                   tcGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGgt                                                                                                                                        >  1:705144/1‑239 (MQ=255)
                                                                                                  ttCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTc                                                                                                                         >  1:746174/1‑239 (MQ=255)
                                                                                                    cACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCAc                                                                                                                       <  2:746174/239‑1 (MQ=255)
                                                                                                          gcCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAg                                                                                                                                                                            <  1:512584/180‑1 (MQ=255)
                                                                                                          gcCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAg                                                                                                                                                                            >  2:512584/1‑180 (MQ=255)
                                                                                                                 tCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCg                                                                                                          <  1:413208/239‑1 (MQ=255)
                                                                                                                   agagGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCgg                                                                                                                                                                                                                           <  1:662561/124‑1 (MQ=255)
                                                                                                                   agagGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCgg                                                                                                                                                                                                                           >  2:662561/1‑124 (MQ=255)
                                                                                                                   agagGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCa                                                                                                                <  1:804763/231‑1 (MQ=255)
                                                                                                                   agagGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCa                                                                                                                >  2:804763/1‑231 (MQ=255)
                                                                                                                      ggCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGa                                                                                                        <  1:455310/236‑1 (MQ=255)
                                                                                                                        cgcgggaCTGCGGGCTCTCACCAGCGTTCACGATGCGTTTTTGCGGGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGt                                                                                                    <  1:616527/232‑1 (MQ=255)
                                                                                                                          cTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTa                                                                                                                                                                                                        <  2:282173/136‑1 (MQ=255)
                                                                                                                          cTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTa                                                                                                                                                                                                        >  1:282173/1‑136 (MQ=255)
                                                                                                                                tGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATg                                                                                           <  2:866083/239‑1 (MQ=255)
                                                                                                                                  cAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAg                                                                                                                                                                  <  1:1005964/166‑1 (MQ=255)
                                                                                                                                  cAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAg                                                                                                                                                                  >  2:1005964/1‑166 (MQ=255)
                                                                                                                                  cAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATgc                                                                                          >  1:339278/1‑238 (MQ=255)
                                                                                                                                       tttCGCCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCa                                                                                    <  1:265114/239‑1 (MQ=255)
                                                                                                                                       tttCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGt                                                                                                                                                                                                             >  2:577068/1‑118 (MQ=255)
                                                                                                                                       tttCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGt                                                                                                                                                                                                             <  1:577068/118‑1 (MQ=255)
                                                                                                                                                     aCGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAgggg                                                                                                                                                               >  2:144077/1‑150 (MQ=255)
                                                                                                                                                     aCGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAgggg                                                                                                                                                               <  1:144077/150‑1 (MQ=255)
                                                                                                                                                          gCGTTTTTGCGCGGGGGCACTGGCGAGAAGCCTTTCCACAGCGTGCTGCTGCACGTCAAGGTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGGGCGGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTGCGCGGTAGTAATGGCGATTTCACCGTACCAGGGCCGGACGGTGCGCTGATGCGCATCAAACCGCCCATCCGCCgcgc                                                                 <  2:1144330/239‑1 (MQ=255)
                                                                                                                                                          gCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGa                                                                                                                                                                          >  1:559036/1‑134 (MQ=255)
                                                                                                                                                          gCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGa                                                                                                                                                                          <  2:559036/134‑1 (MQ=255)
                                                                                                                                                                 tGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGt                                                          >  1:656887/1‑239 (MQ=255)
                                                                                                                                                                                   aggaGCTCTTCCAGCGCGTGCTGCTGCTCGGAAAGCTGCG‑TGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCCGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGCCGGTGCGCTGAGGCGCATCAAAACGCCCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAggg                                         <  2:593553/235‑1 (MQ=255)
                                                                                                                                                                                     aaGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCt                                                                                                                                                                                                               <  2:203287/70‑1 (MQ=255)
                                                                                                                                                                                     aaGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCt                                                                                                                                                                                                               >  1:203287/1‑70 (MQ=255)
                                                                                                                                                                                      aGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCCCCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAg                                       <  1:815928/237‑1 (MQ=255)
                                                                                                                                                                                                     tgctgctgcACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGGCACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGtt                          >  1:829447/1‑235 (MQ=255)
                                                                                                                                                                                                     tgctgctgcACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCCGGGAGCCATCATCCAGtt                          >  2:720376/1‑235 (MQ=255)
                                                                                                                                                                                                       ctgctgcACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGAc                    >  1:361351/1‑239 (MQ=255)
                                                                                                                                                                                                       ctgctgcACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGGATTCAGCGTGAGCAGGGTGCCATCTTCCAGGTGCTGAc                    >  2:326177/1‑239 (MQ=255)
                                                                                                                                                                                                        tgctgcACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCa                                                                                                                                                                                                   >  1:469922/1‑63 (MQ=255)
                                                                                                                                                                                                        tgctgcACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCa                                                                                                                                                                                                   <  2:469922/63‑1 (MQ=255)
                                                                                                                                                                                                         gctgcACGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACgg                  <  2:296384/239‑1 (MQ=255)
                                                                                                                                                                                                           tgcACGTCCAG‑TGCGTTGAAGGTACCCGGCTGGCGGACGGTAAGTCGTGTCCCTAAAGCAGTGAGGGGGCCCTGACGGGTGAGGACGCGGAAGGGGGGTGGCGGGGCATCTGGCGAGGTGGTAATGGCGATTTCACCGTGCGAGGGCCGGACGGTGCGCTGAGGCGCATCAAAACGCGCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCt                <  2:263001/239‑1 (MQ=255)
                                                                                                                                                                                                              aCGTCAAGCTGCG‑TGAAAGTATCCTGCGGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGGGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGAc             <  1:144182/239‑1 (MQ=255)
                                                                                                                                                                                                               cGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAg                                                                                                                                                                            <  1:15224/79‑1 (MQ=255)
                                                                                                                                                                                                               cGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAg                                                                                                                                                                            >  2:15224/1‑79 (MQ=255)
                                                                                                                                                                                                               cGTCAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACg                   >  1:1124630/1‑232 (MQ=255)
                                                                                                                                                                                                                  cAAGCTGCG‑TGAAAGTATCCTGCTGGCGGACGTTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACtgt          <  1:894249/238‑1 (MQ=255)
                                                                                                                                                                                                                       tGCG‑TGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGCATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACGGTGCCttt    >  2:314896/1‑239 (MQ=255)
                                                                                                                                                                                                                          g‑TGAAAGTAGCCGGCTGGCGGACGTTAAATTCTGTCCCGAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAACCGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTgg  <  2:48978/238‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
TGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCG‑TGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGG  >  NC_000913/4516922‑4517377

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑