Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,908,447 C→T 6.2% L118L (CTG→TTG)  PP_3450 → TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,908,4470CT 6.2% 241.0 / 4.9 80L118L (CTG→TTG) PP_3450TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (38/37);  new base T (3/2);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

TGCAGCCGCTGGCCGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTG  >  NC_002947/3908304‑3908596
                                                                                                                                               |                                                                                                                                                     
ttcaaccGCTGGCCGACACAGGCCATTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCtggcca                                                                                                                                                 >  2:2251533/6‑150 (MQ=255)
 gCAGCCGCTGGCCGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCtggccat                                                                                                                                                <  2:2375783/150‑1 (MQ=255)
      cGCTGGCCGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCtggccatggcca                                                                                                                                           >  1:1737292/1‑150 (MQ=255)
             cGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCtggccatggc                                                                                                                                             >  1:1288367/1‑141 (MQ=255)
             cGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCtggccatggc                                                                                                                                             <  2:1288367/141‑1 (MQ=255)
                 acaGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCAtt                                                                                                                                <  2:1861539/150‑1 (MQ=255)
                   aGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGc                                                                                                                              >  1:110516/1‑150 (MQ=255)
                    ggCCAGTTGCCTGCGCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGCCGCGGCTGGCCATGGCCAACCGCGCCATTGcc                                                                                                                             <  2:11100/150‑1 (MQ=255)
                     gCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCg                                                                                                                            >  1:1316063/1‑150 (MQ=255)
                     gCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGAATGTGCTGCTTGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCg                                                                                                                            >  1:1142342/1‑150 (MQ=255)
                              ctgccCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccacc                                                                                                                   >  1:2077053/1‑150 (MQ=255)
                               tgccCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccacca                                                                                                                  >  2:1509465/1‑150 (MQ=255)
                                 cccAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGTAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccaccacc                                                                                                                >  1:2008389/1‑150 (MQ=255)
                                  ccAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccaccacca                                                                                                               <  2:110516/150‑1 (MQ=255)
                                  ccAGGACTGGCAACGCCCGGCCAATGGCTATATCCAGCAATTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccaccacca                                                                                                               <  1:2221566/150‑1 (MQ=255)
                                      gACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGACAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCACAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACcagcag                                                                                                           <  2:1737292/150‑1 (MQ=255)
                                         tGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGcc                                                                                                        >  2:2267076/1‑150 (MQ=255)
                                           gCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCa                                                                                                      >  1:628806/1‑150 (MQ=255)
                                             aaCGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccaccacca                                                                                                               >  2:478285/1‑139 (MQ=255)
                                             aaCGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccaccacca                                                                                                               <  1:478285/139‑1 (MQ=255)
                                             aaCGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAAc                                                                                                    >  2:1675202/1‑150 (MQ=255)
                                                gCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgct                                                                                                 >  2:640847/1‑150 (MQ=255)
                                                gCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgct                                                                                                 <  1:2251533/150‑1 (MQ=255)
                                                gCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAAGGCCTGACTGCGCTGCTGGCCCGCGCCCCCTCCCCGGCCGCGCCTAACAGCTCTCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgct                                                                                                 >  1:782269/1‑150 (MQ=255)
                                                 ccTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgcgg                                                                                                <  2:628806/150‑3 (MQ=255)
                                                  cTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCgcgc                                                                                                                                    >  2:784210/1‑113 (MQ=255)
                                                  cTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCgcgc                                                                                                                                    <  1:784210/113‑1 (MQ=255)
                                                  cTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctgg                                                                                               <  1:39443/150‑1 (MQ=255)
                                                  cTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctgg                                                                                               <  2:1316063/150‑1 (MQ=255)
                                                  cTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGACCACGCCCCCTCCAAACAATCGCGCAGCCGGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctgg                                                                                               >  1:1497026/1‑150 (MQ=255)
                                                     gCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctggctg                                                                                            <  2:261508/150‑1 (MQ=255)
                                                      ccAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGc                                                                                           >  1:178010/1‑150 (MQ=255)
                                                        aaTGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAg                                                                                         <  1:1778028/150‑1 (MQ=255)
                                                          tGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGGCTACCTCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTc                                                                                       <  2:2392881/150‑1 (MQ=255)
                                                             cTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCtggccatgg                                                                                                                                              <  2:2041078/92‑1 (MQ=255)
                                                             cTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCtggccatgg                                                                                                                                              >  1:2041078/1‑92 (MQ=255)
                                                                  tCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGc                                                                               <  1:2138596/150‑1 (MQ=255)
                                                                   ccAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCGCGCCGCGGCACACCAGACGAGGCTAGCCGTGGCCAACCGCGCCAGCGCCGTCGGCCACCTCCAGCAGGCCCTACGCTGGCTGCCTTTGCTACCGcc                                                                              >  2:175825/1‑150 (MQ=255)
                                                                   ccAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGcc                                                                              <  1:640847/150‑1 (MQ=255)
                                                                       cAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGAt                                                                          <  1:1675202/150‑1 (MQ=255)
                                                                          gTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGTTGGCCATGGCCAGCcgc                                                                                                                                      >  1:694728/1‑87 (MQ=255)
                                                                          gTCGACGAAAAGGGCCGGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGTTGGCCATGGCCAGCcgc                                                                                                                                      <  2:694728/87‑1 (MQ=39)
                                                                             gacgaAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTg                                                                    >  2:506149/1‑150 (MQ=255)
                                                                               cgaAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTgca                                                                  <  2:2077053/150‑1 (MQ=255)
                                                                                gaAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTgcag                                                                 >  2:1155153/1‑150 (MQ=255)
                                                                                  aaaGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCa                                                                                          <  2:229466/123‑1 (MQ=255)
                                                                                  aaaGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCa                                                                                          >  1:229466/1‑123 (MQ=255)
                                                                                       gCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCg                                                          >  1:614865/1‑150 (MQ=255)
                                                                                           gACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagca                                                      <  1:927333/150‑1 (MQ=255)
                                                                                             cTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagcagc                                                    <  1:954490/150‑1 (MQ=255)
                                                                                              tGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCa                                                                                          <  2:988675/111‑1 (MQ=255)
                                                                                              tGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCa                                                                                          >  1:988675/1‑111 (MQ=255)
                                                                                              tGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGACATTGCCGTCGGCCACCACAAGCGGAATGCAGGCATGCGGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAgcg                                                   >  2:63569/1‑150 (MQ=255)
                                                                                                  ctgctgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAg                                               <  1:1877905/150‑1 (MQ=255)
                                                                                                     ctgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCACGCCCAACGCTGGCTGCAGCCGCCACCGCCTGACCAACCGCAGCAACCCGAGCAGCGCCCGCaa                                            >  2:249877/1‑150 (MQ=255)
                                                                                                      tgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAAc                                           >  2:1985518/1‑150 (MQ=255)
                                                                                                        gCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGGTCAACTGCAGCAACCCGAGCAGCGCCAGCAACtt                                         <  2:1425772/150‑1 (MQ=255)
                                                                                                          ccACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTg                                       <  1:1960803/150‑1 (MQ=255)
                                                                                                          ccACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTg                                       >  1:1924283/1‑150 (MQ=255)
                                                                                                          ccACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTg                                       <  2:1142342/150‑1 (MQ=255)
                                                                                                          ccACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTg                                       >  2:1960803/1‑150 (MQ=255)
                                                                                                           cACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTgg                                      <  2:31081/150‑1 (MQ=255)
                                                                                                            aCGCCCCCTCCCCCGCAGCGTCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGc                                     >  1:292189/1‑150 (MQ=255)
                                                                                                            aCGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagc                                                       <  1:470473/132‑1 (MQ=255)
                                                                                                            aCGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagc                                                       >  2:470473/1‑132 (MQ=255)
                                                                                                            aCGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGc                                     >  1:1924597/1‑150 (MQ=255)
                                                                                                            aCGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGc                                     >  2:2302323/1‑150 (MQ=255)
                                                                                                                ccccTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctgg                                                                                               >  1:2333813/1‑88 (MQ=255)
                                                                                                                ccccTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctgg                                                                                               <  2:2333813/88‑1 (MQ=255)
                                                                                                                    tCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGcc                             >  2:2370654/1‑150 (MQ=255)
                                                                                                                    tCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGcc                             >  1:2384772/1‑150 (MQ=255)
                                                                                                                        ccGCAGCGGCCAACCAGCCGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCcgcg                         <  1:819041/150‑1 (MQ=255)
                                                                                                                             gCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGc                    <  2:1924283/150‑1 (MQ=255)
                                                                                                                                 ccAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGcgacg                <  1:1509465/150‑1 (MQ=255)
                                                                                                                                  cAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCTGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGcgacga               >  2:1748237/1‑150 (MQ=255)
                                                                                                                                     aatggcgcgGTTGGCCATGGCCAGCCGCGTCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGc            >  2:1282418/6‑150 (MQ=255)
                                                                                                                                           gcgGTTGGCCATGGCCAGCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagca                                                      >  2:1030388/1‑102 (MQ=255)
                                                                                                                                           gcgGTTGGCCATGGCCAGCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagca                                                      <  1:1030388/102‑1 (MQ=255)
                                                                                                                                            cgGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTg     <  2:782269/150‑1 (MQ=255)
                                                                                                                                              gCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGt   <  1:1985518/150‑1 (MQ=255)
                                                                                                                                               cTGGCCATGGCCAACCGGGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTg  >  2:1235410/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
TGCAGCCGCTGGCCGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTG  >  NC_002947/3908304‑3908596

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: