Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,865,288 C→G 9.8% L119L (CTG→CTC PP_3414 ← methyl‑accepting chemotaxis transducer/sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,865,2880CG 9.8% 191.7 / 9.4 70L119L (CTG→CTCPP_3414methyl‑accepting chemotaxis transducer/sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (27/36);  new base G (2/5);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCTCACAGAA  >  NC_002947/3865146‑3865434
                                                                                                                                              |                                                                                                                                                  
cgcgtccTTGGGGGGTGAACTCGATCACCGCAATGTAGCTGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACTCCACCTTGATTACGCTGGTGACTTGCCTTTGCTCGTCCAGCAc                                                                                                                                               <  7:27683/143‑1 (MQ=255)
cccGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAg                                                                                                                                                  <  5:320809/145‑1 (MQ=255)
cccGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAg                                                                                                                                                  >  6:320809/1‑145 (MQ=255)
 ccGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACgg                                                                                                                                             <  6:296454/149‑1 (MQ=255)
   ggCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCCGTGCTGGACGAGCACcgg                                                                                                                                            >  6:139215/1‑148 (MQ=255)
     ccTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCCGTGCTGGACGAGCACCGGCaa                                                                                                                                         <  6:268651/149‑2 (MQ=255)
     ccTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCCGTGCTGGACGAGCACCGGCaa                                                                                                                                         <  1:74126/149‑2 (MQ=255)
     ccTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCCGTGCTGGACGAGCACCGGCaa                                                                                                                                         <  5:139215/149‑2 (MQ=255)
      cTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATg                                                                                                                                        <  2:320395/149‑1 (MQ=255)
         ggggggTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTa                                                                                                                                      <  2:24899/148‑1 (MQ=255)
           ggggTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCCCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCt                                                                                                                                   <  4:54807/149‑1 (MQ=255)
            gggTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTg                                                                                                                                  >  1:374851/1‑149 (MQ=255)
                 aaCTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTc                                                                                                                             <  2:374851/149‑1 (MQ=255)
                   cTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTcca                                                                                                                           <  3:133694/149‑1 (MQ=255)
                      gATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTccagc                                                                                                                         >  8:203436/1‑148 (MQ=255)
                       aTCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTccagcca                                                                                                                       <  2:171383/149‑1 (MQ=255)
                            cGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCt                                                                                                                  >  7:8201/1‑149 (MQ=255)
                            cGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCt                                                                                                                  >  2:272545/1‑149 (MQ=255)
                            cGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCt                                                                                                                  >  3:99243/1‑149 (MQ=255)
                              ccATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTcc                                                                                                                            >  2:181716/1‑137 (MQ=255)
                              ccATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTcc                                                                                                                            <  1:181716/137‑1 (MQ=255)
                              ccATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTcg                                                                                                                >  2:259748/1‑149 (MQ=255)
                               cATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTcgc                                                                                                               <  2:90/149‑1 (MQ=255)
                                      cGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGc                                                                                                        <  4:55982/149‑1 (MQ=255)
                                         cTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCtt                                                                                                     <  1:272545/149‑1 (MQ=255)
                                          tGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTg                                                                                                    >  1:213570/1‑149 (MQ=255)
                                            aTGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTc                                                                                                  <  1:259748/149‑1 (MQ=255)
                                             tGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCg                                                                                                 <  4:280890/149‑1 (MQ=255)
                                               gCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGAt                                                                                               <  1:81735/149‑1 (MQ=255)
                                                        aGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGaa                                                                                      >  2:94908/1‑149 (MQ=255)
                                                         ggCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAAc                                                                                     <  6:267397/149‑1 (MQ=255)
                                                          gCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACg                                                                                    <  8:8201/149‑1 (MQ=255)
                                                           cTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGt                                                                                                   >  8:248646/1‑133 (MQ=255)
                                                            tttCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTgg                                                                                                                                 <  4:228707/102‑1 (MQ=255)
                                                            tttCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTgg                                                                                                                                 >  3:228707/1‑102 (MQ=255)
                                                             ttCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTAc                                                                                 >  7:94227/1‑149 (MQ=255)
                                                                gCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTg                                                                                                                                  >  6:292285/1‑97 (MQ=255)
                                                                gCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTg                                                                                                                                  <  5:292285/97‑1 (MQ=255)
                                                                  tctcGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCt                                                                             >  2:283199/1‑148 (MQ=255)
                                                                   ctcGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGa                                                                           >  5:319415/1‑149 (MQ=255)
                                                                          tCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCtt                                                                    >  4:198239/1‑149 (MQ=255)
                                                                               caccacCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTtgcc                                                                 <  2:213570/147‑1 (MQ=255)
                                                                                 ccaccCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTtgccctgc                                                             <  6:10125/149‑1 (MQ=255)
                                                                                  caccCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTtgccctgcc                                                            >  8:178248/1‑149 (MQ=255)
                                                                                                tGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCcac                                              <  4:99011/149‑1 (MQ=255)
                                                                                                     gACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTg                                         <  7:203436/149‑1 (MQ=255)
                                                                                                         ccACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCGAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTgg                                     <  5:72611/149‑1 (MQ=255)
                                                                                                         ccACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTgg                                     <  1:283199/149‑1 (MQ=255)
                                                                                                         ccACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTgg                                     <  6:288692/149‑1 (MQ=255)
                                                                                                           aCCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGt                                    >  8:175852/1‑148 (MQ=255)
                                                                                                            ccTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTAc                                  >  6:326849/1‑149 (MQ=255)
                                                                                                             cTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACt                                 >  6:410930/1‑149 (MQ=255)
                                                                                                              ttGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTtgccctgcc                                                            >  4:150013/1‑121 (MQ=255)
                                                                                                              ttGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTtgccctgcc                                                            <  3:150013/121‑1 (MQ=255)
                                                                                                                gATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCg                               <  8:170714/148‑1 (MQ=255)
                                                                                                                     cGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCtt                         <  6:319415/149‑1 (MQ=255)
                                                                                                                          gTGACTTGCCGGTGCTCGTCCAGCACGTGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTg                     <  4:99243/148‑1 (MQ=255)
                                                                                                                          gTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACTGTTGCTGGTACTCGGCGCttttt                      <  8:316732/147‑1 (MQ=255)
                                                                                                                          gTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCttttt                      <  1:94908/147‑1 (MQ=255)
                                                                                                                          gTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGc                    <  7:175852/149‑1 (MQ=255)
                                                                                                                           tGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCa                   >  4:141761/1‑149 (MQ=255)
                                                                                                                           ctacatGCCCGTGCTGGACGAGCACCGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCttttt                      <  7:192755/141‑1 (MQ=255)
                                                                                                                           ctacatGCCCGTGCTGGACGAGCACCGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCttttt                      >  8:192755/6‑146 (MQ=255)
                                                                                                                             aCTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATa                 >  4:31765/1‑149 (MQ=255)
                                                                                                                               ttGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGt               >  1:63205/1‑149 (MQ=255)
                                                                                                                                        cTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCTca       <  3:198239/148‑1 (MQ=255)
                                                                                                                                           gTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCTca       >  7:197876/1‑145 (MQ=255)
                                                                                                                                           gTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCTca       <  8:197876/145‑1 (MQ=255)
                                                                                                                                            tCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCTcac      >  6:330816/1‑145 (MQ=255)
                                                                                                                                            tCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCTcac      <  5:330816/145‑1 (MQ=255)
                                                                                                                                             ccAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCTCACAGaa  >  4:219208/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCTCACAGAA  >  NC_002947/3865146‑3865434

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: