Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 534,067 T→G 5.6% intergenic (+8/‑37) PP_t04 → / → secE tRNA‑Trp/protein translocase subunit SecE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947534,0670TG 5.6% 227.6 / 3.8 71intergenic (+8/‑37)PP_t04/secEtRNA‑Trp/protein translocase subunit SecE
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base G (2/2);  total (35/36)
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.21e-01

GAATAAGTAGTTGATTTACTTGATCAGGCCGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGTGTTGTGGGTAATCAATATTACTCCGGCTCTC  >  NC_002947/533924‑534215
                                                                                                                                               |                                                                                                                                                    
gAATAAGAAGATGATTTACTTGATCAGGCCGGTAAAATGGCCGCCCTGATACAGCGTTACACGACAGTAGATCAATAGGCAGAGCGACGGACTCCAAAACCGTAGGTACGGGGATCGATTCCCTCCTCACCTGCCATTTCACCTCGgtg                                                                                                                                                 >  8:121683/1‑149 (MQ=255)
     agtagtTTACTTACTTGATCAGGCCGGTATAATGGCCGGCCTGATACATCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACTGTCTCCAAAACCGTAGTTTGGTGGTTCGATTCCCTCCTTACCTGCCATTTCACCTCGTTGTGAtt                                                                                                                                            <  8:362590/149‑2 (MQ=255)
     agtagtTGATTTACTTGATCAGGCCGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATg                                                                                                                                            <  6:166011/149‑1 (MQ=255)
            gATTTACTTGATCAGGCGGGTATAATGGCCGGCCTGATACAGCGTTTCAGGGGAGTAGCTCAATTGGCAGAGCGACGGTCTCCTAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGAGCTGCCATTTCACCTCGGTGTGATGGCTttc                                                                                                                                      <  7:121683/148‑1 (MQ=255)
              tttACTTGATCAGGCCGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTttcttc                                                                                                                                   >  2:102688/1‑149 (MQ=255)
              tttACTTGATCAGGCCGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTttcttc                                                                                                                                   <  2:13266/149‑1 (MQ=255)
                 aCTTGATCAGGCCGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTa                                                                                                                                >  6:69276/1‑149 (MQ=255)
                 aCTTGATCAGGCCGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCGAGGTGAAATGGca                                                                                                                                         >  1:363064/1‑139 (MQ=255)
                 aCTTGATCAGGCCGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCGAGGTGAAATGGca                                                                                                                                         <  2:363064/140‑2 (MQ=255)
                   ttGATCAGGCCGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTGTTCTTa                                                                                                                                >  3:81997/1‑147 (MQ=255)
                       tCAGGCCGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGAt                                                                                                                          >  7:358153/1‑149 (MQ=255)
                           gCCGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATCACACCGAGGTGAAATGGcagg                                                                                                                                       >  4:34158/1‑129 (MQ=255)
                           gCCGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATCACACCGAGGTGAAATGGcagg                                                                                                                                       <  3:34158/132‑4 (MQ=255)
                             ggttataATGGCGGGCCTGATACAGCGTTACAGGTCAGTAGGTCAATTGGCAGAGCGACGTTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATgg                                                                                                                                           <  6:214653/123‑1 (MQ=255)
                             cGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACATCGGTGTGATgg                                                                                                                                           >  5:214653/1‑126 (MQ=255)
                                      ggccggccTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGAccc                                                                                                           >  6:59981/1‑149 (MQ=255)
                                         cggccTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGAccccc                                                                                                         >  4:223093/1‑148 (MQ=255)
                                              tGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACt                                                                                                   >  8:86769/1‑149 (MQ=255)
                                                 tACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGaa                                                                                                <  3:272087/149‑1 (MQ=255)
                                                 tACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGa                                                                                                 <  5:6058/148‑1 (MQ=255)
                                                    aGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGcc                                                                                             >  5:213406/1‑149 (MQ=255)
                                                       gTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCaa                                                                                          <  5:146796/149‑1 (MQ=255)
                                                        ttACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTa                                                                                                                                <  2:80698/110‑1 (MQ=255)
                                                        ttACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTa                                                                                                                                >  1:80698/1‑110 (MQ=255)
                                                              gTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCGCAAGAATggc                                                                                   >  5:386043/1‑149 (MQ=255)
                                                                  gTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCg                                                                                                                  <  8:247903/114‑1 (MQ=255)
                                                                  gTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCg                                                                                                                  >  7:247903/1‑114 (MQ=255)
                                                                  gTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGa                                                                                                 <  2:168896/131‑1 (MQ=255)
                                                                  gTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGa                                                                                                 >  1:168896/1‑131 (MQ=255)
                                                                  gTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGttt                                                                               <  3:188779/149‑1 (MQ=255)
                                                                   tAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGtttt                                                                              <  1:150545/149‑1 (MQ=255)
                                                                        cAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCt                                                                         >  3:171802/1‑149 (MQ=255)
                                                                         aaTTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTg                                                                        >  4:361424/1‑149 (MQ=255)
                                                                          aTTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGt                                                                       >  7:336568/1‑149 (MQ=255)
                                                                                agagCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCGTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTCTTGATCTGTTCAAGt                                                                 >  8:50981/1‑149 (MQ=255)
                                                                                agagCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGt                                                                 <  7:86769/149‑1 (MQ=255)
                                                                                    cGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGtggct                                                             <  3:361424/149‑1 (MQ=255)
                                                                                    cGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGtggct                                                             <  3:223093/149‑1 (MQ=255)
                                                                                      aCGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGtggctgg                                                           <  1:102688/149‑1 (MQ=255)
                                                                                      aCGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGtggctgg                                                           >  4:346598/1‑149 (MQ=255)
                                                                                         gTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGgctggctg                                                        <  5:59981/149‑1 (MQ=255)
                                                                                            tcCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCaa                                                                                          <  7:282181/112‑1 (MQ=255)
                                                                                            tcCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGaa                                                                                       <  3:188297/115‑1 (MQ=255)
                                                                                            tcCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGaa                                                                                       >  4:188297/1‑115 (MQ=255)
                                                                                            tcCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCATCGATGACCCCCAAAACTGAAGCCCaa                                                                                          >  8:282181/1‑112 (MQ=255)
                                                                                             ccAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGt                                                    >  4:241112/1‑149 (MQ=255)
                                                                                              ccaaaCCGTAGGTTGGGGGGTCGAGTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTg                                                   <  8:358153/147‑1 (MQ=255)
                                                                                              cAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTg                                                   >  1:102470/1‑149 (MQ=255)
                                                                                               aaaaCCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTggctggct                                                         >  8:318820/1‑142 (MQ=255)
                                                                                               aaaaCCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTggctggct                                                         <  7:318820/142‑1 (MQ=255)
                                                                                                   ccGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGGTTTGCTCTGTTCAAGTGGCTGGATGTAGTGGCttt                                              >  6:65397/1‑149 (MQ=255)
                                                                                                    cGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGGTCAAGTGGCTTGCTGTATTGGCTTtg                                             >  2:365631/1‑149 (MQ=255)
                                                                                                      tAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTtggt                                           <  4:171802/149‑1 (MQ=255)
                                                                                                        ggTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTtggtg                                          <  3:346598/148‑1 (MQ=255)
                                                                                                        ggTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTtggtg                                          <  5:65397/148‑1 (MQ=255)
                                                                                                        ggTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTtggtg                                          <  8:336568/148‑1 (MQ=255)
                                                                                                             ggggTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGgttgt                                     <  3:241112/148‑1 (MQ=255)
                                                                                                              gggTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGttgttg                                   <  5:69276/149‑1 (MQ=255)
                                                                                                                gTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGgt                                 >  6:147671/1‑149 (MQ=255)
                                                                                                                    gATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGtgttg                             >  3:303564/1‑149 (MQ=255)
                                                                                                                    gATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGGTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGtgtt                              <  6:213406/148‑1 (MQ=255)
                                                                                                                     aTTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGAt                                                                                                                <  1:201753/65‑1 (MQ=255)
                                                                                                                     aTTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGAt                                                                                                                >  2:201753/1‑65 (MQ=255)
                                                                                                                         cctcctGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGTGTTGTggg                         <  2:102470/148‑1 (MQ=255)
                                                                                                                          ctcctGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGTGTTGTGGGt                        >  5:42515/1‑148 (MQ=255)
                                                                                                                          ctcctGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGTGTTGTGGGt                        >  4:313899/1‑148 (MQ=255)
                                                                                                                             ctGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCGTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTCTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGTGTTGTGGGTAAt                     <  7:50981/148‑1 (MQ=255)
                                                                                                                                      cATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGTGTTGTgg                          >  6:16014/1‑134 (MQ=255)
                                                                                                                                      cATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGTGTTGTgg                          <  5:16014/134‑1 (MQ=255)
                                                                                                                                          tCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGTGTTGTGGGTAATCAATATTACTCCg        <  5:147671/148‑1 (MQ=255)
                                                                                                                                             ccTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGTGTTGTGGGTAATCAATATTACTCCGGctc    >  7:113009/1‑149 (MQ=255)
                                                                                                                                               tCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTCAAGCCAAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGTGTTGTGGGTAATCAAAATTACTCCGGCtctc  >  6:84723/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GAATAAGTAGTTGATTTACTTGATCAGGCCGGTATAATGGCCGGCCTGATACAGCGTTACAGGTCAGTAGCTCAATTGGCAGAGCGACGGTCTCCAAAACCGTAGGTTGGGGGTTCGATTCCCTCCTGACCTGCCATTTCACCTCGGTGTGATGGCTTTCTTCTTACAGGATCCTCCTCGATGACCCCCAAAACTGAAGCCCAAGAATCGCGTTTTGATCTGTTCAAGTGGCTGGCTGTAGTGGCTTTGGTGGTTGTTGGTGTTGTGGGTAATCAATATTACTCCGGCTCTC  >  NC_002947/533924‑534215

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: