Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,822,188 G→T 6.4% intergenic (+30/+28) PP_mr22 → / ← fdxA PrrB_RsmZ ncRNA/ferredoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,822,1880GT 6.4% 191.2 / 3.7 63intergenic (+30/+28)PP_mr22/fdxAPrrB_RsmZ ncRNA/ferredoxin 1
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base T (2/2);  total (28/35)
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.44e-01

ACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAG  >  NC_002947/1822046‑1822329
                                                                                                                                               |                                                                                                                                             
acaGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCGCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGc                                                                                                                                         >  5:337574/1‑149 (MQ=255)
        tCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATTATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCTTCGGGCCttttt                                                                                                                                   <  4:227258/147‑1 (MQ=255)
        tCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTc                                                                                                                                  <  2:1267/148‑1 (MQ=255)
          aGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCTTCGGGCCttttttgtt                                                                                                                               <  7:239123/149‑5 (MQ=255)
           ggATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGt                                                                                                                              <  6:192566/149‑1 (MQ=255)
            gATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTc                                                                                                                             >  1:76736/1‑149 (MQ=255)
            gATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTc                                                                                                                             <  2:76736/149‑1 (MQ=255)
             aTGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGCACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAATGGCCCGAAGGGCCTTTTTCGCGTCt                                                                                                                            <  4:331792/149‑1 (MQ=255)
               gAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTgg                                                                                                                          >  1:3371/1‑149 (MQ=255)
                aaGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTggg                                                                                                                         <  5:245231/149‑1 (MQ=255)
                aaGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGCAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAACGGCCCGACGGGCCTTTTTCGCGTCTgg                                                                                                                          <  2:191087/148‑1 (MQ=255)
                 aGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACACAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTggg                                                                                                                         <  1:53015/149‑1 (MQ=255)
                        ggACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGCACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCa                                                                                                                 <  2:20108/149‑1 (MQ=255)
                         gACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGGAGGGCCTTTTTCGCGTCTGGGCGCTATCAg                                                                                                                >  7:172192/1‑149 (MQ=255)
                             tCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGGAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCt                                                                                                            <  7:113610/149‑1 (MQ=255)
                              cGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTc                                                                                                           >  2:182296/1‑149 (MQ=255)
                               gCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTgg                                                                                                                          <  4:195525/133‑1 (MQ=255)
                                 aagaagCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAg                                                                                                        >  8:296232/1‑149 (MQ=255)
                                    aagCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTc                                                                                                     <  2:344987/149‑1 (MQ=255)
                                     agCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCg                                                                                                    <  2:22205/149‑1 (MQ=255)
                                         aTTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGCCGGGTCAAACCGACCCTTTTTTTGCCAGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGAGTCTGGGCGCTATCACCGCTCCAGGTCGGCaa                                                                                                >  8:226829/1‑149 (MQ=255)
                                              atcaGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTg                                                                                           >  8:172573/1‑149 (MQ=255)
                                               tcaGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGc                                                                                          >  3:131773/1‑149 (MQ=255)
                                                 aGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCa                                                                                        <  8:257988/149‑1 (MQ=255)
                                                        tgTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTg                                                                                  <  3:256015/148‑1 (MQ=255)
                                                               gACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCa                                                                          >  2:80030/1‑149 (MQ=255)
                                                                  aaGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAttc                                                                       >  7:349961/1‑149 (MQ=255)
                                                                    gCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCGCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAttctt                                                                     <  6:337574/149‑1 (MQ=255)
                                                                    gCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAttctt                                                                     >  5:297343/1‑149 (MQ=255)
                                                                     cAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAttctt                                                                     >  5:34958/1‑148 (MQ=255)
                                                                      aGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCa                                                                   <  3:138300/149‑1 (MQ=255)
                                                                       gggACTACGGAAAAAACTGTGGGCGGGGCAAACCGCCCCGTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGcc                                                                                <  6:225856/135‑1 (MQ=255)
                                                                         gACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCg                                                                >  7:22550/1‑149 (MQ=255)
                                                                            tACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTcc                                                             >  1:293703/1‑149 (MQ=255)
                                                                                gAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGAAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCgg                                                           >  5:276271/1‑147 (MQ=255)
                                                                                gAAAAAAATGTGGGCGGGTCAAACAGCCCCTTTTTTTTCCCGTAGGAAAACAAAAAAGCCCCGAAGGGCCTTTTTCGCGTATGGGAGCTATCAGAGCTCCAGGTCGGCCATCATGCCAAGCTTGACCTACCATTATTAAGCGtaccg                                                           >  4:361729/1‑143 (MQ=255)
                                                                                      aaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGt                                                   >  8:110161/1‑149 (MQ=255)
                                                                                        tgtgGGCGGGTCAATCCGCCCCTTTTTTTGCCCGTAGGAAAAGAAAAAAGGCCCGTAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTcc                                                 <  7:226829/149‑1 (MQ=255)
                                                                                            ggCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTAc                                             >  8:200816/1‑149 (MQ=255)
                                                                                                  tCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCgg                                       <  6:34958/149‑1 (MQ=255)
                                                                                                    aaaCCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGGAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGt                                      <  8:172192/148‑1 (MQ=255)
                                                                                                    aaaCCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGt                                      >  2:68196/1‑148 (MQ=255)
                                                                                                    aaaCCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTg                                     <  7:296232/149‑1 (MQ=255)
                                                                                                    aaaCCGCCCCTTTTTTTGCCCCAAGCAAAACAAAAAAGGCCCGAAGGGCATTTTTCGAGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTTCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGAAGGGCGTCCTTACGCTCGGt                                      >  8:238417/1‑148 (MQ=255)
                                                                                                       ccGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGata                                  <  2:3371/149‑1 (MQ=255)
                                                                                                        cGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGatat                                 >  8:5794/1‑149 (MQ=255)
                                                                                                         gCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATAtt                                <  7:110161/149‑1 (MQ=255)
                                                                                                          ccccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGGAGGGCGTCCTTACGCTCGGTGATATTa                               <  8:342715/149‑1 (MQ=255)
                                                                                                           cccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAg                              <  8:349961/149‑1 (MQ=255)
                                                                                                            ccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAgg                             >  8:304682/1‑149 (MQ=255)
                                                                                                              tttttttGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGc                            >  7:40657/1‑148 (MQ=255)
                                                                                                                tttttGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCa                          <  4:131773/148‑1 (MQ=255)
                                                                                                                 ttttGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAg                         <  2:293703/148‑1 (MQ=255)
                                                                                                                      cccGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTcc                   <  5:241911/149‑1 (MQ=255)
                                                                                                                      cccGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTcc                   <  8:40657/149‑1 (MQ=255)
                                                                                                                          tAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCa               <  1:80030/149‑1 (MQ=255)
                                                                                                                               aaaCAAAAAAGGCCCGAAGGGACTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCACTCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCg          <  8:284415/149‑1 (MQ=255)
                                                                                                                                  gcgaaaaaGGCCCTTCGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGc        >  6:94971/4‑149 (MQ=255)
                                                                                                                                  cAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCg       >  3:37709/1‑149 (MQ=255)
                                                                                                                                    gaaaaaGGCCCTTCGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGt      >  7:197908/2‑149 (MQ=255)
                                                                                                                                        aaGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAg  <  1:182296/149‑1 (MQ=255)
                                                                                                                                              cGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTg                                     >  5:116047/1‑108 (MQ=255)
                                                                                                                                              cGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTg                                     <  6:116047/108‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
ACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAG  >  NC_002947/1822046‑1822329

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