Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,938,664 T→G 8.0% S409A (TCG→GCG)  rbbA → ribosome‑associated ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,938,6640TG 8.0% 186.1 / 5.5 61S409A (TCG→GCG) rbbAribosome‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base G (3/2);  total (33/28)
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.00e+00

AACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTCCACCCA  >  NC_002947/5938520‑5938812
                                                                                                                                                |                                                                                                                                                    
aaCCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGACAGGCGCt                                                                                                                                                  >  4:319346/1‑149 (MQ=255)
                cccGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCt                                                                                                                                   <  7:49833/148‑1 (MQ=255)
                       gTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGcc                                                                                                                           >  2:158264/1‑149 (MQ=255)
                       gTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGcc                                                                                                                           >  1:54094/1‑149 (MQ=255)
                           gACTTGCCCAAAGCCGAGAGCGCCCAGCACATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCgg                                                                                                                        >  8:36730/1‑148 (MQ=255)
                                   cAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAATCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGc                                                                                                                                 >  5:95242/1‑131 (MQ=255)
                                   cAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAATCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGc                                                                                                                                 <  6:95242/131‑1 (MQ=255)
                                    aaaGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGc                                                                                                                            <  6:270571/135‑1 (MQ=255)
                                    aaaGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGa                                                                                                               <  6:362104/148‑1 (MQ=255)
                                    aaaGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGa                                                                                                               >  5:362104/1‑148 (MQ=255)
                                        ccGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTc                                                                                                          <  2:54094/149‑1 (MQ=255)
                                               cgcCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGACAGGCGCTGGcgcaa                                                                                                                                           <  1:314320/109‑5 (MQ=255)
                                               cgcCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGACAGGCGCGGGcgcaa                                                                                                                                           >  2:314320/1‑105 (MQ=255)
                                                gcCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAAcc                                                                                                  <  1:158264/149‑1 (MQ=255)
                                                  ccAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGACAGGCGCTGGcg                                                                                                                                              <  3:41434/103‑2 (MQ=255)
                                                  ccAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGACAGGCGCTGGcg                                                                                                                                              >  4:41434/1‑102 (MQ=255)
                                                   cAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTg                                                                                                                <  4:194185/132‑1 (MQ=255)
                                                   cAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTg                                                                                                                >  3:194185/1‑132 (MQ=255)
                                                           cGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCg                                                                                       >  5:305710/1‑149 (MQ=255)
                                                             aTGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTc                                                                                     <  7:36730/149‑1 (MQ=255)
                                                               gAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCggagatag                                                                                   >  3:283722/1‑143 (MQ=255)
                                                               gAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCg                                                                                          <  4:283722/142‑1 (MQ=255)
                                                                 gCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGAcc                                                                                 >  1:152852/1‑149 (MQ=255)
                                                                  cTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGa                                                                                                                       <  5:167518/110‑1 (MQ=255)
                                                                  cTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGa                                                                                                                       >  6:167518/1‑110 (MQ=255)
                                                                     aTCGAGCGCTTCGCCCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCccgc                                                                             >  3:137587/1‑149 (MQ=255)
                                                                        gAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCcgccgc                                                                          >  4:178797/1‑149 (MQ=255)
                                                                          gcgcTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGccgccc                                                                        >  5:305634/1‑149 (MQ=255)
                                                                           cgcTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGccgcccg                                                                       >  5:283298/1‑149 (MQ=255)
                                                                               tCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGcgac                                                                   >  1:347236/1‑149 (MQ=255)
                                                                                 gACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGcgacga                                                                 <  3:178797/149‑1 (MQ=255)
                                                                                  aCCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCgacgac                                                                <  6:283298/149‑1 (MQ=255)
                                                                                       gccgcTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCt                                                            <  6:305634/148‑1 (MQ=255)
                                                                                        ccgcTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTgg                                                          <  4:137587/149‑1 (MQ=255)
                                                                                              tCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGACTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCctg                                                    >  1:253558/1‑149 (MQ=255)
                                                                                                gCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCctgc                                                   <  2:152852/148‑1 (MQ=255)
                                                                                                      cAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAg                                             >  8:303273/1‑148 (MQ=255)
                                                                                                       aaCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGa                                                                                   <  5:206802/109‑1 (MQ=255)
                                                                                                       aaCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGa                                                                                   >  6:206802/1‑109 (MQ=255)
                                                                                                       aaCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGaa                                           >  1:57024/1‑149 (MQ=255)
                                                                                                        aCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAAc                                          >  5:322012/1‑149 (MQ=255)
                                                                                                          cGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTg                                        >  8:154319/1‑149 (MQ=255)
                                                                                                            tCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTc                                      <  6:305710/149‑1 (MQ=255)
                                                                                                            tCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTc                                      >  1:138497/1‑149 (MQ=255)
                                                                                                              tGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTccc                                    >  7:31712/1‑149 (MQ=255)
                                                                                                                gcgcCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTcc                                     <  6:322012/146‑1 (MQ=255)
                                                                                                                  gcCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGa                                <  2:347236/149‑1 (MQ=255)
                                                                                                                        ccccTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGAc                                                                                               <  7:114030/80‑1 (MQ=255)
                                                                                                                        ccccTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGAc                                                                                               >  8:114030/1‑80 (MQ=255)
                                                                                                                        ccccTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAggg                          >  7:138795/1‑149 (MQ=255)
                                                                                                                        ccccTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAgg                           <  7:303273/148‑1 (MQ=255)
                                                                                                                         cccTTGGTTTGCGCCAGCGCCTGTCGCTGGGAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGc                         >  1:27855/1‑149 (MQ=255)
                                                                                                                          ccTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTg                                                                                                                >  1:111424/1‑61 (MQ=255)
                                                                                                                          ccTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTg                                                                                                                <  2:111424/61‑1 (MQ=255)
                                                                                                                          ccTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGACCCGCGAGCAGGGCg                        >  2:249910/1‑149 (MQ=255)
                                                                                                                           cTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGt                       <  8:31712/149‑1 (MQ=255)
                                                                                                                              ggTTTGCGCCAGCGCCTGTCGCTGGGAGTGGCTGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGa                     <  2:27855/148‑1 (MQ=255)
                                                                                                                               gTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGAc                    <  8:138795/148‑1 (MQ=255)
                                                                                                                                     gcCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTc             <  7:154319/149‑1 (MQ=255)
                                                                                                                                       cAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTcc            >  1:309217/1‑148 (MQ=255)
                                                                                                                                             cTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTccac     >  5:100672/1‑149 (MQ=255)
                                                                                                                                                tCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTccaccca  <  6:100672/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTCCACCCA  >  NC_002947/5938520‑5938812

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: