Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,472,541 T→C 7.9% H93R (CAC→CGC) ‡ rsfS ← ribosomal silencing factor RsfS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,472,5410TC 7.9% 172.6 / 4.6 63H93R (CAC→CGC) ‡rsfSribosomal silencing factor RsfS
Reads supporting (aligned to +/- strand):  ref base T (27/31);  new base C (3/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

TCCTTATTCGCGGTCTTTGAGCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGTTGGAGGTACCGGTGGCAAT  >  NC_002947/5472396‑5472689
                                                                                                                                                 |                                                                                                                                                    
tCCTTATTCGCGGTCTTTGAGCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCa                                                                                                                                                   >  2:267642/1‑149 (MQ=255)
 ccTTATTCGCGGTCTTTGAGCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCa                                                                                                                                                   <  5:176554/148‑1 (MQ=255)
            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)
                    gCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCAtgtg                                                                                                                                                     >  4:124637/1‑127 (MQ=255)
                    gCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCAtgtg                                                                                                                                                     <  3:124637/127‑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                                                                                                                      >  5:14865/1‑142 (MQ=255)
                                    gTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCAcc                                                                                                                      <  6:14865/142‑1 (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                                                                   <  7:84567/144‑1 (MQ=255)
                                                                                     tCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGcc                                                                   >  8:84567/1‑144 (MQ=255)
                                                                                            ccTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCa                                                       >  1:229370/1‑149 (MQ=255)
                                                                                            ccTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCa                                                       >  6:10113/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                                      >  1:330700/1‑132 (MQ=255)
                                                                                                                              ggCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCAt                                      <  2:330700/132‑1 (MQ=255)
                                                                                                                              ggCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCAt                                      <  2:335899/132‑1 (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                        <  1:57359/136‑1 (MQ=255)
                                                                                                                                 gaacgacgTCATCGTGCACATGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGc                        >  2:57359/8‑143 (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)
                                                                                                                                                 |                                                                                                                                                    
TCCTTATTCGCGGTCTTTGAGCTTGTCGGAGTACTCGTGGGCATGGTCCGGGCTGTGGTGCTTGCCATCGGCGGCACGGCTTTGCTCGGCACCCTGCCACAGACGCTCAAGGTCGTAGAACTGGCGGGCAGCGGCGGTCATCATGTGCACGATGACGTCGTTCAGGTCCAGCAGCACCCAGTCGCTGTCGCCCTTGCCTTCTTCACCCAGTGGCTGGGCGCCCTTGGCCTTGACTGCTTCACGGACCTTTTCGGCCATCGCGTTGATCTGGCGGTTGGAGGTACCGGTGGCAAT  >  NC_002947/5472396‑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: