Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,257,697 G→C 5.8% R35R (CGC→CGG)  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,6970GC 5.8% 280.5 / 6.5 86R35R (CGC→CGG) PP_4634alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base G (37/44);  new base C (3/2);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

GAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTT  >  NC_002947/5257550‑5257844
                                                                                                                                                    |                                                                                                                                                   
gAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGgc                                                                                                                                                    >  3:57679/1‑149 (MQ=255)
     cACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTgg                                                                                                                                               <  3:201439/149‑1 (MQ=255)
     cACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTgg                                                                                                                                               >  2:117174/1‑149 (MQ=255)
     cACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTgg                                                                                                                                               <  2:299639/149‑1 (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:95323/1‑149 (MQ=255)
                   ggccggGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGc                                                                                                                                 <  7:79959/149‑1 (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                                                                                                                       >  5:163476/1‑148 (MQ=255)
                              gCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCcag                                                                                                                       >  7:317712/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                                                                                                        <  2:229679/148‑1 (MQ=255)
                                             gTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAgg                                                                                                        <  3:147336/148‑1 (MQ=255)
                                             gTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAgg                                                                                                        <  5:107101/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                                                                                           <  5:2267/148‑1 (MQ=255)
                                                          acaTCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCgg                                                                                           >  6:2267/1‑148 (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                                                                       <  1:346090/149‑1 (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)
                                                                              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:117174/149‑1 (MQ=255)
                                                                                       ggCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGtt                                                             <  1:269737/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                                          <  2:256394/149‑1 (MQ=255)
                                                                                                          cGGTAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGgc                                          >  1:256394/1‑149 (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                                                                         <  6:259578/100‑1 (MQ=255)
                                                                                                                             gCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATg                                                                         >  5:259578/1‑100 (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)
                                                                                                                                             cgaGCGGGCGGTGGGGATCAGGTTGGGGGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCt                                 <  8:244733/122‑1 (MQ=255)
                                                                                                                                             ccAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGGTACTCATGGCGCATGTCCt                                 >  7:244733/1‑124 (MQ=255)
                                                                                                                                              cAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCa          >  6:17718/1‑146 (MQ=255)
                                                                                                                                                gCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGAt     <  4:236077/149‑1 (MQ=255)
                                                                                                                                                gCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGAt     <  3:75186/149‑1 (MQ=255)
                                                                                                                                                 cGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGAtt    <  6:194521/149‑1 (MQ=255)
                                                                                                                                                   ggCGCTGGGGATCAGGTTGGGCGTTGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGAtttt  <  8:257510/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAAACCACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCG‑AAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTT  >  NC_002947/5257550‑5257844

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