Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,865,286 T→A 8.3% D120V (GAC→GTC)  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,2860TA 8.3% 202.4 / 7.5 72D120V (GAC→GTC) PP_3414methyl‑accepting chemotaxis transducer/sensory box protein
Reads supporting (aligned to +/- strand):  ref base T (29/37);  new base A (2/4);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

CTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCTCAC  >  NC_002947/3865138‑3865430
                                                                                                                                                    |                                                                                                                                                
cTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGAACCAACCGCTGGGTGATGTCCGACGCCACCTTGATTACGCGTGTGACTTTCCGGTGCTCGt                                                                                                                                                  >  3:384305/1‑149 (MQ=255)
  tGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTc                                                                                                                                                 >  5:296454/1‑148 (MQ=255)
   gATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTc                                                                                                                                                 <  8:206262/147‑1 (MQ=255)
   gATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTc                                                                                                                                                 >  5:10125/1‑147 (MQ=255)
   gATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTc                                                                                                                                                 >  6:244338/1‑147 (MQ=255)
        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                                                                                                                                     <  5:139215/149‑2 (MQ=255)
             ccTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCCGTGCTGGACGAGCACCGGCaa                                                                                                                                     <  1:74126/149‑2 (MQ=255)
             ccTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCCGTGCTGGACGAGCACCGGCaa                                                                                                                                     <  6:268651/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                                                                                                              >  3:99243/1‑149 (MQ=255)
                                    cGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCt                                                                                                              >  2:272545/1‑149 (MQ=255)
                                    cGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCt                                                                                                              >  7:8201/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                                 <  6:288692/149‑1 (MQ=255)
                                                                                                                 ccACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTgg                                 <  1:283199/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                                                        <  3:150013/121‑1 (MQ=255)
                                                                                                                      ttGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTtgccctgcc                                                        >  4:150013/1‑121 (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)
                                                                                                                                                    |                                                                                                                                                
CTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGTCCGGCTCAC  >  NC_002947/3865138‑3865430

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