Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,380,292 2 bp→A coding (225‑226/735 nt) frdB ← fumarate reductase iron‑sulfur protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,380,2920G.94.7% 158.7 / 3.3 38coding (226/735 nt)frdBfumarate reductase iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (22/14);  minor base A (1/1);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
*NC_0009134,380,2930TA97.2% 127.0 / NA 36K75N (AAA→AAT) frdBfumarate reductase iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (21/14);  minor base C (1/0);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTC  >  NC_000913/4380045‑4380531
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
gcggcgTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAAttt                                                                                                                                                                                                                                              <  2:131440/251‑1 (MQ=255)
   gcgTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAAtttggc                                                                                                                                                                                                                                           >  1:150672/1‑248 (MQ=255)
              aCCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCGTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAAc                                                                                                                                                                                                                               <  2:268240/251‑1 (MQ=255)
                                  agaACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACcac                                                                                                                                                                                                           <  1:252579/251‑1 (MQ=255)
                                       tGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTGCGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATc                                                                                                                                                                                                      >  2:142709/1‑251 (MQ=255)
                                            attcTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCAt                                                                                                                                                                                                 <  2:216479/248‑1 (MQ=255)
                                                           cgcCGGGGTCTGGATGTTAGTACCCTGATCCGCGGAGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTCTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCg                                                                                                                                                                                  >  1:255159/1‑251 (MQ=255)
                                                            gcCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCgt                                                                                                                                                                                 <  2:248330/251‑2 (MQ=255)
                                                             ccGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCACCATACGGCAGGACCAGCggg                                                                                                                                                                                >  1:7104/1‑250 (MQ=255)
                                                              cGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGt                                                                                                                                                                                >  2:203685/1‑250 (MQ=255)
                                                                     tGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAgg                                                                                                                                                                        >  1:78802/1‑251 (MQ=255)
                                                                       gATGTTATTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGGTGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTc                                                                                                                                                                      <  2:55499/251‑1 (MQ=255)
                                                                            tAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt                                                                                                                                                                  <  2:78802/250‑1 (MQ=255)
                                                                                     gATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGt                                                                                                                                                        >  1:107844/1‑251 (MQ=255)
                                                                                                     gtgtgcGATGATGTACGGTTTGATCCCTTCCAGACGTGCGATGAAGTGGGTCATATCGATCACCAGATCGCGTGCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCCTCGTTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGTAGGAACCACAAAGCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTTATGTAGGCCAg                                                                                                                                        <  2:145534/246‑1 (MQ=255)
                                                                                                                  tACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAAt                                                                                                                           >  1:184664/1‑251 (MQ=255)
                                                                                                                    cGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGa                                                                                                                         <  2:5103/251‑1 (MQ=255)
                                                                                                                                   cAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATa                                                                                                          <  2:210174/251‑1 (MQ=255)
                                                                                                                                                    gggTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATg                                                                                         <  1:82228/251‑1 (MQ=255)
                                                                                                                                                     ggTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGGAGTTGCGTCATAAGGCCCTTCATAGAATgc                                                                                        >  2:238909/1‑251 (MQ=255)
                                                                                                                                                       tCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATgcgc                                                                                      >  2:19720/1‑251 (MQ=255)
                                                                                                                                                         atatCGACCACCAGATCGCGTTCAATCGGGAAGTTCGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCCTCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCt                                                                                     >  2:18028/1‑250 (MQ=255)
                                                                                                                                                          tatCGACCCCCCGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTAt                                                                                   >  1:12118/1‑251 (MQ=255)
                                                                                                                                                                    caGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATgcggtgcggt                                                                         <  1:203685/251‑1 (MQ=255)
                                                                                                                                                                          gcgTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGAc                                                                   >  2:258529/1‑251 (MQ=255)
                                                                                                                                                                           cgTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACt                                                                  >  1:12174/1‑251 (MQ=255)
                                                                                                                                                                             ttCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTc                                                                >  1:83734/1‑251 (MQ=255)
                                                                                                                                                                                     ggaagggaGCTAACGCTTCAACCTTCATACCGTCTGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAAGCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAAGTCGCAATGCGGTGCGGTATCGACTTCCGGGtgat                                                        <  2:165062/244‑4 (MQ=255)
                                                                                                                                                                                     ggAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTtat                                                        >  1:162505/1‑251 (MQ=255)
                                                                                                                                                                                            aGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGTTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGcac                                                 >  1:231869/1‑251 (MQ=255)
                                                                                                                                                                                                 aCGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCt                                            <  2:105347/251‑1 (MQ=255)
                                                                                                                                                                                                  cGCTTCAACCGTCATACCGTCGGTGTAAGCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACTGTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATCCGCTATGCGGTGCGGTATCGCCTTCCGGGTTATAGCGCACCACCTc                                           <  2:205204/251‑1 (MQ=255)
                                                                                                                                                                                                         aCCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCTA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATTTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCACTTTTc                                    >  2:245523/1‑251 (MQ=255)
                                                                                                                                                                                                           cTTCATACCGTCGGTGTAATCACGCAGGGAGGTTTTACCTGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCCGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTc                                    >  1:177433/1‑249 (MQ=255)
                                                                                                                                                                                                                   cGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTc                           >  1:255482/1‑250 (MQ=255)
                                                                                                                                                                                                                          gtAATCACGCAGGAAGGTTTTACATGCCA‑CTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGGGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGc                   >  2:60053/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ggAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGccttc        >  1:53260/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ttttACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGGTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGccttctccttc  >  2:232490/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
GCGGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTC  >  NC_000913/4380045‑4380531

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