Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,257,689 A→C 5.9% V38G (GTG→GGG)  PP_4634 ← alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,257,6890AC 5.9% 268.5 / 7.2 85V38G (GTG→GGG) PP_4634alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base A (40/40);  new base C (3/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

GAGGCCTGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGC  >  NC_002947/5257542‑5257835
                                                                                                                                                    |                                                                                                                                                  
gaggCCTGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCAc                                                                                                                                                   >  1:185819/1‑149 (MQ=255)
   gCCTGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAg                                                                                                                                                >  6:42151/1‑149 (MQ=255)
      tGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCgg                                                                                                                                             >  4:147336/1‑149 (MQ=255)
      tGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCgg                                                                                                                                             >  8:79959/1‑149 (MQ=255)
        gAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGgc                                                                                                                                           >  3:57679/1‑149 (MQ=255)
             cACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTgg                                                                                                                                      <  3:201439/149‑1 (MQ=255)
             cACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTgg                                                                                                                                      <  2:299639/149‑1 (MQ=255)
             cACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTgg                                                                                                                                      >  2:117174/1‑149 (MQ=255)
              aCGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTAGCAGCGTTGCTGATCACCAGCGGGCGCgggg                                                                                                                                     >  8:247632/1‑149 (MQ=255)
                 cAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCCCCAGCGCCCGCTGGTGAt                                                                                                                                  <  2:218934/149‑1 (MQ=255)
                  aGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCCCCAGCGCCCGCTGGTGATc                                                                                                                                 <  6:123277/149‑1 (MQ=255)
                   gCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATTGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCa                                                                                                                                <  3:251290/149‑1 (MQ=255)
                     gAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTCGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAgg                                                                                                                              <  7:247632/149‑1 (MQ=255)
                        gccggccgGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGCAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGtt                                                                                                                            <  1:71120/149‑1 (MQ=255)
                           ggccggGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGc                                                                                                                        <  2:188736/149‑1 (MQ=255)
                           ggccggGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGc                                                                                                                        <  7:79959/149‑1 (MQ=255)
                           ggccggGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGc                                                                                                                        >  7:95323/1‑149 (MQ=255)
                              cggGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTg                                                                                                                     >  2:346090/1‑149 (MQ=255)
                                 gATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGgcg                                                                                                                  >  7:66028/1‑149 (MQ=255)
                                   ttCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGcc                                                                                                                <  5:42151/149‑1 (MQ=255)
                                      gCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCcag                                                                                                              >  7:317712/1‑148 (MQ=255)
                                      gCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCcag                                                                                                              >  5:163476/1‑148 (MQ=255)
                                          aTGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGTGCGTAATAGCTGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGaa                                                                                                         <  3:336466/149‑1 (MQ=255)
                                          aTGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGaa                                                                                                         <  6:282562/149‑1 (MQ=255)
                                              cACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAgccg                                                                                                     >  2:109802/1‑149 (MQ=255)
                                                     gTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAgg                                                                                               <  3:147336/148‑1 (MQ=255)
                                                     gTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAgg                                                                                               <  5:107101/148‑1 (MQ=255)
                                                     gTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAgg                                                                                               <  2:229679/148‑1 (MQ=255)
                                                      tCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCt                                                                                             <  1:109802/149‑1 (MQ=255)
                                                               aGCACATCGCGGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTc                                                                                    >  7:325967/1‑149 (MQ=255)
                                                               aGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGACCACCAGCGGGCGCTGGGGATGAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTc                                                                                    <  1:241734/149‑1 (MQ=255)
                                                                 cacaTCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCgg                                                                                  <  3:202485/149‑1 (MQ=255)
                                                                 cacaTCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGACAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCgg                                                                                  >  6:324513/1‑149 (MQ=255)
                                                                 cacaTCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGACAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCgg                                                                                  >  6:329461/1‑149 (MQ=255)
                                                                  acaTCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCgg                                                                                  >  6:2267/1‑148 (MQ=255)
                                                                  acaTCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCgg                                                                                  <  5:2267/148‑1 (MQ=255)
                                                                  acaTCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGc                                                                                 <  6:87051/149‑1 (MQ=255)
                                                                    aTCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCcg                                                                               >  2:269737/1‑149 (MQ=255)
                                                                         aGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCgggg                                                                          >  5:335515/1‑149 (MQ=255)
                                                                            cGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTgg                                                                       <  6:163476/149‑1 (MQ=255)
                                                                                 gCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGaaa                                                                  >  5:281166/1‑149 (MQ=255)
                                                                                  cGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCTTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAAt                                                                 <  8:317712/149‑1 (MQ=255)
                                                                                  cGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAAt                                                                 >  3:78997/1‑149 (MQ=255)
                                                                                  cGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAAt                                                                 >  3:213553/1‑149 (MQ=255)
                                                                                   gAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATg                                                                >  3:236077/1‑149 (MQ=255)
                                                                                     aCAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATggg                                                              >  4:75186/1‑149 (MQ=255)
                                                                                     aCAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATggg                                                              <  4:78997/149‑1 (MQ=255)
                                                                                     aCAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATggg                                                              <  2:185819/149‑1 (MQ=255)
                                                                                     aCAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATggg                                                              <  1:346090/149‑1 (MQ=255)
                                                                                      cAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATgggg                                                             <  6:335515/149‑1 (MQ=255)
                                                                                        ggTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATgggg                                                             >  3:184085/1‑147 (MQ=255)
                                                                                         gTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAg                                                          >  8:254056/1‑149 (MQ=255)
                                                                                          tAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGc                                                         <  6:281166/149‑1 (MQ=255)
                                                                                             tCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCgg                                                                            <  6:296295/127‑1 (MQ=255)
                                                                                             tCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGCCTGTATCCGTCGGCCGCgg                                                                            >  5:296295/1‑127 (MQ=255)
                                                                                               ggCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGtt                                                    <  1:269737/149‑1 (MQ=255)
                                                                                               ggCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGtt                                                    <  1:117174/149‑1 (MQ=255)
                                                                                                gCAAAGCGGTCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATg                                                                <  8:304049/136‑1 (MQ=255)
                                                                                                gCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATg                                                                >  7:304049/1‑136 (MQ=255)
                                                                                                gCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGtt                                                    <  4:169199/148‑1 (MQ=255)
                                                                                                gCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGtt                                                    >  3:169199/1‑148 (MQ=255)
                                                                                                  aaaGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGa                                                 >  7:257510/1‑149 (MQ=255)
                                                                                                              aTAGCGGCAGCGTACCG‑AAGTTGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTc                                      <  4:57679/148‑1 (MQ=255)
                                                                                                                 gCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATgg                                  <  8:95323/149‑1 (MQ=255)
                                                                                                                  cGGTAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGgc                                 >  1:256394/1‑149 (MQ=255)
                                                                                                                  cGGTAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGgc                                 <  2:256394/149‑1 (MQ=255)
                                                                                                                      aGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCAt                             >  5:300336/1‑149 (MQ=255)
                                                                                                                       gCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATg                            >  5:194521/1‑149 (MQ=255)
                                                                                                                         gTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGGGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTc                          <  6:213977/149‑1 (MQ=255)
                                                                                                                             cG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCtgt                      >  3:322432/1‑149 (MQ=255)
                                                                                                                                  gTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCtgttgtt                  <  4:184085/149‑1 (MQ=255)
                                                                                                                                     gCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATg                                                                >  5:259578/1‑100 (MQ=255)
                                                                                                                                     gCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATg                                                                <  6:259578/100‑1 (MQ=255)
                                                                                                                                        gCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTgcgc            <  5:324513/149‑1 (MQ=255)
                                                                                                                                        gCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTgcgc            <  5:329461/149‑1 (MQ=255)
                                                                                                                                          gTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGcc           >  8:69904/1‑148 (MQ=255)
                                                                                                                                          caacctgatCCCCAGCGCCCGCTGGTGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGcc           >  4:85023/6‑148 (MQ=255)
                                                                                                                                           aacctgatCCCCAGCGCCCGCTGGTGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCa          >  3:229407/5‑148 (MQ=255)
                                                                                                                                             gatgatCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCaaaa       >  6:192371/1‑149 (MQ=255)
                                                                                                                                               tgatCCCCAGCGCCCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAAc      >  2:153521/1‑148 (MQ=255)
                                                                                                                                                gatCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCaaaa       <  4:315853/146‑1 (MQ=255)
                                                                                                                                                gatCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCaaaa       >  3:315853/1‑146 (MQ=255)
                                                                                                                                                gatCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAAc      >  7:66118/1‑147 (MQ=255)
                                                                                                                                                 atCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAg   <  2:299589/149‑1 (MQ=255)
                                                                                                                                                   cACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGc  >  2:43639/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAGGCCTGGAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGC  >  NC_002947/5257542‑5257835

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