Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,908,460 A→G 6.8% N122S (AAC→AGC)  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,4600AG 6.8% 256.3 / 6.7 88N122S (AAC→AGC) PP_3450TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (40/42);  new base G (4/2);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

ACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAG  >  NC_002947/3908321‑3908607
                                                                                                                                           |                                                                                                                                                   
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                                                                                                                                               <  1:784210/113‑1 (MQ=255)
                                 cTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCgcgc                                                                                                                                               >  2:784210/1‑113 (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)
                                                 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                                                                                                     >  1:229466/1‑123 (MQ=255)
                                                                 aaaGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCa                                                                                                     <  2:229466/123‑1 (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                                                  <  2:1142342/150‑1 (MQ=255)
                                                                                         ccACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTg                                                  >  2:1960803/1‑150 (MQ=255)
                                                                                         ccACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTg                                                  >  1:1924283/1‑150 (MQ=255)
                                                                                          cACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTgg                                                 <  2:31081/150‑1 (MQ=255)
                                                                                           aCGCCCCCTCCCCCGCAGCGTCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGc                                                >  1:292189/1‑150 (MQ=255)
                                                                                           aCGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagc                                                                  >  2:470473/1‑132 (MQ=255)
                                                                                           aCGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagc                                                                  <  1:470473/132‑1 (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                                        >  1:2384772/1‑150 (MQ=255)
                                                                                                   tCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGcc                                        >  2:2370654/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                                                                 <  1:1030388/102‑1 (MQ=255)
                                                                                                                          gcgGTTGGCCATGGCCAGCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagca                                                                 >  2:1030388/1‑102 (MQ=255)
                                                                                                                           cgGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTg                <  2:782269/150‑1 (MQ=255)
                                                                                                                             gCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGt              <  1:1985518/150‑1 (MQ=255)
                                                                                                                              cTGGCCATGGCCAACCGGGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTg             >  2:1235410/1‑150 (MQ=255)
                                                                                                                               tggccatggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTgc            <  2:1497026/150‑1 (MQ=255)
                                                                                                                                ggccatggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTgcg           <  1:1155153/150‑1 (MQ=255)
                                                                                                                                 gccatggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTgcgg          >  1:1247802/1‑150 (MQ=255)
                                                                                                                                  ccatggccaGCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAACGCCAGCAACTTTGGCATCGTGCCCGCGAAGTCGACGATGCAGCGCTGGTccgga         >  2:1600148/1‑149 (MQ=255)
                                                                                                                                   catggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACttt                                                   >  2:1385378/1‑107 (MQ=255)
                                                                                                                                   catggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACttt                                                   <  1:1385378/107‑1 (MQ=255)
                                                                                                                                   catggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTgcggc         <  1:821001/149‑1 (MQ=255)
                                                                                                                                   catggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTgcggc         >  2:821001/1‑149 (MQ=255)
                                                                                                                                    atggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGgcgc       <  1:1235410/150‑1 (MQ=255)
                                                                                                                                        ccaACCGCGCCATTGCCGTTGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGa   <  1:406907/150‑1 (MQ=255)
                                                                                                                                        ccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGa   <  2:1924597/150‑1 (MQ=255)
                                                                                                                                        ccaACCGCGCCATTGCCGGCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCTCCAGCAACTTTGGCAACTGTCCCGCGAAGACGACAATGCAGCGCTGGTGCGGCGCCTGa   >  2:496993/1‑150 (MQ=255)
                                                                                                                                         caACCGCGCCATTGCCGTCGGCCACCCCCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAg  <  1:249877/150‑1 (MQ=255)
                                                                                                                                         caACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGACCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAg  >  1:1465974/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
ACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAG  >  NC_002947/3908321‑3908607

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