Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,110,656 A→T 6.5% L285Q (CTA→CAA)  PP_5361 ← protein P47K

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,110,6560AT 6.5% 185.9 / 2.6 61L285Q (CTA→CAA) PP_5361protein P47K
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base T (2/2);  total (30/31)
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.80e-01

TTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTGGTACATGTTCACCGCGCAGCTCCTGCAACCAGCCTGGT  >  NC_002947/6110509‑6110802
                                                                                                                                                   |                                                                                                                                                  
ttATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAg                                                                                                                                                   >  6:55666/1‑149 (MQ=255)
   tCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCt                                                                                                                                                >  1:135464/1‑149 (MQ=255)
     cacaTCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTc                                                                                                                                              >  8:225707/1‑149 (MQ=255)
                gCCACTGCTCGCGTGGCACGAAGCGCCACCAGCTGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAg                                                                                                                                   <  2:57228/149‑1 (MQ=255)
                gCCACTGCTCGCGTGGCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCa                                                                                                                                    <  7:189169/148‑1 (MQ=255)
                  cACTGCTCGCGTGGCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTc                                                                                                                                 <  4:55621/149‑1 (MQ=255)
                   aCTGCTCGCGTGGCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTcc                                                                                                                                >  8:144871/1‑149 (MQ=255)
                       cTCGCGTGGCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTtgt                                                                                                                            <  8:46655/149‑1 (MQ=255)
                       cTCGCGTGGCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTtgt                                                                                                                            <  3:68776/149‑1 (MQ=255)
                       cTCGCGTGGCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTtgt                                                                                                                            >  8:180163/1‑149 (MQ=255)
                        tCGCGTGGCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTtgtg                                                                                                                           >  8:68748/1‑149 (MQ=255)
                              ggCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATa                                                                                                                       <  8:262073/147‑1 (MQ=255)
                               gCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAg                                                                                                                    <  7:150609/149‑1 (MQ=255)
                               gCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAg                                                                                                                    <  7:68748/149‑1 (MQ=255)
                                 aCGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTc                                                                                                                  <  7:144871/149‑1 (MQ=255)
                                   gAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGt                                                                                                                <  7:225707/149‑1 (MQ=255)
                                       cgcCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTaa                                                                                                              >  1:97903/1‑147 (MQ=255)
                                        gcgagcaGTGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGtttt                                                                                                                                        <  6:125788/115‑1 (MQ=255)
                                        gcgagcaGTGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGtttt                                                                                                                                        >  5:125788/6‑120 (MQ=255)
                                          caccaGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAgcg                                                                                                         >  8:108117/1‑149 (MQ=255)
                                           accaGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAgcgc                                                                                                        >  5:71485/1‑149 (MQ=255)
                                            ccaGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATTTTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCATTCCGTTTTCCCAGTCCTTCTGGATAAACTCGTAAAAGCGCt                                                                                                       >  7:109223/1‑149 (MQ=255)
                                                  gccagccaGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGAt                                                                                                 >  5:190748/1‑149 (MQ=255)
                                                   ccagccagACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATg                                                                                                <  2:97903/149‑1 (MQ=255)
                                                      gccagACCATGGCGCATCATGCCGCCCACTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAg                                                                                             >  8:201925/1‑149 (MQ=255)
                                                           aCCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGc                                                                                        <  2:135464/149‑1 (MQ=255)
                                                                 gcgcATCATGCCGCCCGCTTGCGACCAGCGGCCAGCTTCCTGATCTTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCGGTGGATGAAGCGGGCGTCGa                                                                                   <  4:163772/148‑1 (MQ=255)
                                                                  cgcATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTtgt                                                                                                                            >  8:12784/1‑106 (MQ=255)
                                                                  cgcATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTtgt                                                                                                                            <  7:12784/106‑1 (MQ=255)
                                                                   gcATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGcc                                                                                <  6:190748/149‑1 (MQ=255)
                                                                    catcatGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTtgt                                                                                                                            >  7:196800/1‑104 (MQ=255)
                                                                    catcatGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTtgt                                                                                                                            <  8:196800/104‑1 (MQ=255)
                                                                    catcatGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCt                                                                               <  6:71485/149‑1 (MQ=255)
                                                                     atcatGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTg                                                                              >  2:130189/1‑149 (MQ=255)
                                                                      tcatGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTTCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGc                                                                             <  1:246787/149‑1 (MQ=255)
                                                                            cgcccacTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTc                                                                       <  7:201925/149‑1 (MQ=255)
                                                                               ccgcTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATg                                                                    <  5:55666/149‑1 (MQ=255)
                                                                                     gCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCcgc                                                               >  5:191484/1‑148 (MQ=255)
                                                                                             cTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTa                                                      <  1:130189/149‑1 (MQ=255)
                                                                                                 cAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTActct                                                  >  1:10197/1‑149 (MQ=255)
                                                                                                  aGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCtt                                                 <  2:192399/149‑1 (MQ=255)
                                                                                                   cattCTTGCTATTTCCAGGCCAACCATAAAAAACCTTTCGAACCCAGTACCCTTCCGTTTTAGCCTTTCTTGTTGTTAAATTCGTAAAACCGCTGTGGATTAAGCGGGCGTCGTGCCTGCGAGGTCATGGCCGCGATACCGTACTCTTc                                                <  6:53467/147‑1 (MQ=255)
                                                                                                         tGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTActc                                                   <  2:236239/140‑1 (MQ=255)
                                                                                                         tGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTActc                                                   >  1:236239/1‑140 (MQ=255)
                                                                                                         tGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTc                                          >  3:300317/1‑149 (MQ=255)
                                                                                                          gATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCTTActctgc                                                >  5:53467/1‑140 (MQ=255)
                                                                                                            tatTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTgg                                       <  6:191484/149‑1 (MQ=255)
                                                                                                            tatTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTgg                                       <  7:108117/149‑1 (MQ=255)
                                                                                                             atTTGCTGGCCAACCATAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTTGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTgg                                       <  8:109223/148‑1 (MQ=255)
                                                                                                                 gCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTGGTACAt                                  >  7:182591/1‑149 (MQ=255)
                                                                                                                     gCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTGGTACATGTTc                              >  6:83188/1‑149 (MQ=255)
                                                                                                                       cAACCAGAAAAACCCTTTGGAACGCAGTTGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGGGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTc                                                <  5:191807/129‑1 (MQ=255)
                                                                                                                       cAACCAGAAAAACCCTTTGGAACGCAGTTGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTc                                                <  5:196321/129‑1 (MQ=255)
                                                                                                                       cAACCAGAAAAACCCTTTGGAACGCAGTTGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTc                                                >  6:196321/1‑129 (MQ=255)
                                                                                                                       cAACCAGAAAAACCCTTTGGAACGCAGTTGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTc                                                >  6:191807/1‑129 (MQ=255)
                                                                                                                        aaCCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGATGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTGGTACATGTTCAcc                           >  5:98189/1‑149 (MQ=255)
                                                                                                                         aCCAGAAAAACCCTTTGGAACGCGGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTGGTACATGTTCACcg                          >  2:70827/1‑149 (MQ=255)
                                                                                                                                                aGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTGGTACATGTTCACCGCGCAGCTCCTGCAACCAGCCTgg   >  3:235978/1‑149 (MQ=255)
                                                                                                                                                aGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTGGTACATGTTCACCGCGCAGCTCCTGCAACCAGCCTgg   >  5:268450/1‑149 (MQ=255)
                                                                                                                                                aGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTGGTACATGTTCACCGCGCAGCTCCTGCAACCAGCCTgg   >  5:81000/1‑149 (MQ=255)
                                                                                                                                                 gTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTGGTACATGTTCACCGCGCAGCTCCTGCAACCAGCCTGGt  <  7:180163/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCACCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGATGAAGCGGGCGTCGAGCCTGCCAGGTCATGGCCGCGATACCGTACTCTTCGGTTTCTGGTACATGTTCACCGCGCAGCTCCTGCAACCAGCCTGGT  >  NC_002947/6110509‑6110802

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: ‑