Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,551,285 A→C 9.4% intergenic (+44/+82) groL → / ← pykA chaperonin GroEL/pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,551,2850AC 9.4% 177.8 / 4.7 61intergenic (+44/+82)groL/pykAchaperonin GroEL/pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base A (29/26);  new base C (2/4);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTA  >  NC_002947/1551144‑1551432
                                                                                                                                             |                                                                                                                                                    
tCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCgg                                                                                                                                              >  3:80931/1‑149 (MQ=255)
 caccacCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCggg                                                                                                                                             >  2:106922/1‑149 (MQ=255)
  accaccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCgggg                                                                                                                                            <  7:335731/149‑1 (MQ=255)
        gAAGCCATGGTTGCCGCCCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCttttt                                                                                                                                      >  5:107134/1‑149 (MQ=255)
            ccATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAtt                                                                                                                                  >  6:143832/1‑149 (MQ=255)
             cATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTg                                                                                                                                 >  8:171027/1‑149 (MQ=255)
               tGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGCGGCATGGGGGGCATGGGCGGCATGATGTAAGCCCGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCt                                                                                                                               <  6:309447/149‑1 (MQ=255)
                     ccGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCggg                                                                                                                         <  7:193828/149‑1 (MQ=255)
                           tGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCctat                                                                                                                   >  3:245991/1‑147 (MQ=255)
                             ccGGAAGACAATCAAGCTGCTGGCATGCCTGATATGGGGTGCATGCGGCGCATGGGCGGCATGATGTAAGCCAGCCTTACACCCCTGCACCATCAAAAAGCCCCGGTTCGCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCctct                                                                                                                   <  4:213142/147‑1 (MQ=255)
                                gAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgca                                                                                                              <  4:29870/149‑1 (MQ=255)
                                         ccAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGc                                                                                                     >  2:203970/1‑149 (MQ=255)
                                         ccAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGc                                                                                                     <  5:143832/149‑1 (MQ=255)
                                           aGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCgg                                                                                                                                              >  3:125074/1‑106 (MQ=255)
                                           aGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCgg                                                                                                                                              <  4:125074/106‑1 (MQ=255)
                                                  ggaaTGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                            <  7:270344/146‑1 (MQ=255)
                                                   gCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGATCTACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                            <  1:106922/149‑1 (MQ=255)
                                                      tGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCt                                                                                        >  1:354320/1‑149 (MQ=255)
                                                          tGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCggg                                                                                                                         <  4:234616/112‑1 (MQ=255)
                                                          tGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCggg                                                                                                                         >  3:234616/1‑112 (MQ=255)
                                                             tatGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtg                                                                                 <  5:353085/149‑1 (MQ=255)
                                                             tatGGGTGGCATGGGAGGCATGCGCGGCATGAAGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCCGTTCGCGACAACCGGGGCTGTTTTATAGCTCTCGGGGCCTCTTCGCACCACAAGGCTGCGCTTGCGCCTAGTGCtg                                                                                 >  6:309075/1‑149 (MQ=255)
                                                                  gtggcatgggtggcatggGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGaaa                                                                            <  8:238098/149‑1 (MQ=255)
                                                                   tggcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                            <  6:384164/132‑1 (MQ=255)
                                                                   tggcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGa                                                                                            >  5:384164/1‑132 (MQ=255)
                                                                     gcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAAc                                                                           >  4:44734/1‑147 (MQ=255)
                                                                       atgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                       <  6:107134/149‑1 (MQ=255)
                                                                       atgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                       <  3:15462/149‑1 (MQ=255)
                                                                          ggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                                     <  7:171027/148‑1 (MQ=255)
                                                                                atgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAAcc                                                              >  4:11258/1‑149 (MQ=255)
                                                                                   ggcggcATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGAc                                                           <  3:11258/149‑1 (MQ=255)
                                                                                    gcggcATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGAc                                                           >  6:359397/1‑148 (MQ=255)
                                                                                       gcATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGc                                                       <  1:343234/149‑1 (MQ=255)
                                                                                           gatgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGAGCAACCGACCAGCGGAt                                                   >  1:271962/1‑149 (MQ=255)
                                                                                            atgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGa                                                                    >  6:311451/1‑131 (MQ=255)
                                                                                            atgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGa                                                                    <  5:311451/131‑1 (MQ=255)
                                                                                               tAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCAcc                                               >  4:31387/1‑149 (MQ=255)
                                                                                                aaGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATc                                                                  >  8:171353/1‑129 (MQ=255)
                                                                                                aaGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATc                                                                  <  7:171353/129‑1 (MQ=255)
                                                                                                 agccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGGCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGa                                             <  2:354320/149‑1 (MQ=255)
                                                                                                  gccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGAc                                            >  4:123518/1‑149 (MQ=255)
                                                                                                  gccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGAc                                            >  8:202925/1‑149 (MQ=255)
                                                                                                  gccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGAc                                            >  5:137097/1‑149 (MQ=255)
                                                                                                   ccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGCCCAGCGGATCACCGACg                                           >  8:160462/1‑149 (MQ=255)
                                                                                                        cTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAg                                      <  3:31387/149‑1 (MQ=255)
                                                                                                               cccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGAt                                                   >  3:60592/1‑129 (MQ=255)
                                                                                                               cccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGAt                                                   <  4:60592/129‑1 (MQ=255)
                                                                                                               cccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTc                               <  6:137097/149‑1 (MQ=255)
                                                                                                                 cTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCATGTGCTGTGAAACGAGAGGATCAACCGACCAGCGGATCACCGAAGTGCAGGAACTTCAt                             >  8:57035/1‑149 (MQ=255)
                                                                                                                     aCCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCa                         >  3:354551/1‑149 (MQ=255)
                                                                                                                     aCCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGAGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGCCCAGCGGATCACCGACGTGCAGGATCTTCATGCCa                         >  6:55255/1‑149 (MQ=255)
                                                                                                                        aTCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTg                      <  7:160462/149‑1 (MQ=255)
                                                                                                                            aaaaGCCCCGGATCTCGCGACCGGGGGCTTTTTTATTGATGTCGGGGCCTGTTCCCAGCACAAGGGTGCGCTTGAGCGTAGTGCTGTGAAACGAGTGGATCACCGGACGAGCGGATCACCGCCGTGCAGGATGTTCATGCCATTGGtcc                  >  8:122467/1‑147 (MQ=255)
                                                                                                                               aGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc               >  6:139177/1‑149 (MQ=255)
                                                                                                                                 ccccGGTTCGCGAGAACCGGGGCTTTTTTAATGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAcc              >  1:36965/1‑148 (MQ=255)
                                                                                                                                  cccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGa            >  1:169530/1‑149 (MQ=255)
                                                                                                                                    cGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg          <  5:359397/149‑1 (MQ=255)
                                                                                                                                      gTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgt        <  7:202925/149‑1 (MQ=255)
                                                                                                                                       ttCTCGCGAACCGGGGCTTTTTGATGGTGCAGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgtgt      <  6:374338/149‑1 (MQ=38)
                                                                                                                                          gcgATAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTgg    <  4:80931/149‑1 (MQ=255)
                                                                                                                                          gcgAGTACCGGGGCTTTTTTATTGCTGTCGGGGCGTCTTCGGAGGACAAGGCTGCGCTTGAGCTTTGTGCTGTGAAACGAGTGGATCAAGCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGtgtg     <  3:15958/148‑1 (MQ=255)
                                                                                                                                            gagaACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTa  <  1:152937/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
TCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTA  >  NC_002947/1551144‑1551432

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