Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,421,334 C→A 6.1% D174Y (GAC→TAC)  PP_4761 ← haloacid dehalogenase‑like family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,421,3340CA 6.1% 210.2 / 2.6 65D174Y (GAC→TAC) PP_4761haloacid dehalogenase‑like family hydrolase
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base A (2/2);  total (32/33)
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.68e-01

CCACAGGCCCCGAGACGTGCGCAATACCTGTGGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTGGCAATCTTCAGCAGCCCATCCGGGCTCGGCTT  >  NC_002947/5421186‑5421475
                                                                                                                                                    |                                                                                                                                              
ccacAGGCCCCGAGACGTGCGCAATACCTGTGGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTc                                                                                                                                                >  1:122634/1‑149 (MQ=255)
  acaGGCCCCGAGACGTGCGCAATACCTGTGGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCg                                                                                                                                              <  6:260482/149‑1 (MQ=255)
      gCCCCGAGACGTGCGCAATACCTGTGGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAg                                                                                                                                          <  8:97481/149‑1 (MQ=255)
       ccccGAGACGTGCGCAATACCTGTGGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAgg                                                                                                                                         <  1:49802/149‑1 (MQ=255)
             gaCGTGCGCAATACCTGTGGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCa                                                                                                                                           <  2:72755/141‑1 (MQ=255)
             gaCGTGCGCAATACCTGTGGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCa                                                                                                                                           >  1:72755/1‑141 (MQ=255)
              aCGTGCGCAATACCTGTGGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCa                                                                                                                                  >  5:218859/1‑149 (MQ=255)
                     cAATACCTGTGGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGt                                                                                                                           <  2:122634/149‑1 (MQ=255)
                      aaTACCTGTGGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTg                                                                                                                          >  2:33767/1‑149 (MQ=255)
                            tgtgGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGcc                                                                                                                    >  7:111898/1‑149 (MQ=255)
                             gtgGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGc                                                                                                                                            >  2:103542/1‑124 (MQ=255)
                             gtgGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGc                                                                                                                                            <  1:103542/124‑1 (MQ=255)
                             gtgGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGccc                                                                                                                   >  7:153084/1‑149 (MQ=255)
                             gtgGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGccc                                                                                                                   >  2:137528/1‑149 (MQ=255)
                                 aaGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTACGGCAGGTTCACCAGCACCGTGCGCGCCCccgc                                                                                                               >  2:130760/1‑149 (MQ=255)
                                 aaGCGGGCATGCCCGCGAAGAGGCCAGCATTGCCAACACACGCGTCAGCCCAGCATCACCATCAACGCCCGGCAATACGCCGCATCCCAATCCACCAGCTCTGGCCAGCGGTTGTACGGACGGTTCACCAGCACCCTGCGCGccccccc                                                                                                               >  8:151354/1‑147 (MQ=255)
                                   gCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCccgccg                                                                                                             >  8:237847/1‑149 (MQ=255)
                                     gggCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCgccgcc                                                                                                           >  8:53173/1‑149 (MQ=255)
                                       gCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGccgcccg                                                                                                         <  8:111898/149‑1 (MQ=255)
                                        cATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGccgcccgc                                                                                                        <  7:237847/149‑1 (MQ=255)
                                          tGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGccgcccgccc                                                                                                      <  7:59312/149‑1 (MQ=255)
                                            cccGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGccgcccgcccgc                                                                                                    <  6:278797/149‑1 (MQ=255)
                                            cccGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGccgcccgcccgc                                                                                                    >  3:247408/1‑149 (MQ=255)
                                                    gagGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAgg                                                                                            <  4:247408/149‑1 (MQ=255)
                                                      ggCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGt                                                                                           >  2:208396/1‑148 (MQ=255)
                                                       gCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTc                                                                                          >  1:106965/1‑148 (MQ=255)
                                                          cGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCa                                                                                                   <  3:233367/136‑1 (MQ=255)
                                                          cGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCa                                                                                                   >  4:233367/1‑136 (MQ=255)
                                                          cGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAg                                                                                      <  1:208396/149‑1 (MQ=255)
                                                          cGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAg                                                                                      <  2:169494/149‑1 (MQ=255)
                                                                 ccAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAAt                                                                               >  1:24071/1‑149 (MQ=255)
                                                                   aaCACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATcc                                                                             <  2:106965/149‑1 (MQ=255)
                                                                    acacacGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATccc                                                                            <  2:24071/149‑1 (MQ=255)
                                                                      acacGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATccccc                                                                          <  7:156256/149‑1 (MQ=255)
                                                                          gcgTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCccacca                                                                      <  8:196248/149‑1 (MQ=255)
                                                                            gTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATc                                                                    <  7:53173/149‑1 (MQ=255)
                                                                                            cTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCg                                                   >  8:82290/1‑149 (MQ=255)
                                                                                            cTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCg                                                   >  1:230600/1‑149 (MQ=255)
                                                                                             ttCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGt                                                  >  3:100204/1‑149 (MQ=255)
                                                                                                    cccGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTaa                                                                                <  5:273727/113‑1 (MQ=255)
                                                                                                    cccGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTaa                                                                                >  6:273727/1‑113 (MQ=255)
                                                                                                       ggCAATCCGCCGCATGCCAATCCACCCGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCCCCCGCAGGCCCGGTCAAAGCGGTAATCCCCCACCCTCCCCAGTTCGCTCCGCGTCACGcccccg                                        >  2:305337/1‑147 (MQ=255)
                                                                                                        gCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAg                                                                                      >  1:216665/1‑103 (MQ=255)
                                                                                                        gCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAg                                                                                      <  2:216665/103‑1 (MQ=255)
                                                                                                           aTCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTACGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGGCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCAGGCGCt                                   <  7:151354/149‑1 (MQ=255)
                                                                                                            tCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCt                                   >  4:136273/1‑149 (MQ=255)
                                                                                                            tCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTACGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCt                                   <  1:130760/149‑1 (MQ=255)
                                                                                                              cgccgcATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCAACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTg                                  <  1:33767/149‑1 (MQ=255)
                                                                                                              cgccgcATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTgg                                 <  1:137528/149‑1 (MQ=255)
                                                                                                              cgccgcATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTg                                  <  4:100204/148‑1 (MQ=255)
                                                                                                                     tGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTGGCAATCtt                          <  2:230600/149‑1 (MQ=255)
                                                                                                                     tGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTGGCAATCtt                          <  8:153084/149‑1 (MQ=255)
                                                                                                                      gCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTGGCAATCtt                          >  4:193716/1‑148 (MQ=255)
                                                                                                                       ccAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGccgccc                                                                                                          >  8:288574/1‑68 (MQ=255)
                                                                                                                       ccAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGccgccc                                                                                                          <  7:288574/68‑1 (MQ=255)
                                                                                                                        cAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCCCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCCGGCGCGGGCAAACTTcag                       >  6:100898/1‑149 (MQ=255)
                                                                                                                              agtaGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTgg                                 <  4:46870/130‑1 (MQ=255)
                                                                                                                              accaGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCACTCCAGGTCAAAGCGGTAATCCCCCACCAACAACAGTTCCCTCGGCGTCACGCCCCCGGCGCTGGCAATCTTCAGAAGcc                  >  7:127879/1‑148 (MQ=255)
                                                                                                                               ccaGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCCCCGCCCGCCCGCAGTGCAGGTCCAAGCGGTAATCCCCCCCCATCACCAGTTCGCTCGGCGTCACGCGCCAGGCGCGGGGAATCTTCAGCCGCCCa                >  6:314361/1‑149 (MQ=255)
                                                                                                                                caGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGt                                                  >  1:37817/1‑114 (MQ=255)
                                                                                                                                caGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGt                                                  <  2:37817/114‑1 (MQ=255)
                                                                                                                                  gCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTGGCAATCTTCAGCAGCCCATcc             >  3:190179/1‑149 (MQ=255)
                                                                                                                                          cAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTGGCAATCTTCAGCAGCCCATCCGGGCTCgg     >  4:248877/1‑149 (MQ=255)
                                                                                                                                           aGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTGGCAATCtt                          <  8:195157/127‑1 (MQ=255)
                                                                                                                                           aGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTGGCAATCtt                          >  7:195157/1‑127 (MQ=255)
                                                                                                                                             gggTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTGGCAATCTTCAGCAGCCCATCCGGGCTCGGCtt  <  3:193716/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CCACAGGCCCCGAGACGTGCGCAATACCTGTGGAAGCGGGCATGCCCGCGAAGAGGCCCGCATTGCCAACACACGCGTCAGCCCAGCATCACCTTCAACGCCCGGCAATCCGCCGCATGCCAATCCACCAGCTCTGGCCAGGGGTTGTCCGGCAGGTTCACCAGCACCGTGCGCGCCCCCGCCGCCCGCCCGCAGTCCAGGTCAAAGCGGTAATCCCCCACCATCACCAGTTCGCTCGGCGTCACGCCCCAGGCGCTGGCAATCTTCAGCAGCCCATCCGGGCTCGGCTT  >  NC_002947/5421186‑5421475

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