Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,472,545 C→T 8.0% V92M (GTG→ATG)  rsfS ← ribosomal silencing factor RsfS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,472,5450CT 8.0% 178.5 / 4.4 63V92M (GTG→ATG) rsfSribosomal silencing factor RsfS
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base T (3/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GTCTTTGAGCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGTTGGAGGTACCGGTGGCAAT  >  NC_002947/5472408‑5472689
                                                                                                                                         |                                                                                                                                                
gTCTTTGAGCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGAcgtcgt                                                                                                                                       >  3:173217/1‑149 (MQ=255)
 tCTTTGAGCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGtt                                                                                                                                      <  3:246957/149‑1 (MQ=255)
   tttGAGCTTGTCTGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGa                                                                                                                                             <  4:63256/140‑1 (MQ=255)
   tttGAGCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGa                                                                                                                                             >  3:63256/1‑140 (MQ=255)
   tttGAGCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGtt                                                                                                                                      >  6:363769/1‑147 (MQ=255)
     tGAGCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAg                                                                                                                                   <  2:106601/148‑1 (MQ=255)
      gAGCTTGTCGGAGTACTTGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAgg                                                                                                                                  <  2:354577/148‑1 (MQ=255)
            gTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCcag                                                                                                                             <  1:267642/147‑1 (MQ=255)
                     cTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCAAGATGACGTCGTTCAGGTCCAGCAGCACCCAGt                                                                                                                  >  4:216675/1‑149 (MQ=255)
                        gTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCGTGCACATGATGAcc                                                                                                                                           <  6:167921/121‑2 (MQ=255)
                        gTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCGTGCACATGATGAcc                                                                                                                                           >  5:167921/1‑120 (MQ=255)
                        gTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAgg                                                                                                                                  >  8:117170/1‑130 (MQ=255)
                        gTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCAcc                                                                                                                      <  6:14865/142‑1 (MQ=255)
                        gTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCAcc                                                                                                                      >  5:14865/1‑142 (MQ=255)
                            gCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTc                                                                                                                                     >  1:206118/1‑123 (MQ=255)
                            gCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTc                                                                                                                                     <  2:206118/123‑1 (MQ=255)
                                              tgCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCAc                                                                                          >  6:347973/1‑148 (MQ=255)
                                                 ttGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAg                                                                                      >  1:179521/1‑149 (MQ=255)
                                                   gCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGt                                                                                     <  3:20422/148‑1 (MQ=255)
                                                    ccATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTgg                                                                                   >  5:156880/1‑149 (MQ=255)
                                                      aTCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCt                                                                                 <  4:173217/149‑1 (MQ=255)
                                                                cGGCGTTGCTCGGCACCCTGCCACAGACGCTCACGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCtt                                                                       <  5:198857/149‑1 (MQ=255)
                                                                    tttGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGcc                                                                   >  8:190391/1‑149 (MQ=255)
                                                                       gCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCtt                                                                 <  6:353809/148‑1 (MQ=255)
                                                                       gCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCtt                                                                 >  1:350953/1‑148 (MQ=255)
                                                                        cTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGa                                                               <  2:350953/149‑1 (MQ=255)
                                                                         tCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGcc                                                                   >  8:84567/1‑144 (MQ=255)
                                                                         tCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGcc                                                                   <  7:84567/144‑1 (MQ=255)
                                                                                ccTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCa                                                       >  6:10113/1‑149 (MQ=255)
                                                                                ccTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCa                                                       >  1:229370/1‑149 (MQ=255)
                                                                                 cTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCAc                                                      >  6:368683/1‑149 (MQ=255)
                                                                                    ccACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCtt                                                                 >  5:175716/1‑135 (MQ=255)
                                                                                    ccACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACg                                                     <  2:275653/147‑1 (MQ=255)
                                                                                    ccACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTAGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCtt                                                                 <  6:175716/135‑1 (MQ=255)
                                                                                      acaGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCtt                                                                 <  2:2230/133‑1 (MQ=255)
                                                                                        aGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCtt                                               >  7:32659/1‑149 (MQ=255)
                                                                                                      gTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATcgcg                                  <  4:93538/148‑1 (MQ=255)
                                                                                                      gTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGt                                 <  8:32659/149‑1 (MQ=255)
                                                                                                       tAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGtt                                <  2:229370/149‑1 (MQ=255)
                                                                                                        aGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGtt                                >  8:266068/1‑148 (MQ=255)
                                                                                                            cttgCGGGAAGCGGGTGTCATCATGTGCACGATGACTTCGTTCATGTCCAGCAGCACCCAGTCTCTTTCGCCATTGCCTTCTTCACCCAGTTGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTTATCt                           <  3:216675/146‑1 (MQ=255)
                                                                                                            cTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCt                           <  5:363769/149‑1 (MQ=255)
                                                                                                             tGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTg                          >  7:120192/1‑149 (MQ=255)
                                                                                                               gCGGGCAGCGGCGGTCCCCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACgg                                                    <  1:298086/121‑1 (MQ=255)
                                                                                                               gCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACgg                                                    >  2:298086/1‑121 (MQ=255)
                                                                                                                 gggCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCgg                      >  4:351116/1‑149 (MQ=255)
                                                                                                                  ggCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCAt                                      <  2:335899/132‑1 (MQ=255)
                                                                                                                  ggCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCAt                                      <  2:330700/132‑1 (MQ=255)
                                                                                                                  ggCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCAt                                      >  1:330700/1‑132 (MQ=255)
                                                                                                                  ggCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCAt                                      >  1:335899/1‑132 (MQ=255)
                                                                                                                  ggCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGt                     <  5:10113/149‑1 (MQ=255)
                                                                                                                   gCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGtt                    >  6:197245/1‑149 (MQ=255)
                                                                                                                    cAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGTTg                   <  5:347973/149‑1 (MQ=255)
                                                                                                                     gaacgacgTCATCGTGCACATGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGc                        >  2:57359/8‑143 (MQ=255)
                                                                                                                     gaacgacgTCATCGTGCACATGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGc                        <  1:57359/136‑1 (MQ=255)
                                                                                                                       cggcggTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGTTGGAg                <  5:368683/149‑1 (MQ=255)
                                                                                                                        ggcggTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGTTGGAgg               <  6:156880/149‑1 (MQ=255)
                                                                                                                         gacgTCATCGTGCACATGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGTTGGAgg               >  5:117888/4‑148 (MQ=255)
                                                                                                                           ggTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGTTGGAGGTAc            <  7:190391/149‑1 (MQ=255)
                                                                                                                                tcatGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATc                            >  3:35131/1‑128 (MQ=255)
                                                                                                                                tcatGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATc                            <  4:35131/128‑1 (MQ=255)
                                                                                                                                    gtgCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGTTGGAGGTACCGGTGGCaa   >  4:346390/1‑149 (MQ=255)
                                                                                                                                     tgCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGTTGGAGGTACCGGTGGCAAt  >  6:285420/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                
GTCTTTGAGCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGTTGGAGGTACCGGTGGCAAT  >  NC_002947/5472408‑5472689

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: