Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 323,062 A→G 7.8% E242G (GAG→GGG)  ykgF → ferridoxin‑like LutB family protein; putative electron transport chain YkgEFG component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913323,0620AG 7.8% 184.4 / 4.0 63E242G (GAG→GGG) ykgFferridoxin‑like LutB family protein; putative electron transport chain YkgEFG component
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base G (3/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACCGCGCGAAGCTGGGCACGTTGATGGTCCTGAAGAGTTTCATCT  >  NC_000913/322916‑323209
                                                                                                                                                  |                                                                                                                                                   
cgAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGc                                                                                                                                                   <  8:367275/149‑1 (MQ=255)
   aGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTa                                                                                                                                                <  7:360320/149‑1 (MQ=255)
     aTTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTAtt                                                                                                                                              >  1:226759/1‑149 (MQ=255)
       ttCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGGGATGGGAATGGAGCGTATTGc                                                                                                                                            >  1:343687/1‑149 (MQ=255)
        tCCTCAGTGCTTAAATAGGTATTACCTGCTGTAATTTCGCGTTGGCAGAGACCTGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGc                                                                                                                                            <  6:163440/148‑1 (MQ=255)
          cTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGGGCGTATTGcccc                                                                                                                                         >  6:361859/1‑149 (MQ=255)
          cTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGc                                                                                                                                                   >  7:412317/1‑139 (MQ=255)
          cTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGAGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGc                                                                                                                                                   <  8:412317/139‑1 (MQ=255)
             aGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCAc                                                                                                                                      >  1:107550/1‑149 (MQ=255)
                 cTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGttt                                                                                                                                  >  5:309433/1‑149 (MQ=255)
                      aTAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGcc                                                                                                                               >  4:370137/1‑147 (MQ=255)
                      aTAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGa                                                                                                                             <  3:186514/149‑1 (MQ=255)
                      aTAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGa                                                                                                                             <  7:36918/149‑1 (MQ=255)
                       tAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAg                                                                                                                            <  6:309433/149‑1 (MQ=255)
                         ggTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGt                                                                                                                          <  8:38871/149‑1 (MQ=255)
                         ggTATTACCGGCTGTAAATTCGCGGTGGCAGAGACCTGTTCGGTATGCCTGGGGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGt                                                                                                                          <  5:238854/149‑1 (MQ=255)
                          gTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAgg                                                                                                                           >  1:286163/1‑147 (MQ=255)
                           tATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAgg                                                                                                                           >  3:380061/1‑146 (MQ=255)
                            aTTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAgg                                                                                                                           >  7:141329/1‑145 (MQ=255)
                             ttACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGAt                                                                                                                      <  7:80817/149‑1 (MQ=255)
                             ttACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGGGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGAt                                                                                                                      <  2:337350/149‑1 (MQ=255)
                                     tGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGATGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGc                                                                                                                                                   >  8:333481/1‑112 (MQ=255)
                                     tGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGATGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGc                                                                                                                                                   <  7:333481/112‑1 (MQ=255)
                                       tAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGCATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGAt                                                                                                              <  5:394239/147‑1 (MQ=255)
                                       tAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCa                                                                                                            <  8:35131/149‑1 (MQ=255)
                                          tttCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCa                                                                                                         >  1:207135/1‑149 (MQ=255)
                                                               ttCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTg                                                                                    <  3:370137/149‑1 (MQ=255)
                                                               ttCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTg                                                                                    <  3:156638/149‑1 (MQ=255)
                                                                 cGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTgg                                                                                   <  2:286163/148‑1 (MQ=255)
                                                                  ggTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTg                                                                                 >  5:93389/1‑149 (MQ=255)
                                                                      tGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCAc                                                                              <  4:187090/148‑1 (MQ=255)
                                                                      tGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCAc                                                                              <  7:222814/148‑1 (MQ=255)
                                                                         cTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCa                                                                                             >  2:153112/1‑130 (MQ=255)
                                                                         cTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCa                                                                                             <  1:153112/130‑1 (MQ=255)
                                                                            gTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGAc                                                                       >  3:191130/1‑149 (MQ=255)
                                                                             tGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACg                                                                      <  6:179371/149‑1 (MQ=255)
                                                                             tGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACg                                                                      <  2:107550/149‑1 (MQ=255)
                                                                              gACCAATGCACGTAATTCACGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATTGGAATGGAGCGTATTGCCCCTACGTTTGCCGAGGTAGATGGATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACgg                                                                     <  6:322260/149‑1 (MQ=255)
                                                                               aCCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACggg                                                                    <  6:347235/149‑1 (MQ=255)
                                                                                 cAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGGGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGa                                                                   >  6:325859/1‑148 (MQ=255)
                                                                                 cAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGa                                                                   >  8:180764/1‑148 (MQ=255)
                                                                                   aTGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGGGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATac                                                                <  5:361859/149‑1 (MQ=255)
                                                                                   aTGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTa                                                                                                                   <  8:17676/98‑1 (MQ=255)
                                                                                   aTGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTa                                                                                                                   >  7:17676/1‑98 (MQ=255)
                                                                                        ggTAATGCGCGAATTTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAacac                                                           <  2:239415/149‑1 (MQ=255)
                                                                                          tAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGGGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACg                                                                      <  1:270319/136‑1 (MQ=255)
                                                                                          tAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGGGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACg                                                                      >  2:270319/1‑136 (MQ=255)
                                                                                               cgcgAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCt                                                     >  2:423042/1‑148 (MQ=255)
                                                                                                  gAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGa                                                   >  8:404776/1‑147 (MQ=255)
                                                                                                  gAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGa                                                   >  5:261603/1‑147 (MQ=255)
                                                                                                           ccACGCTGCCTAAAACGCATATTGAAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGCTCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTGTGACGGGATACAACACCTGGCTGACAGGACcgcgc                                        >  6:392958/1‑149 (MQ=255)
                                                                                                              cGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACcgcgcg                                       >  8:183098/1‑147 (MQ=255)
                                                                                                                   ccTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACCGCGCGAAGCTgg                                >  3:137614/1‑149 (MQ=255)
                                                                                                                     tAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACCGCGCGAAGCTGGGc                              >  1:82421/1‑149 (MQ=255)
                                                                                                                          cGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACCGCGCGAAGCTGGGCACGtt                         >  4:306255/1‑149 (MQ=255)
                                                                                                                          cGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATAACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACCACACCTGGCAGACAGGACCGCGCGACCCTGGGCACGtt                         >  2:76909/1‑149 (MQ=255)
                                                                                                                              tatTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACCGCGCGAAGCTGGGCACGTTGATg                     <  1:423042/149‑1 (MQ=255)
                                                                                                                                  gCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACCGCGCGAAGCTGGGCACGTTGATGGTcc                 >  3:78232/1‑149 (MQ=255)
                                                                                                                                  gCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACCGCGCGAAGCTGGGCACGTTGATGGTc                  <  2:207135/148‑1 (MQ=255)
                                                                                                                                          gggAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACCGCGCGAAGCTGGGCACGTTGATGGTCCTGAAGAGt         >  7:275659/1‑149 (MQ=255)
                                                                                                                                          gggAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACCGCGCGAAGCTGGGCACGTTGATGGTCCTGAAGAGt         >  7:232887/1‑149 (MQ=255)
                                                                                                                                                 gAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACCGCGCGAAGCTGGGCACGTTGATGGTCCTGAAGAGTTTCATCt  <  4:191130/149‑1 (MQ=255)
                                                                                                                                                 gAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACCGCGCGAAGCTGGGCACGTTGATGGTCCTGAAGAGTTTCATCt  <  6:191824/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACCTGGCTGACAGGACCGCGCGAAGCTGGGCACGTTGATGGTCCTGAAGAGTTTCATCT  >  NC_000913/322916‑323209

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