Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,041,016 A→C 6.3% Q192H (CAA→CAC)  PP_2652 → hydratase/decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,041,0160AC 6.3% 260.7 / 6.2 94Q192H (CAA→CAC) PP_2652hydratase/decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (45/42);  new base C (3/3);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TCGACCGCAAGGCTGTCGGCTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGC  >  NC_002947/3040868‑3041164
                                                                                                                                                    |                                                                                                                                                    
tCGACCGCAAGGCTGTCGGCTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAAc                                                                                                                                                     <  2:2049217/150‑1 (MQ=255)
  gACCGCAAGGCTGTCGGCTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCACCgg                                                                                                                                                   >  1:557855/1‑150 (MQ=255)
  gACCGCAAGGCTGTCGGCTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGAGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAAc                                                                                                                                                     <  2:2045980/148‑1 (MQ=255)
  gACCGCAAGGCTGTCGGCTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGAGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAAc                                                                                                                                                     >  1:2045980/1‑148 (MQ=255)
       cAAGGCTGTCGGCTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGgcgcgc                                                                                                                                              >  1:389558/1‑150 (MQ=255)
         aGGCTGTCGGCTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGACGCGCCg                                                                                                                                            >  1:2380881/1‑150 (MQ=255)
           gCTGTCGGCTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGt                                                                                                                                          <  1:2358229/150‑1 (MQ=255)
            cTGTCGGCTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTg                                                                                                                                         >  2:2014037/1‑150 (MQ=255)
            cTGTCGGCTGGCCAAGTTTTGTGGCTGAAAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCACCGGCGCGCCGTTg                                                                                                                                         <  2:557855/150‑1 (MQ=255)
             tGTCGGCTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCACCGGCGCGCCGTTgg                                                                                                                                        <  1:1770813/150‑1 (MQ=255)
             tGTCGGCTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTCAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTgg                                                                                                                                        <  2:1486026/150‑1 (MQ=255)
                   cTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTc                                                                                                                                  <  2:389558/150‑1 (MQ=255)
                      gCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACAGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTgg                                                                                                                               >  2:857871/1‑150 (MQ=255)
                       ccAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGc                                                                                                                              >  1:2098932/1‑150 (MQ=255)
                       ccAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGc                                                                                                                              <  2:2098932/150‑1 (MQ=255)
                       ccAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGc                                                                                                                              >  2:142693/1‑150 (MQ=255)
                        cAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCg                                                                                                                             >  1:277214/1‑150 (MQ=255)
                         aaGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGa                                                                                                                            >  1:320658/1‑150 (MQ=255)
                         aaGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGATTTTTGCTCGCCGGTGGCTTGTCTGGCGa                                                                                                                            >  1:602420/1‑150 (MQ=255)
                          aGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGAt                                                                                                                           <  1:1384702/150‑1 (MQ=255)
                            ttttGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGc                                                                                                                         >  2:1260065/1‑150 (MQ=255)
                                 tgGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAAGAGCCTGGGCGCGGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCa                                                                                                                    <  2:631445/150‑1 (MQ=255)
                                  ggCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCAc                                                                                                                   >  2:1647896/1‑150 (MQ=255)
                                   gCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCCGCcccc                                                                                                                  >  2:449227/1‑150 (MQ=255)
                                          aCGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGc                                                                                                           <  1:1260065/150‑1 (MQ=255)
                                            gTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTg                                                                                                         <  2:1259638/150‑1 (MQ=255)
                                                ggTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGACACCGCGGATGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGgtgt                                                                                                     >  2:1775953/1‑150 (MQ=255)
                                                 gTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTc                                                                                                    >  1:1720234/1‑150 (MQ=255)
                                                  tttGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCg                                                                                                   >  1:1347819/1‑150 (MQ=255)
                                                   ttGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGc                                                                                                  <  2:320658/150‑1 (MQ=255)
                                                    tGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCt                                                                                                 <  1:327200/150‑1 (MQ=255)
                                                      aaGCATTCGTGTTCGGTAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTgg                                                                                               <  2:1995098/150‑1 (MQ=255)
                                                         cATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGGGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTc                                                                                            >  2:2192166/1‑150 (MQ=255)
                                                          aTTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCa                                                                                           >  2:1968167/1‑150 (MQ=255)
                                                           ttGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCa                                                                                           <  2:1245360/149‑1 (MQ=255)
                                                           ttGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCa                                                                                           >  1:1245360/1‑149 (MQ=255)
                                                               tggtCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTgg                                                                                               <  2:1026892/141‑1 (MQ=255)
                                                               tggtCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTgg                                                                                               >  1:1026892/1‑141 (MQ=255)
                                                               tggtCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGt                                                                                      >  2:1065576/1‑150 (MQ=255)
                                                                ggtCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTa                                                                                     >  1:1938984/1‑150 (MQ=255)
                                                                ggtCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTa                                                                                     <  1:1647896/150‑1 (MQ=255)
                                                                 gtCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTAc                                                                                    >  1:1648285/1‑150 (MQ=255)
                                                                  tCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCa                                                                                                                    <  2:1560456/117‑1 (MQ=255)
                                                                  tCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCa                                                                                                                    >  1:1560456/1‑117 (MQ=255)
                                                                    ggCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGca                                                                                 >  1:1532985/1‑150 (MQ=255)
                                                                       aaTAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCacac                                                                              <  1:1826574/150‑1 (MQ=255)
                                                                            ccTGGGCGATGGTATAAATGTGCCGCTGCTCGCCGACCCGGCCGACACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCATGTCTGGCGATGCGGCCACCGACCCGCTGGTGTTGCTGGATCACATGTGGGCACACGCTGc                                                                         <  2:532500/150‑1 (MQ=255)
                                                                                    cTGGTCTAAATGTGTTTTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTg                                                                 <  2:661400/150‑1 (MQ=255)
                                                                                        tataAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCa                                                             <  2:2188669/150‑1 (MQ=255)
                                                                                          taAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGa                                                           >  1:285395/1‑150 (MQ=255)
                                                                                                 gCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGaa                                                    <  1:1065576/150‑1 (MQ=255)
                                                                                                    gTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCt                                                                                        >  2:2311802/1‑111 (MQ=255)
                                                                                                    gTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCt                                                                                        <  1:2311802/111‑1 (MQ=255)
                                                                                                    gTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGc                                                 >  2:1595668/1‑150 (MQ=255)
                                                                                                     ttGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTg                                                      >  2:823453/1‑144 (MQ=255)
                                                                                                     ttGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTg                                                      <  1:823453/144‑1 (MQ=255)
                                                                                                      tGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCg                                               <  2:964522/150‑1 (MQ=255)
                                                                                                       gCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCgg                                              <  2:952569/150‑1 (MQ=255)
                                                                                                         tCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCgggg                                            >  1:185738/1‑150 (MQ=255)
                                                                                                         tCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCgggg                                            >  2:2280619/1‑150 (MQ=255)
                                                                                                                ccTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGt                                     <  1:1595668/150‑1 (MQ=255)
                                                                                                                ccTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGACGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGt                                     <  2:2380881/150‑1 (MQ=255)
                                                                                                                 cTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTg                                    >  1:1939300/1‑150 (MQ=255)
                                                                                                                 cTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTg                                    <  1:1281559/150‑1 (MQ=255)
                                                                                                                  tGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGa                                   <  2:602420/150‑1 (MQ=255)
                                                                                                                  tGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGa                                   <  2:2129422/150‑1 (MQ=255)
                                                                                                                  tGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGa                                   <  1:1914201/150‑1 (MQ=255)
                                                                                                                   ggCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGac                                  <  1:142693/150‑1 (MQ=255)
                                                                                                                    gCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGacc                                 >  2:2005673/1‑150 (MQ=255)
                                                                                                                      cGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGaccac                               >  1:2303300/1‑150 (MQ=255)
                                                                                                                      cGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTATGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGaccac                               >  1:1352625/1‑150 (MQ=255)
                                                                                                                        aaaCCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCg                             >  1:1720998/1‑150 (MQ=255)
                                                                                                                          accaccGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCggg                           <  1:2014037/150‑1 (MQ=255)
                                                                                                                            caccGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGc                         >  2:342663/1‑150 (MQ=255)
                                                                                                                             accGAGGTGCATGTGGATGGCCACCGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGc                                                                            <  2:2343662/98‑1 (MQ=255)
                                                                                                                             accGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCAc                                                                                          >  1:2359187/1‑84 (MQ=255)
                                                                                                                             accGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCAc                                                                                          <  2:2359187/84‑1 (MQ=255)
                                                                                                                             accGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGc                                                                            >  1:2343662/1‑98 (MQ=255)
                                                                                                                              ccGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGa                                           >  2:1625693/1‑130 (MQ=255)
                                                                                                                              ccGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGa                                           <  1:1625693/130‑1 (MQ=255)
                                                                                                                               cGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCAt                      <  1:440368/150‑1 (MQ=255)
                                                                                                                                 aGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCAtgt                    <  2:1532985/150‑1 (MQ=255)
                                                                                                                                      cATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAg               <  1:521220/150‑1 (MQ=255)
                                                                                                                                         gtgGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGc                                                                         <  2:1128054/89‑1 (MQ=255)
                                                                                                                                         gtgGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGc                                                                         >  1:1128054/1‑89 (MQ=255)
                                                                                                                                           ggATGGCCAACGGCGTGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACTCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGtt          <  2:221600/150‑1 (MQ=255)
                                                                                                                                            gATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGCTGCGGCCACCGACCCGCTGGTGTCGCTGCCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCg                                               <  2:2171115/112‑1 (MQ=255)
                                                                                                                                            gATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCg                                               >  1:2171115/1‑112 (MQ=255)
                                                                                                                                            gATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTc         >  1:778601/1‑150 (MQ=255)
                                                                                                                                            gATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTc         >  2:1425416/1‑150 (MQ=255)
                                                                                                                                              aagCCACCGGCGCGCCGTTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGa       >  1:912301/3‑150 (MQ=255)
                                                                                                                                              aagCCACCGGCGCGCCGTTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGa       >  2:1121082/3‑150 (MQ=255)
                                                                                                                                                 ccAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGAttt    >  1:346661/1‑150 (MQ=255)
                                                                                                                                                 ccAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTCCCGAGCGTGGCCAGACCTTGAAGGCCGGGGCTGTCGTGACCACCGGGGCCCTGTGCCAGCCGTTCGAttt    >  2:74970/1‑150 (MQ=255)
                                                                                                                                                   aaCGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGc  <  2:277214/150‑1 (MQ=255)
                                                                                                                                                   aaCGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGc  <  2:185738/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGACCGCAAGGCTGTCGGCTGGCCAAGTTTTGTGGCTGACAACGTCGGGTTTGAAGCATTGGTGGTCGGCAATAGCCTGGGCGCTGGTATAAATGTGCCGTTGCTCGCCGACCTGGCCGAAACCACCGAGGTGCATGTGGATGGCCAACGGCGCGCCGGTGGCTTGTCTGGCGATGCGGCCACCGACCCGCTGGTGTCGCTGGCTCACCTGTACGCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGC  >  NC_002947/3040868‑3041164

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑