Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,865,278 G→C 8.2% R123G (CGG→GGG)  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,2780GC 8.2% 220.7 / 6.5 72R123G (CGG→GGG) PP_3414methyl‑accepting chemotaxis transducer/sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (30/35);  new base C (2/4);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

GTTGGCCTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGT  >  NC_002947/3865132‑3865421
                                                                                                                                                  |                                                                                                                                               
gTTGGCCTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGt                                                                                                                                               >  3:280890/1‑149 (MQ=255)
 ttGGCCTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTg                                                                                                                                              <  2:62266/149‑1 (MQ=255)
    gCCTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCg                                                                                                                                                 <  5:161480/143‑1 (MQ=255)
    gCCTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCg                                                                                                                                                 >  6:161480/1‑143 (MQ=255)
      cTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGAACCAACCGCTGGGTGATGTCCGACGCCACCTTGATTACGCGTGTGACTTTCCGGTGCTCGt                                                                                                                                         >  3:384305/1‑149 (MQ=255)
        tGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTc                                                                                                                                        >  5:296454/1‑148 (MQ=255)
         gATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTc                                                                                                                                        <  8:206262/147‑1 (MQ=255)
         gATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTc                                                                                                                                        >  6:244338/1‑147 (MQ=255)
         gATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTc                                                                                                                                        >  5:10125/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                                                                                                                            <  1:74126/149‑2 (MQ=255)
                   ccTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCCGTGCTGGACGAGCACCGGCaa                                                                                                                            <  6:268651/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                                                                                                     >  2:272545/1‑149 (MQ=255)
                                          cGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCt                                                                                                     >  7:8201/1‑149 (MQ=255)
                                          cGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCt                                                                                                     >  3:99243/1‑149 (MQ=255)
                                            ccATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTcc                                                                                                               <  1:181716/137‑1 (MQ=255)
                                            ccATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTcc                                                                                                               >  2:181716/1‑137 (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                                                                                                                     <  5:292285/97‑1 (MQ=255)
                                                                              gCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTg                                                                                                                     >  6:292285/1‑97 (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)
                                                                                                                                                  |                                                                                                                                               
GTTGGCCTTGATCACCCGGCCTTGGGGGGTGAACTCGATCACCGCCATGGAGCGGCTGATGGCCTTGAGCAGGCTTTCGCTCTCGTGCTCCTGCACCACCCGCTGGGTGATGTCCGACGCCACCTTGATTACGCTGGTGACTTGCCGGTGCTCGTCCAGCACGGGCATGTAGCTGGCTTCCAGCCAGACCTCGCGCCCCGCCTTGTCGATACGCTCGAACGTACCGCTGATAGCCTTGCCCTGCCCCAACTCGCGCCACAGTTGCTGGTACTCGGCGCTTTTTGCATAGT  >  NC_002947/3865132‑3865421

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: