Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,897,594 G→C intergenic (‑5372/‑3225) RPPX_08470 ← / → RPPX_08480 tryptophan synthase subunit beta/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099741,897,5940GC83.0% 5317.9 / inf 3386intergenic (‑5372/‑3225)RPPX_08470/RPPX_08480tryptophan synthase subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (249/248);  new base C (1223/1199);  total (1716/1670)
Fisher's exact test for biased strand distribution p-value = 8.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

Only 100 of 3401 total aligned reads displayed.

GCAACCGCCCCGGACAACGTCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGG  >  CP009974/1897445‑1897739
                                                                                                                                                     |                                                                                                                                                 
gtcactgccccggacaacgtctacgtcggtagcaacgacccggttgtgatgtcgatcgcaaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGc                                                                                                                                                  >  1:5067163‑M2/129‑151 (MQ=255)
     tgCCCCGGACAACGTCTATGTCGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCCGTTAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGc                                                                                                                                             >  2:1283469/2‑151 (MQ=37)
      gcaccagACAACGTCTATGTCGGCACCAATGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCa                                                                                                                                            >  2:3564923/7‑151 (MQ=37)
      gCCCCGGACAACGTCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCa                                                                                                                                            >  1:1210531/1‑151 (MQ=255)
        cccGGACAACGTCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAAc                                                                                                                                          <  1:4704143/151‑1 (MQ=255)
        cccGGACAACGTCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAAc                                                                                                                                          >  1:4080542/1‑151 (MQ=255)
        accagACAACGTCTATGTCGGCACCAATGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAAc                                                                                                                                          <  2:3777795/147‑1 (MQ=37)
        accagACAACGTCTATGTCGGCACCAATGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAAc                                                                                                                                          <  1:3401554/147‑1 (MQ=37)
         ccGGACAACGTCTATGTCGGCACTAACGCGCCGGTGGTCAATGCCATCGATGCCGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACg                                                                                                                                         >  1:3053096/1‑151 (MQ=255)
         ccGGACAACGTCTATGTCGGCACTAACCAGCCGGTGGGCAATGCCATCGACGCCGTTAGCGGCGCTGATGCGGGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACg                                                                                                                                         <  2:4094078/151‑1 (MQ=37)
              cAACGTCTATGTCGGCACTAACGCGCCGGTGGTCAATGCCATCGATGCCGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGc                                                                                                                                    >  1:4461571/1‑151 (MQ=255)
               aaCGTCTATGTCGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCCGTTAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCg                                                                                                                                   <  2:2432447/151‑1 (MQ=25)
                 cGTCTATGTCGGTACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGCGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGac                                                                                                                                 >  2:2795811/1‑150 (MQ=255)
                 cGTCTATGTCGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGat                                                                                                                                 <  1:382955/151‑1 (MQ=255)
                 cGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGac                                                                                                                                 >  2:4406310/1‑150 (MQ=25)
                   tCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCACCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGacaa                                                                                                                               >  2:1533456/1‑148 (MQ=255)
                   tCTATGTCGGCACCAATGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGGTCGAAAACCTGAACCTGGACCAGACCCCGGCCAGCACCGAGGTCACCGACGAGCGCGGCACCCCGCGCAACGAAGGCCatct                                                                                                                               >  2:3479609/1‑149 (MQ=255)
                   tCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAAACCCCCGTCAGCACCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGatat                                                                                                                               >  1:607723/1‑151 (MQ=255)
                    cTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACatc                                                                                                                              <  1:1124593/151‑1 (MQ=255)
                      aTGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACatcgt                                                                                                                            <  1:3207963/151‑3 (MQ=255)
                      aTGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACatcgt                                                                                                                            >  1:3983548/1‑149 (MQ=255)
                                tAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCAct                                                                                                                  >  2:3909788/1‑150 (MQ=255)
                                 aaCCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCActa                                                                                                                 >  2:2034251/1‑149 (MQ=255)
                                 aaCCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCActa                                                                                                                 <  1:3507324/151‑3 (MQ=255)
                                 aaCCAGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCa                                                                                                                 <  1:3289269/151‑1 (MQ=255)
                                  atgcgCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCAt                                                                                                                >  1:171604/5‑151 (MQ=38)
                                     cgCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCAc                                                                                                             <  2:1162154/150‑1 (MQ=39)
                                     cgCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCAc                                                                                                             >  2:2901965/2‑151 (MQ=38)
                                     cgCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGAGCAACGAAGGCGATCTGGCCAAGGTCACCATCAc                                                                                                             >  1:1005410/2‑151 (MQ=25)
                                     aGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTAc                                                                                                             >  2:2532169/1‑151 (MQ=255)
                                     aGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCAc                                                                                                             >  2:332256/1‑151 (MQ=255)
                                      gCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCActattact                                                                                                            <  1:531114/151‑8 (MQ=255)
                                      gCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCAcc                                                                                                            >  1:4355150/1‑150 (MQ=255)
                                       ccGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTAttactg                                                                                                           >  1:3571727/1‑146 (MQ=38)
                                         ggtggtCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGcc                                                                                                         <  2:3219956/151‑1 (MQ=38)
                                         ggtggtCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGct                                                                                                         >  2:3382747/1‑150 (MQ=255)
                                               aaaggCCAGCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAg                                                                                                   <  2:2330346/147‑1 (MQ=37)
                                                aacgCCATCGACGCGGTCAGCGGCGCAGATGCGGGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACAGCCGACCAGa                                                                                                  <  2:2139837/148‑1 (MQ=25)
                                                 aTGCCATCGACGCGGTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGAc                                                                                                 >  1:4962065/1‑151 (MQ=37)
                                                     cATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCg                                                                                             <  1:680942/151‑1 (MQ=25)
                                                      aTCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCgg                                                                                            >  2:525282/1‑151 (MQ=21)
                                                       tCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGCCCTCGGt                                                                                           >  1:152275/1‑151 (MQ=21)
                                                        cGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTg                                                                                          >  2:4792236/1‑151 (MQ=37)
                                                           cgcgGTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGcc                                                                                       <  1:4859321/151‑2 (MQ=25)
                                                            gccgTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGccg                                                                                      >  2:1064872/4‑149 (MQ=25)
                                                              ggTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCgag                                                                                    >  2:4306997/1‑151 (MQ=25)
                                                              ggTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCgag                                                                                    <  2:739388/151‑1 (MQ=37)
                                                                 cAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAAc                                                                                 <  2:3687097/151‑1 (MQ=21)
                                                                 cAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAAACCCCCGTCAGCACCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAAc                                                                                 <  2:812873/151‑1 (MQ=25)
                                                                 cAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAAc                                                                                 <  2:1620774/151‑1 (MQ=18)
                                                                  aGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACg                                                                                <  1:3655672/151‑1 (MQ=18)
                                                                  aGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACg                                                                                >  2:136267/1‑151 (MQ=25)
                                                                  aGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACg                                                                                <  1:2838531/151‑1 (MQ=18)
                                                                   gcggcgCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGt                                                                               >  2:459230/1‑151 (MQ=18)
                                                                     ggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCa                                                                             <  2:1909859‑M2/92‑1 (MQ=255)
                                                                     ggcgCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGCCa                                                                             >  2:4008246/1‑151 (MQ=25)
                                                                     ggcgCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCa                                                                             >  2:4213416/1‑151 (MQ=21)
                                                                        gcAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTAGCCGAGAACGTCaaac                                                                          >  1:2472089/1‑149 (MQ=14)
                                                                        gcAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCaaac                                                                          >  2:4875849/1‑149 (MQ=12)
                                                                         cAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAAcc                                                                         <  2:2230921/151‑1 (MQ=12)
                                                                          tgATGCATGGAAGTTCGAAAACCTGAACCTGGACAAAACCCCCGTCAGCACCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAAccg                                                                        >  1:2170921/2‑150 (MQ=21)
                                                                           gatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCAGAGAacgtcaaacaga                                                                       >  2:1841241‑M2/54‑146 (MQ=255)
                                                                           gATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAAccga                                                                       <  1:13959/151‑3 (MQ=12)
                                                                           gATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAAccga                                                                       >  1:226170/1‑149 (MQ=17)
                                                                             agcgtGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGAcc                                                                     <  1:2006462/147‑1 (MQ=255)
                                                                               cgtGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCtt                                                                   >  2:4960621/3‑151 (MQ=18)
                                                                                  ggAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAAc                                                                                 <  1:1819842/134‑1 (MQ=25)
                                                                                  ggAAGTTCGAAAACCGGACCCTGGACAGGCCCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAACCCGACCTTCAc                                                                <  2:1454531/151‑1 (MQ=17)
                                                                                   gAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCAcc                                                               >  1:2557977/1‑151 (MQ=21)
                                                                                    aagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCg                                                              >  2:4684900‑M2/45‑151 (MQ=255)
                                                                                        tCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCa                                                          >  2:900382/1‑151 (MQ=2)
                                                                                        tCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCa                                                          <  1:85101/151‑1 (MQ=17)
                                                                                         cGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCAc                                                         >  1:3900126/1‑151 (MQ=2)
                                                                                         cGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCAc                                                         >  1:4273319/1‑151 (MQ=14)
                                                                                          gAAAACCTGAACCTGGACAAAACCCCCGTCAGCACCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACg                                                        <  1:2368779/151‑1 (MQ=18)
                                                                                            aaaCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTc                                                      >  1:2744965/1‑151 (MQ=18)
                                                                                             aacctgaacctgGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCa                                                     >  2:3105012/1‑151 (MQ=2)
                                                                                             aacctgaacctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAAc                                                                                                                                          >  1:2652280/1‑66 (MQ=14)
                                                                                                ctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTGCACATCAaca                                                  <  1:479033‑M2/119‑2 (MQ=255)
                                                                                                 tgaacctgGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGCACGTCAAACCGACCTTCACCGTGCACGTCAACCa                                                 >  1:968566/1‑151 (MQ=17)
                                                                                                  gaacctgGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAg                                                >  1:2954284/1‑151 (MQ=0)
                                                                                                  gaacctgGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAg                                                <  1:2094704/151‑1 (MQ=17)
                                                                                                     cctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGccc                                             >  1:5158430/1‑150 (MQ=14)
                                                                                                      ctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTGCacatcaacaccgccc                                            >  1:4523466‑M2/27‑144 (MQ=255)
                                                                                                                ccccggTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGa                                  >  1:1287116/6‑151 (MQ=17)
                                                                                                                    aGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTg                              <  1:3733384/151‑1 (MQ=14)
                                                                                                                    aGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTg                              >  2:4506250/1‑151 (MQ=14)
                                                                                                                    aGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGACAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTg                              >  2:2719006/1‑151 (MQ=17)
                                                                                                                    aGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTg                              >  2:627824/1‑151 (MQ=14)
                                                                                                                     gTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTga                             <  1:77654/151‑2 (MQ=14)
                                                                                                                       cAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATc                           >  1:4180244/1‑151 (MQ=17)
                                                                                                                            ccGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGAc                      <  2:2619485/151‑1 (MQ=11)
                                                                                                                            ccGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGAc                      <  1:163131/151‑1 (MQ=14)
                                                                                                                              gAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGAccc                    <  1:1913749/151‑1 (MQ=17)
                                                                                                                                ggTGACCGAGGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTg                  >  1:3125074/1‑151 (MQ=14)
                                                                                                                                 gTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGa                 <  1:4616062/151‑1 (MQ=14)
                                                                                                                                       gacgaGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaaca           >  2:3004500/1‑151 (MQ=11)
                                                                                                                                         cgaGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaacaat         >  2:5153241/1‑150 (MQ=14)
                                                                                                                                          gaGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaacaatg        >  2:1751787/1‑149 (MQ=14)
                                                                                                                                            gCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGct      <  2:5062416/151‑2 (MQ=17)
                                                                                                                                                     |                                                                                                                                                 
GCAACCGCCCCGGACAACGTCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGG  >  CP009974/1897445‑1897739

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