Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,257,704 G→T 5.2% P33H (CCC→CAC)  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,7040GT 5.2% 298.1 / 3.2 93P33H (CCC→CAC) PP_4634alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base G (46/42);  new base T (3/2);  total (49/44)
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.65e-01

ACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCT  >  NC_002947/5257556‑5257851
                                                                                                                                                    |                                                                                                                                                   
aCGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTggg                                                                                                                                                     >  8:312791/1‑149 (MQ=255)
aCGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTggg                                                                                                                                                     >  3:119995/1‑149 (MQ=255)
aCGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTggg                                                                                                                                                     <  4:119995/149‑1 (MQ=255)
   cAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGAt                                                                                                                                                  <  6:332264/149‑1 (MQ=255)
     gCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCCCCAGCGCCCGCTGGTGATCa                                                                                                                                                <  7:8803/149‑1 (MQ=255)
     gCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCa                                                                                                                                                <  8:325806/149‑1 (MQ=255)
       gAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAgg                                                                                                                                              <  1:162106/149‑1 (MQ=255)
           ccggccggGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTgg                                                                                                                                          >  5:305178/1‑149 (MQ=255)
                     ttCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGcc                                                                                                                                <  4:401088/149‑1 (MQ=255)
                       cGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCca                                                                                                                               <  4:195006/148‑1 (MQ=255)
                        gCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTgg                                                                                                                                          >  6:281071/1‑136 (MQ=255)
                        gCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTgg                                                                                                                                          <  5:281071/136‑1 (MQ=255)
                        gCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCcag                                                                                                                              <  3:102911/148‑1 (MQ=255)
                        gCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCcag                                                                                                                              >  5:174989/1‑148 (MQ=255)
                        gCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCcag                                                                                                                              >  1:7648/1‑148 (MQ=255)
                        gCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCcag                                                                                                                              >  1:393200/1‑148 (MQ=255)
                        gCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCcag                                                                                                                              <  8:281903/148‑1 (MQ=255)
                            aTGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGaa                                                                                                                         <  8:86053/149‑1 (MQ=255)
                                cACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAgccg                                                                                                                     >  5:115892/1‑149 (MQ=255)
                                cACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAgccg                                                                                                                     >  7:379541/1‑149 (MQ=255)
                                    gTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGa                                                                                                                 <  2:7648/149‑1 (MQ=255)
                                     tAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCATAGCGGGCGTAATAGCGGCAGCGTACCGAACTGGCAGCGTTGATGATCACCAGCGGGCTCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAg                                                                                                                <  3:60459/149‑1 (MQ=255)
                                       gTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAgg                                                                                                               <  6:70885/148‑1 (MQ=255)
                                        tCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGa                                                                                                                 >  4:398247/1‑145 (MQ=255)
                                        tCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGa                                                                                                                 <  3:398247/145‑1 (MQ=255)
                                        tCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCt                                                                                                             >  2:208502/1‑149 (MQ=255)
                                        tCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCt                                                                                                             >  5:330284/1‑149 (MQ=255)
                                          gTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTg                                                                                                            >  4:394421/1‑148 (MQ=255)
                                              gTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATcc                                                                                                       >  3:30354/1‑149 (MQ=255)
                                              gTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATcc                                                                                                       >  6:260563/1‑149 (MQ=255)
                                                  gCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCg                                                                                                   >  1:24120/1‑149 (MQ=255)
                                                  gCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCg                                                                                                   >  3:151633/1‑149 (MQ=255)
                                                      aTCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTc                                                                                                    <  3:210410/144‑1 (MQ=255)
                                                      aTCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTc                                                                                                    >  4:210410/1‑144 (MQ=255)
                                                         gCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCgg                                                                                            <  2:24120/149‑1 (MQ=255)
                                                          cAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCggg                                                                                           >  5:329404/1‑149 (MQ=255)
                                                          cAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCggg                                                                                           <  5:260563/149‑1 (MQ=255)
                                                              cGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTgg                                                                                       <  3:176131/149‑1 (MQ=255)
                                                               gTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGc                                                                                      >  1:225625/1‑149 (MQ=255)
                                                                tGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTTATGATCACGAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCg                                                                                     <  6:305178/149‑1 (MQ=255)
                                                                   gCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCgg                                                                                            >  7:363426/1‑139 (MQ=255)
                                                                   gCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCgg                                                                                            <  8:363426/139‑1 (MQ=255)
                                                                    cGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAAt                                                                                 <  6:115892/149‑1 (MQ=255)
                                                                     gAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATg                                                                                >  7:192731/1‑149 (MQ=255)
                                                                       aCAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATggg                                                                              <  7:312791/149‑1 (MQ=255)
                                                                         aGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGAGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATgggg                                                                             >  4:172797/1‑148 (MQ=255)
                                                                         aGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGCCAGCGTTGATGATAACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGACGCCGAGGATGTATCCGTCGGCCGCGGGGTGGCGAAATgggg                                                                             >  5:273400/1‑148 (MQ=255)
                                                                           gTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAg                                                                          <  6:329404/149‑1 (MQ=255)
                                                                           gTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAg                                                                          >  5:207937/1‑149 (MQ=255)
                                                                           gTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGAACACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAg                                                                          >  6:111650/1‑149 (MQ=255)
                                                                            ttttCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGTCGCCAGCATAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGc                                                                         <  6:273400/146‑1 (MQ=255)
                                                                            tAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGc                                                                         <  8:379541/149‑1 (MQ=255)
                                                                            tAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGAGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGc                                                                         >  7:361124/1‑149 (MQ=255)
                                                                                cGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGt                                                                     >  8:359034/1‑149 (MQ=255)
                                                                                 ggCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGtt                                                                    >  2:112484/1‑149 (MQ=255)
                                                                                  gCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTg                                                                   <  7:359034/149‑1 (MQ=255)
                                                                                    aaaGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCg                                                                                     <  2:245280/129‑1 (MQ=255)
                                                                                    aaaGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCg                                                                                     >  1:245280/1‑129 (MQ=255)
                                                                                     atgCGGGCGTAATTGCGGCAGCGTACCGAAGTTGCAGCGTTGATTATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGAt                                                                <  6:99865/147‑1 (MQ=255)
                                                                                       gCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGAt                                                                >  3:380338/1‑147 (MQ=255)
                                                                                          ggCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGc                                                           >  6:349704/1‑149 (MQ=255)
                                                                                           gCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGc                                                           >  4:75965/1‑148 (MQ=255)
                                                                                                aTAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGTGTTGGATGAGGCTACTc                                                      <  6:207937/148‑1 (MQ=255)
                                                                                                 tAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCAt                                                    >  4:47368/1‑149 (MQ=255)
                                                                                                    cGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGgc                                                 <  2:144245/149‑1 (MQ=255)
                                                                                                    cGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGgc                                                 >  7:209837/1‑149 (MQ=255)
                                                                                                      gCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGgcg                                                >  4:310462/1‑148 (MQ=255)
                                                                                                       cAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCa                                              <  8:361124/149‑1 (MQ=255)
                                                                                                               cGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGACAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCtgt                                      >  2:351245/1‑149 (MQ=255)
                                                                                                                gAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCtgtt                                     <  4:151633/149‑1 (MQ=255)
                                                                                                                gAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCtgtt                                     <  8:209837/149‑1 (MQ=255)
                                                                                                                 aaGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTgg                                                                                       >  1:169949/1‑98 (MQ=255)
                                                                                                                 aaGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTgg                                                                                       <  2:169949/98‑1 (MQ=255)
                                                                                                                 aaGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCtgttg                                    <  3:75965/149‑1 (MQ=255)
                                                                                                                 aaGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCtgttg                                    <  1:112484/149‑1 (MQ=255)
                                                                                                                   gTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCtgttgtt                                  <  8:192731/149‑1 (MQ=255)
                                                                                                                   gTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCtgttgtt                                  <  6:330284/149‑1 (MQ=255)
                                                                                                                            ttGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCa                          >  7:205242/1‑148 (MQ=255)
                                                                                                                              gatgatCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCaaaa                       >  4:279431/1‑149 (MQ=255)
                                                                                                                               ctgatCCCCAGCGCCCGCTGGTGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTc                                          >  8:159855/2‑129 (MQ=255)
                                                                                                                               ctgatCCCCAGCGCCCGCTGGTGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTc                                          <  7:159855/128‑1 (MQ=255)
                                                                                                                               atgatCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCaaaa                       <  4:30354/148‑1 (MQ=255)
                                                                                                                                tgatCCCCAGCGCCCGCTGGTGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAAc                      >  7:69232/1‑148 (MQ=255)
                                                                                                                                tgatCCCCAGCGCCCGCTGGTGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAAc                      >  5:405596/1‑148 (MQ=255)
                                                                                                                                tgatCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAAc                      <  4:380338/148‑1 (MQ=255)
                                                                                                                                 gatCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAATGGGGTAGAGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAg                   >  1:302472/1‑149 (MQ=255)
                                                                                                                                  atCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAg                   >  3:356112/1‑149 (MQ=255)
                                                                                                                                  atCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAg                   >  6:302061/1‑149 (MQ=255)
                                                                                                                                   tCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAg                   >  7:312552/1‑148 (MQ=255)
                                                                                                                                    cACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTa                    >  1:318581/1‑146 (MQ=255)
                                                                                                                                    cACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTa                    <  2:318581/146‑1 (MQ=255)
                                                                                                                                     aCCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCa                 >  7:327906/1‑148 (MQ=255)
                                                                                                                                       cAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTa                    >  7:241403/1‑143 (MQ=255)
                                                                                                                                       cAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTa                    <  8:241403/143‑1 (MQ=255)
                                                                                                                                         gCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGAt            <  5:349704/149‑1 (MQ=255)
                                                                                                                                               gcTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTg        >  2:278464/1‑147 (MQ=255)
                                                                                                                                                   gggATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCt  <  6:174989/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGCAGCGAAGCCGGCCGGGATTCGCCGATGCCACGGTAGTCGTAGGTCAGCACATCGCAGCCGTGGGCGAACAGGTAGTCGGCAAAGCGGGCGTAATAGCGGCAGCGTACCGAAGTGGCAGCGTTGATGATCACCAGCGGGCGCTGGGGATCAGGTTGGGCGTGGCGCCAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCT  >  NC_002947/5257556‑5257851

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: