Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,540,058 A→G 5.3% N249D (AAT→GAT)  PP_3127 → exopolysaccharide transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,540,0580AG 5.3% 223.3 / 2.4 75N249D (AAT→GAT) PP_3127exopolysaccharide transport protein
Reads supporting (aligned to +/- strand):  ref base A (32/39);  new base G (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

CCACGCCAATGCCTTGCCCGAGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAGGACGCCCACCCGGATGTCGTGGCGGTCAAGCGC  >  NC_002947/3539916‑3540204
                                                                                                                                               |                                                                                                                                                  
ccACGCCAATGCCTTGCCCGAGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCg                                                                                                                                              <  1:317658/149‑1 (MQ=255)
ccACGCCAATGCCTTGCCCGAGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCg                                                                                                                                              <  5:263872/149‑1 (MQ=255)
   cGCCAATGCCTTGCCCGAGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCg                                                                                                                                              >  1:368615/1‑146 (MQ=255)
   cGCCAATGCCTTGCCCGAGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCg                                                                                                                                              <  2:368615/146‑1 (MQ=255)
      cAATGCCTTGCCCGAGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTg                                                                                                                                        >  7:74225/1‑149 (MQ=255)
       aaTGCCTTGCCCGAGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGc                                                                                                                                       >  8:144270/1‑149 (MQ=255)
            cTTGCCCGAGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGccacc                                                                                                                                   >  4:48181/1‑148 (MQ=255)
                  cGAGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTg                                                                                                                            <  6:155794/149‑1 (MQ=255)
                   gAGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCGATGCCGGCCTTGCCACCAAGGtt                                                                                                                             <  8:3292/147‑1 (MQ=255)
                   gAGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCGATGCCGGCCTTGCCACCAAGGtt                                                                                                                             >  7:3292/1‑147 (MQ=255)
                    aGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGcca                                                                                                                                     <  1:170629/138‑1 (MQ=255)
                    aGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGcca                                                                                                                                     >  2:170629/1‑138 (MQ=255)
                      cACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGa                                                                                                                        <  3:48181/149‑1 (MQ=255)
                        ccAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGGTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCtt                                                                                                                                         >  8:325365/1‑130 (MQ=255)
                        ccAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGGTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCtt                                                                                                                                         <  7:325365/130‑1 (MQ=255)
                        ccAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCGATGCCGGCCTTGccac                                                                                                                                    <  1:38068/135‑1 (MQ=255)
                        ccAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCGATGCAGGCCTTGccac                                                                                                                                    >  2:38068/1‑135 (MQ=255)
                        ccAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAgg                                                                                                                      <  8:249758/149‑1 (MQ=255)
                        ccAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAgg                                                                                                                      >  8:23150/1‑149 (MQ=255)
                        ccAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAgg                                                                                                                      >  8:178915/1‑149 (MQ=255)
                        ccAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAgg                                                                                                                      <  3:924/149‑1 (MQ=255)
                              ccTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCAccc                                                                                                                >  6:322008/1‑149 (MQ=255)
                              ccTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCAccc                                                                                                                >  7:281462/1‑149 (MQ=255)
                                  cgcATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGcc                                                                                                            <  1:352951/149‑1 (MQ=255)
                                     aTGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACGCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCaa                                                                                                          >  1:263469/1‑149 (MQ=255)
                                          cATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTg                                                                                                    >  3:309196/1‑149 (MQ=255)
                                           aTGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGc                                                                                                   <  7:144270/149‑1 (MQ=255)
                                               tGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGCCCAACCCGTGCCGAc                                                                                               >  2:55163/1‑149 (MQ=255)
                                               tGTCGCGGGCGGAGCTGGAATTACGTGAGGTCGACCGTGACTACAAGGCTGCAAAGGAAGAGCTGAGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCAACCAAGGTTGCCGAGGGATCCCGGCCAACCAGTGCCGAc                                                                                               >  2:60633/1‑149 (MQ=255)
                                                gTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGc                                                                                                   >  3:358627/1‑144 (MQ=255)
                                                gTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGc                                                                                                   <  4:358627/144‑1 (MQ=255)
                                                   gcgggcggAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGc                                                                                           <  7:38418/149‑1 (MQ=255)
                                                     gggcggAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGccc                                                                                         <  7:178915/149‑1 (MQ=255)
                                                           aGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGAc                                                                                   <  3:246343/149‑1 (MQ=255)
                                                               ggAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAgg                                                                                                                      <  6:25314/110‑1 (MQ=255)
                                                                gAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAgg                                                                                                                               >  8:162208/1‑100 (MQ=255)
                                                                gAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAgg                                                                                                                               <  7:162208/100‑1 (MQ=255)
                                                                gAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGc                                                                               <  1:100449/148‑1 (MQ=255)
                                                                  aTTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCaa                                                                            <  6:181966/149‑1 (MQ=255)
                                                                     ccGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGc                                                                          <  7:23150/148‑1 (MQ=255)
                                                                      cGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCt                                                                        >  8:292457/1‑149 (MQ=255)
                                                                      cGTGAGGGCGACCGGGACTACAAGGCTGCCCCGGAACACCTGCGCTACCTGGAACTGCAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCt                                                                        <  6:375657/149‑1 (MQ=255)
                                                                          aGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGaaa                                                                    <  8:74225/149‑1 (MQ=255)
                                                                          aGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGaaa                                                                    <  4:309196/149‑1 (MQ=255)
                                                                               gACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCg                                                                <  5:25253/148‑1 (MQ=255)
                                                                                aCCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGaa                                                              >  7:305205/1‑149 (MQ=255)
                                                                                   cGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATa                                                            <  2:263469/149‑1 (MQ=255)
                                                                                       aCTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATAcgcg                                                        <  8:305205/149‑1 (MQ=255)
                                                                                           cAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGaaa                                                                    <  5:86752/133‑1 (MQ=255)
                                                                                           cAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGaaa                                                                    >  6:86752/1‑133 (MQ=255)
                                                                                           cAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCaa                                                     <  5:322008/148‑1 (MQ=255)
                                                                                              ggCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAActgc                                                 >  7:283921/1‑149 (MQ=255)
                                                                                                cTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAActgctg                                               <  7:292457/149‑1 (MQ=255)
                                                                                                 tGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGa                                              >  5:320873/1‑149 (MQ=255)
                                                                                                   cccAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGAcc                                            <  8:281462/149‑1 (MQ=255)
                                                                                                    ccAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGAccc                                           >  4:160270/1‑149 (MQ=255)
                                                                                                     cAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCg                                          >  1:349504/1‑149 (MQ=255)
                                                                                                        gAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTa                                       >  7:39491/1‑149 (MQ=255)
                                                                                                             gCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGaaa                                                                    >  2:172780/1‑115 (MQ=255)
                                                                                                             gCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGaaa                                                                    <  1:172780/115‑1 (MQ=255)
                                                                                                             gCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAg                                   >  6:119946/1‑148 (MQ=255)
                                                                                                                      ccTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAActgc                                                 <  3:170311/125‑1 (MQ=255)
                                                                                                                      ccTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAActgc                                                 >  4:170311/1‑125 (MQ=255)
                                                                                                                         ggAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAGGAc                                >  3:354191/1‑139 (MQ=255)
                                                                                                                         ggAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAGGAc                                <  4:354191/139‑1 (MQ=255)
                                                                                                                            aCTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGTCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCTACTGCTGACCCGCTACAAGGACGCCCACCCGGATg                   <  1:60633/149‑1 (MQ=255)
                                                                                                                            aCTGGAGCTGGCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCTACTGCTGACCCGCTACAAGGACGCCCACCCGGATg                   <  1:55163/149‑1 (MQ=255)
                                                                                                                                gAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAGGACGCCCACCCGGATGTCGt               <  8:283921/149‑1 (MQ=255)
                                                                                                                                   cTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAGGACGCCCACCCGGATGTCGTGGc            >  3:51363/1‑149 (MQ=255)
                                                                                                                                   cTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAGGACGCCCACCCGGATGTCGTGGc            <  8:39491/149‑1 (MQ=255)
                                                                                                                                      tCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAGGACGCCCACCCGGATGTCGTGGCgg          <  6:320873/148‑1 (MQ=255)
                                                                                                                                        tGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAGGACGCCCACCCGGATGTCGTGGCGGTc        >  7:31427/1‑148 (MQ=255)
                                                                                                                                        tGCGGCCAATGCCGGCCTTGCCACCAAGCTTGCCGAGGGCACACGGCCAACCCGTGCCGACCTGCCCGAAGACATGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAGGACGCCCACCCCGATGTCTTGGAGGTc        >  4:264061/1‑148 (MQ=255)
                                                                                                                                         gCGGCCACTGCCGGCCCTGCCACCAAGGGGGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATCCGCGCAACTGCTGACCCGCTACAAGGACGCCCACCCGCATGTCGTGGCGGTCaa      <  8:31427/149‑1 (MQ=255)
                                                                                                                                         gCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAGGACGCCCACCCGGATGTCGTGGCGGTc        >  2:166840/1‑147 (MQ=255)
                                                                                                                                               aaTGCCGGCCTTGCAACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAGGACGCCCACCCGGATGTCGTGGCGGTCAAgcgc  >  3:58604/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CCACGCCAATGCCTTGCCCGAGCACCAGACCCTGCGCATGAACATGCTGTCGCGGGCGGAGCTGGAATTCCGTGAGGTCGACCGTGACTACAAGGCTGCCCAGGAAGAGCTGCGCTACCTGGAACTGGAGCTGTCTGCGGCCAATGCCGGCCTTGCCACCAAGGTTGCCGAGGGCACCCGGCCAACCCGTGCCGACCTGCCCGAAGACCTGCCAAGCCTGAAAGCCGAATACGCGCAACTGCTGACCCGCTACAAGGACGCCCACCCGGATGTCGTGGCGGTCAAGCGC  >  NC_002947/3539916‑3540204

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: