Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,588,522 G→T 5.2% R193R (CGG→AGG)  PP_1393 ← transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,588,5220GT 5.2% 244.3 / 3.2 77R193R (CGG→AGG) PP_1393transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (38/35);  new base T (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

GCGACCAATCGCCGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGTTCGGTGCATCGCCGTAG  >  NC_002947/1588378‑1588670
                                                                                                                                                |                                                                                                                                                    
gCGACCAATCGCCGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAg                                                                                                                                                  <  8:38029/149‑1 (MQ=255)
 cGACCAATCGCCGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAgg                                                                                                                                                 <  5:19035/149‑1 (MQ=255)
  gACCAATCGCCGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAggg                                                                                                                                                <  5:336366/149‑1 (MQ=255)
   aCCAATCGCCGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGa                                                                                                                                               <  4:85591/149‑1 (MQ=255)
      aaTCGCCGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATgc                                                                                                                                            <  4:69868/149‑1 (MQ=255)
       accGCCGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATg                                                                                                                                             >  4:37696/3‑147 (MQ=255)
        tCGCCGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATgcag                                                                                                                                          >  5:214051/1‑149 (MQ=255)
          gCCGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCTCGCGGGATgaaca                                                                                                                                         <  1:366735/148‑5 (MQ=255)
          gCCGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCTCGCGGGATgaaca                                                                                                                                         <  4:147992/148‑5 (MQ=255)
           ccGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCCCAACCGCCAGTCGTGTTCATCCCGCGAGGGATgcagcag                                                                                                                                       >  8:291489/1‑149 (MQ=255)
            cGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATgcagca                                                                                                                                        >  6:22421/1‑147 (MQ=255)
            cGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATgcagca                                                                                                                                        <  5:22421/147‑1 (MQ=255)
                  ccACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGc                                                                                                                                <  1:115810/149‑1 (MQ=255)
                  ccACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGc                                                                                                                                >  5:336851/1‑149 (MQ=255)
                   cACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGcc                                                                                                                               >  5:56487/1‑149 (MQ=255)
                          gTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAg                                                                                                                        <  5:219715/149‑1 (MQ=255)
                           tCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAg                                                                                                                        <  6:336851/148‑1 (MQ=255)
                             cccTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTc                                                                                                                     >  8:267464/1‑149 (MQ=255)
                               cTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGATATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGACTGCACCCACCACCGCCAGTCGTGTTCATCCCGCGAGGGATGCCGCAGCCTGTGCCTGGCCAGGTCtt                                                                                                                   >  7:242161/1‑149 (MQ=255)
                               cTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCATGTCtt                                                                                                                   <  6:214051/149‑1 (MQ=255)
                               cTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCtt                                                                                                                   >  1:160256/1‑149 (MQ=255)
                                      gacatggcgatGGCCATGGCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACctg                                                                                                            >  4:228128/3‑149 (MQ=255)
                                       ccatggccatCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACctgc                                                                                                           >  2:10558/1‑149 (MQ=255)
                                             ccatCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAAc                                                                                                     <  1:10558/149‑1 (MQ=255)
                                                   ccATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCAcc                                                                                               >  7:310683/1‑149 (MQ=255)
                                                   ccATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCAcc                                                                                               >  5:69141/1‑149 (MQ=255)
                                                     aTGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACACCCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCgg                                                                                             <  6:63857/149‑1 (MQ=255)
                                                          cAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCtt                                                                                        >  3:161541/1‑149 (MQ=255)
                                                            gtgtCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAACGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCa                                                                                      >  4:349186/1‑149 (MQ=255)
                                                             tgtCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCtt                                                                                        >  4:155901/1‑146 (MQ=255)
                                                             tgtCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCtt                                                                                        <  3:155901/146‑1 (MQ=255)
                                                               tctcGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACg                                                                                    >  3:8789/1‑148 (MQ=255)
                                                                 tcGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAg                                                                                 >  2:333143/1‑149 (MQ=255)
                                                                   gAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACCACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGt                                                                                >  6:101830/1‑148 (MQ=255)
                                                                    aaaTGCGGGTAGGGTCCCTGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCACTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGa                                                                              <  8:242161/149‑1 (MQ=255)
                                                                     aaTGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGAc                                                                             <  4:161541/149‑1 (MQ=255)
                                                                      aTGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACg                                                                            <  1:215283/149‑1 (MQ=255)
                                                                      aTGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACg                                                                            <  3:154623/149‑1 (MQ=255)
                                                                        gCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGa                                                                                                                                               >  8:289788/1‑80 (MQ=255)
                                                                        gCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGa                                                                                                                                               <  7:289788/80‑1 (MQ=255)
                                                                         cTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGt                                                                          <  6:69141/148‑1 (MQ=255)
                                                                            tttGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTgcgc                                                                      <  2:180078/149‑1 (MQ=255)
                                                                               gggTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCaa                                                                                                      >  2:337901/1‑114 (MQ=255)
                                                                                 gTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGcaca                                                                   >  5:25360/1‑147 (MQ=255)
                                                                                   ccATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAggg                                                              <  8:101540/149‑1 (MQ=255)
                                                                                    cATGCGTGGCCAACGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAggg                                                              <  3:349186/149‑1 (MQ=255)
                                                                                         gTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAg                                                         <  3:228128/149‑1 (MQ=255)
                                                                                         gTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAg                                                         >  8:243873/1‑149 (MQ=255)
                                                                                            gCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAgctg                                                      >  5:237227/1‑149 (MQ=255)
                                                                                              cAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAgctgct                                                    >  1:79847/1‑149 (MQ=255)
                                                                                                  gTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATc                                                <  2:160256/149‑1 (MQ=255)
                                                                                                  gTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCAACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATc                                                <  7:267464/149‑1 (MQ=255)
                                                                                                   tCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATc                                                >  7:89999/1‑148 (MQ=255)
                                                                                                      ttCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGt                                                            >  5:76029/1‑133 (MQ=255)
                                                                                                      ttCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGt                                                            <  6:76029/133‑1 (MQ=255)
                                                                                                      ttCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAAc                                            <  7:243873/149‑1 (MQ=255)
                                                                                                      ttCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAAc                                            <  4:15633/149‑1 (MQ=255)
                                                                                                         cAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGc                                         <  7:291489/149‑1 (MQ=255)
                                                                                                            ccgccTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTg                                      >  6:264724/1‑149 (MQ=255)
                                                                                                                   cAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCt                               <  1:333143/149‑1 (MQ=255)
                                                                                                                    aaCCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTg                              >  4:325386/1‑149 (MQ=255)
                                                                                                                     aCCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTg                              >  8:239266/1‑148 (MQ=255)
                                                                                                                      ccACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTg                              >  3:130203/1‑147 (MQ=255)
                                                                                                                       cacaACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAgctgct                                                    >  5:123180/1‑124 (MQ=255)
                                                                                                                       cacaACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAgctgct                                                    <  6:123180/124‑1 (MQ=255)
                                                                                                                          aaCCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGATCATCGAACAGCTTGAGAACCTGCAGccc                        >  4:220181/1‑149 (MQ=255)
                                                                                                                          aaCCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCAAACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGccc                        >  6:257121/1‑149 (MQ=255)
                                                                                                                               ccAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGt                   >  5:30172/1‑149 (MQ=255)
                                                                                                                                  gTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGTTCg                <  3:37696/149‑1 (MQ=255)
                                                                                                                                       ctgcATCCCTCGCGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGTTCGGTGCa           >  2:166743/4‑149 (MQ=255)
                                                                                                                                       ctgcATCCCTCGCGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGTTCGGTGCa           >  6:379397/4‑149 (MQ=255)
                                                                                                                                           aTCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGTTCgg               <  5:117146/141‑1 (MQ=255)
                                                                                                                                           aTCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGTTCGGTGCATCg        <  6:30172/148‑1 (MQ=255)
                                                                                                                                           aTCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGTTCGGTGCATAGc       >  1:185589/1‑149 (MQ=255)
                                                                                                                                           aTCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGGGGTTCgg               >  6:117146/1‑141 (MQ=255)
                                                                                                                                            tCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGTTCGGTGCATCGcc      <  4:8789/149‑1 (MQ=255)
                                                                                                                                              ccGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGTTCGGTGCATCGCCg     >  4:147479/1‑148 (MQ=255)
                                                                                                                                                gcgAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGTTCGGTGCATCGCCg     >  5:75550/1‑146 (MQ=255)
                                                                                                                                                gcgAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGTTCGGTGCATCGCCGTAg  <  5:264724/149‑1 (MQ=255)
                                                                                                                                                gcgAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTCGTTCGGTGCATCGCCg     >  2:210789/1‑146 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCGACCAATCGCCGATCGCCACCCCAGTCCCCTGCGAGGCCATGGCCATCGCCATGTCCAGTGTCTCGAAATGCTGTTTGGGTCCATGCGTGGCCAGCGTCATTCCAGCCGCCTGCAACCACAACCGCCAGTCGTGTTCATCCCGCGAGGGATGCAGCAGCATGTGCCTGGCCAGGTCTTCCACCTGCTGCAACGGCACCGGCCCTTCACGCAGTGACGGTGCGCACACAGGGGTCAGCTGCTCATCGAACAGCTTGCGCACCTGCAGCCCGTGGTTCGGTGCATCGCCGTAG  >  NC_002947/1588378‑1588670

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: