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

TGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGG  >  NC_000913/4516922‑4517377
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
tGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCgg                                                                                                                                                                                                                           <  1:995246/239‑1 (MQ=255)
ggCGGTAACGGGTTAGCGTGGCTAGCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCgg                                                                                                                                                                                                                           <  2:988597/238‑1 (MQ=255)
     tAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGt                                                                                                                                                                                                                      >  1:866083/1‑239 (MQ=255)
            ttAGCGTGCGTATCACCTCCCCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCt                                                                                                                                                                                                               <  2:415338/239‑1 (MQ=255)
             tAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCt                                                                                                                                                                                                               <  2:300855/238‑1 (MQ=255)
                     cTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTaaa                                                                                                                                                                                                      >  1:48665/1‑239 (MQ=255)
                        tccccgcaccCCGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCa                                                                                                                                                                                                   <  1:1151469/233‑1 (MQ=255)
                        tcacctcaccCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCa                                                                                                                                                                                                   <  2:966337/239‑1 (MQ=255)
                          acctcaccCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGt                                                                                                                                                                                                 >  2:455310/1‑239 (MQ=255)
                             tcaccCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAg                                                                                                                                                                                              <  1:295023/239‑1 (MQ=255)
                               accCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCATTGa                                                                                                                                                                                               >  1:296384/1‑236 (MQ=255)
                                            tCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCCCTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCCCTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               <  1:960179/239‑1 (MQ=255)
                                            tCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               <  2:48665/239‑1 (MQ=255)
                                            tCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               <  1:768075/239‑1 (MQ=255)
                                            tCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               <  1:222681/239‑1 (MQ=255)
                                             cGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               <  1:1128961/238‑1 (MQ=255)
                                                  aaGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCt                                                                                                                                                                                                               <  2:1007420/201‑1 (MQ=255)
                                                  aaGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCt                                                                                                                                                                                                               >  1:1007420/1‑201 (MQ=255)
                                                  aaGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGt                                                                                                                                                                                                             >  2:7392/1‑203 (MQ=255)
                                                  aaGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGt                                                                                                                                                                                                             <  1:7392/203‑1 (MQ=255)
                                                            tGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAgggg                                                                                                                                                               <  1:96227/239‑1 (MQ=255)
                                                              tCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               <  2:708335/221‑1 (MQ=255)
                                                              tCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt                                                                                                                                                                               >  1:708335/1‑221 (MQ=255)
                                                               ccTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGGGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGt                                                                                                                                                            >  1:708595/1‑239 (MQ=255)
                                                                             gTACTCTCGTCATCCAGCGGTTTCGCTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCttt                                                                                                                                              <  2:301837/239‑1 (MQ=255)
                                                                                   tcGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGa                                                                                                                                                                          >  2:251188/1‑205 (MQ=255)
                                                                                   tcGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGa                                                                                                                                                                          <  1:251188/205‑1 (MQ=255)
                                                                                   tcGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGgt                                                                                                                                        >  1:705144/1‑239 (MQ=255)
                                                                                                  ttCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTc                                                                                                                         >  1:746174/1‑239 (MQ=255)
                                                                                                    cACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCAc                                                                                                                       <  2:746174/239‑1 (MQ=255)
                                                                                                          gcCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAg                                                                                                                                                                            <  1:512584/180‑1 (MQ=255)
                                                                                                          gcCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAg                                                                                                                                                                            >  2:512584/1‑180 (MQ=255)
                                                                                                                 tCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCg                                                                                                          <  1:413208/239‑1 (MQ=255)
                                                                                                                   agagGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCgg                                                                                                                                                                                                                           <  1:662561/124‑1 (MQ=255)
                                                                                                                   agagGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCgg                                                                                                                                                                                                                           >  2:662561/1‑124 (MQ=255)
                                                                                                                   agagGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCa                                                                                                                <  1:804763/231‑1 (MQ=255)
                                                                                                                   agagGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCa                                                                                                                >  2:804763/1‑231 (MQ=255)
                                                                                                                      ggCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGa                                                                                                        <  1:455310/236‑1 (MQ=255)
                                                                                                                        cgcgggaCTGCGGGCTCTCACCAGCGTTCACGATGCGTTTTTGCGGGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGt                                                                                                    <  1:616527/232‑1 (MQ=255)
                                                                                                                          cTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTa                                                                                                                                                                                                        <  2:282173/136‑1 (MQ=255)
                                                                                                                          cTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTa                                                                                                                                                                                                        >  1:282173/1‑136 (MQ=255)
                                                                                                                                tGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATg                                                                                           <  2:866083/239‑1 (MQ=255)
                                                                                                                                  cAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAg                                                                                                                                                                  <  1:1005964/166‑1 (MQ=255)
                                                                                                                                  cAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAg                                                                                                                                                                  >  2:1005964/1‑166 (MQ=255)
                                                                                                                                  cAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATgc                                                                                          >  1:339278/1‑238 (MQ=255)
                                                                                                                                       tttCGCCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCa                                                                                    <  1:265114/239‑1 (MQ=255)
                                                                                                                                       tttCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGt                                                                                                                                                                                                             >  2:577068/1‑118 (MQ=255)
                                                                                                                                       tttCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGt                                                                                                                                                                                                             <  1:577068/118‑1 (MQ=255)
                                                                                                                                                     aCGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAgggg                                                                                                                                                               >  2:144077/1‑150 (MQ=255)
                                                                                                                                                     aCGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAgggg                                                                                                                                                               <  1:144077/150‑1 (MQ=255)
                                                                                                                                                          gCGTTTTTGCGCGGGGGCACTGGCGAGAAGCCTTTCCACAGCGTGCTGCTGCACGTCAAGGTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGGGCGGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTGCGCGGTAGTAATGGCGATTTCACCGTACCAGGGCCGGACGGTGCGCTGATGCGCATCAAACCGCCCATCCGCCgcgc                                                                 <  2:1144330/239‑1 (MQ=255)
                                                                                                                                                          gCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGa                                                                                                                                                                          >  1:559036/1‑134 (MQ=255)
                                                                                                                                                          gCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGa                                                                                                                                                                          <  2:559036/134‑1 (MQ=255)
                                                                                                                                                                 tGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGt                                                          >  1:656887/1‑239 (MQ=255)
                                                                                                                                                                                   aggaGCTCTTCCAGCGCGTGCTGCTGCTCGGAAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCCGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGCCGGTGCGCTGAGGCGCATCAAAACGCCCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAggg                                         <  2:593553/235‑1 (MQ=255)
                                                                                                                                                                                     aaGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCt                                                                                                                                                                                                               <  2:203287/70‑1 (MQ=255)
                                                                                                                                                                                     aaGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCt                                                                                                                                                                                                               >  1:203287/1‑70 (MQ=255)
                                                                                                                                                                                      aGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCCCCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAg                                       <  1:815928/237‑1 (MQ=255)
                                                                                                                                                                                                     tgctgctgcACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGGCACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGtt                          >  1:829447/1‑235 (MQ=255)
                                                                                                                                                                                                     tgctgctgcACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCCGGGAGCCATCATCCAGtt                          >  2:720376/1‑235 (MQ=255)
                                                                                                                                                                                                       ctgctgcACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGAc                    >  1:361351/1‑239 (MQ=255)
                                                                                                                                                                                                       ctgctgcACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGGATTCAGCGTGAGCAGGGTGCCATCTTCCAGGTGCTGAc                    >  2:326177/1‑239 (MQ=255)
                                                                                                                                                                                                        tgctgcACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCa                                                                                                                                                                                                   >  1:469922/1‑63 (MQ=255)
                                                                                                                                                                                                        tgctgcACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCa                                                                                                                                                                                                   <  2:469922/63‑1 (MQ=255)
                                                                                                                                                                                                         gctgcACGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACgg                  <  2:296384/239‑1 (MQ=255)
                                                                                                                                                                                                           tgcACGTCCAGTGCGTTGAAGGTACCCGGCTGGCGGACGGTAAGTCGTGTCCCTAAAGCAGTGAGGGGGCCCTGACGGGTGAGGACGCGGAAGGGGGGTGGCGGGGCATCTGGCGAGGTGGTAATGGCGATTTCACCGTGCGAGGGCCGGACGGTGCGCTGAGGCGCATCAAAACGCGCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCt                <  2:263001/239‑1 (MQ=255)
                                                                                                                                                                                                              aCGTCAAGCTGCGTGAAAGTATCCTGCGGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGGGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGAc             <  1:144182/239‑1 (MQ=255)
                                                                                                                                                                                                               cGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAg                                                                                                                                                                            <  1:15224/79‑1 (MQ=255)
                                                                                                                                                                                                               cGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAg                                                                                                                                                                            >  2:15224/1‑79 (MQ=255)
                                                                                                                                                                                                               cGTCAAGCTGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACg                   >  1:1124630/1‑232 (MQ=255)
                                                                                                                                                                                                                  cAAGCTGCGTGAAAGTATCCTGCTGGCGGACGTTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACtgt          <  1:894249/238‑1 (MQ=255)
                                                                                                                                                                                                                       tGCGTGAAAGTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGCATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACGGTGCCttt    >  2:314896/1‑239 (MQ=255)
                                                                                                                                                                                                                          gTGAAAGTAGCCGGCTGGCGGACGTTAAATTCTGTCCCGAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAACCGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTgg  <  2:48978/238‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
TGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGG  >  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: