Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,897,666 G→A intergenic (‑5444/‑3153) 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,6660GA77.5% 2886.6 / inf 2639intergenic (‑5444/‑3153)RPPX_08470/RPPX_08480tryptophan synthase subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (307/285);  new base A (1026/1009);  total (1339/1300)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-11

Only 100 of 2646 total aligned reads displayed.

CAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCACAGGGCGATGACGTCTACCAGGACGCAGGC  >  CP009974/1897518‑1897816
                                                                                                                                                    |                                                                                                                                                      
      cATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCtt                                                                                                                                                >  1:4183366/1‑151 (MQ=9)
       gtGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTc                                                                                                                                               <  1:7217185/150‑1 (MQ=21)
       gtGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTc                                                                                                                                               <  2:6574915/150‑1 (MQ=25)
       gtGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTc                                                                                                                                               >  2:2652223/2‑151 (MQ=25)
       gtGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTc                                                                                                                                               <  2:4634429/150‑1 (MQ=17)
       gtGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTc                                                                                                                                               >  1:5378761/2‑151 (MQ=17)
       aTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGCCAAACCGACCTTc                                                                                                                                               >  2:6463157/1‑151 (MQ=12)
        tGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTcc                                                                                                                                              >  2:751997/1‑150 (MQ=17)
           aaGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAAGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCg                                                                                                                                           <  1:78246/151‑1 (MQ=18)
           aaGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCg                                                                                                                                           >  1:6878052/1‑151 (MQ=17)
            aGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCActgt                                                                                                                                          >  1:4315899/1‑148 (MQ=9)
             gttcgaacagctggcgctggacaaggcgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTg                                                                                                                                         <  2:1226754‑M2/109‑1 (MQ=255)
             gTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTg                                                                                                                                         <  1:6343549/151‑1 (MQ=18)
              ttCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGc                                                                                                                                        <  2:5397454/151‑1 (MQ=255)
                cGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCAc                                                                                                                                      >  2:3521879/1‑151 (MQ=14)
                cGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCAc                                                                                                                                      >  2:2545757/1‑151 (MQ=17)
                 gAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACg                                                                                                                                     >  1:2199783/1‑151 (MQ=14)
                 gAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACg                                                                                                                                     >  1:4614255/1‑151 (MQ=21)
                 gAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACg                                                                                                                                     >  2:3892191/1‑151 (MQ=18)
                    accctgaccctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCa                                                                                                                                  <  2:5198755/149‑1 (MQ=14)
                    aacctgaacctgGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGCACGTCAAACCGACCTTCACCGTGCACGTCa                                                                                                                                  >  2:7030519/1‑151 (MQ=2)
                     gcctgaacctgGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCaa                                                                                                                                 <  2:4475126/150‑1 (MQ=2)
                     acctgaacctgGACAAGGCCCCGGTGAGCGCCGAGGTGACCGACGAGCTAGGTACGCCGGGCAGCGAAGGCGATATGGACAAGGTCACCACCACGGCCGACCAGACCTCGGTAGCCGAGAACGACAAACCGACCTTCACCGTGCAcctaac                                                                                                                                 >  2:2263339/1‑146 (MQ=14)
                     acctgaacctgGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCaa                                                                                                                                 <  2:3794301/151‑1 (MQ=17)
                      cctgaacctgGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAAc                                                                                                                                >  2:150796/1‑151 (MQ=11)
                       ctgaacctgGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAAcc                                                                                                                               <  1:2764624/151‑1 (MQ=2)
                          aacctgGACAAAACCCCCGTCAGCACCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGc                                                                                                                            >  1:9068/1‑151 (MQ=17)
                           cgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGGGGCCGAGAACGTCAAACCGACGTTCACCGTGCACatcaacaccgc                                                                                                                           >  2:4882405‑M2/29‑146 (MQ=255)
                           acctgGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGcc                                                                                                                           <  2:873927/151‑1 (MQ=11)
                           acctgGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACTAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGcc                                                                                                                           <  2:6882295/151‑1 (MQ=17)
                            cctgGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTTGCACCCCAGTCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTTACCAGACTTCTGTAGCCGAGAACGCCCAACCGACCTTCACTGTGCACGTCAACCaccac                                                                                                                          >  2:5023465/1‑147 (MQ=255)
                            cctgGAACAGACCCCCGTCAGCACCGAGGTGCCCGACGAGCCACGCACCCTGGGCAACGAAGGCGATATCGTCAAGGTCACCATCGCGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCg                                                                                                                          <  2:4747453/151‑1 (MQ=14)
                               ggACAAGCCCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTg                                                                                                                       <  2:6737954/151‑1 (MQ=11)
                               ggACAAGACTCCAGTCAGCACCGAGGGCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCCTCAAGGTCACCATCACCGCCGCCCAGCCCTCGGTCGCTGAGCACGGCAACCCCACCTCCACCGCGCACGCCCGCCAGCCGcgc                                                                                                                       >  2:7127673/1‑149 (MQ=255)
                               ggACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTg                                                                                                                       >  1:3330274/1‑151 (MQ=11)
                               ggACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTg                                                                                                                       >  1:3877538/1‑151 (MQ=14)
                                   aaGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGccc                                                                                                                   >  1:6506271/1‑151 (MQ=17)
                                   aaGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGctc                                                                                                                   <  1:1002582/151‑3 (MQ=11)
                                    aGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCa                                                                                                                  >  1:5458316/1‑151 (MQ=11)
                                       ccccggTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGa                                                                                                               >  2:5548460/6‑151 (MQ=11)
                                        ccccggcAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGAc                                                                                                              <  2:1355849/145‑1 (MQ=14)
                                         ccAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGAcc                                                                                                             <  1:4042577/151‑1 (MQ=14)
                                          cAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGActt                                                                                                            <  2:4185449/151‑3 (MQ=18)
                                           aGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTg                                                                                                           >  2:6035717/1‑151 (MQ=14)
                                            gTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTga                                                                                                          >  1:3221309/1‑150 (MQ=14)
                                            gTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTgg                                                                                                          >  1:237941/1‑151 (MQ=14)
                                             tCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTgat                                                                                                         >  1:7118155/1‑149 (MQ=17)
                                             tCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGActcggt                                                                                                         >  2:3284748/1‑146 (MQ=17)
                                             tCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGt                                                                                                         <  1:2042395/151‑1 (MQ=17)
                                              cAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTc                                                                                                        <  1:6179724/151‑1 (MQ=11)
                                              cAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTc                                                                                                        <  1:3058348/151‑1 (MQ=11)
                                                gCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGt                                                                                                      <  1:4437232/151‑1 (MQ=14)
                                                gCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGt                                                                                                      >  1:161352/1‑151 (MQ=17)
                                                   ccGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGAc                                                                                                   >  2:2146355/1‑151 (MQ=14)
                                                    cGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGAcc                                                                                                  <  1:3740840/151‑1 (MQ=17)
                                                    cGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGAcc                                                                                                  >  2:2840137/1‑151 (MQ=17)
                                                     gaggagACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGAccc                                                                                                 >  2:4029535/1‑151 (MQ=11)
                                                        gtcaCCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                              >  2:1116305/4‑151 (MQ=11)
                                                        gtaaCCTACGAGCCAGGCACCGCGGGCAAGGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGAGCGTGACCCTGa                                                                                              <  2:6155169/148‑1 (MQ=11)
                                                        gTGACCGACGAGCCAGGCACCCCGAGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGa                                                                                              <  1:1166422/151‑1 (MQ=14)
                                                         tGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAg                                                                                             >  2:5299121/1‑151 (MQ=1)
                                                            ccGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGcaa                                                                                          <  1:5059931/151‑1 (MQ=17)
                                                               acgaGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaacaa                                                                                       >  1:4772684/1‑151 (MQ=14)
                                                                 gaGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACg                                                                                     <  1:5987714/151‑1 (MQ=255)
                                                                 gaGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCGCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaacaatg                                                                                     <  2:1843153/151‑3 (MQ=14)
                                                                    ccAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaacaatgctc                                                                                  >  1:4534477/1‑146 (MQ=14)
                                                                    ccAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGctc                                                                                  >  2:3408366/1‑149 (MQ=14)
                                                                     cAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAAGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCa                                                                                 >  2:3129124/1‑151 (MQ=1)
                                                                     cAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCa                                                                                 <  1:5164094/151‑1 (MQ=17)
                                                                      aGGTACGCCGGGCAACGAAGGCGACTTGGTCAAAGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAg                                                                                >  1:5833932/1‑151 (MQ=11)
                                                                      aGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAg                                                                                >  1:3469388/1‑151 (MQ=17)
                                                                       ggTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAgg                                                                               >  2:2995681/1‑151 (MQ=17)
                                                                              caggCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCAt                                                                        <  2:6311674/149‑1 (MQ=11)
                                                                               gggCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCCCAGGTCACCATc                                                                       <  2:5899258/151‑1 (MQ=11)
                                                                                 gCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCCCAGGTCACCATCaa                                                                     >  1:1284761/1‑151 (MQ=11)
                                                                                 gCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCaa                                                                     >  1:2637943/1‑151 (MQ=11)
                                                                                  cAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAg                                                                    <  1:5254558/151‑1 (MQ=14)
                                                                                     cGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGct                                                                 <  2:4315808/151‑2 (MQ=255)
                                                                                        aGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTggc                                                              >  1:6742724/1‑150 (MQ=11)
                                                                                         ggCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCt                                                                                                                        <  2:1982289/92‑1 (MQ=14)
                                                                                          gCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGCCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAAGCCggcca                                                            >  2:5810788/1‑148 (MQ=11)
                                                                                           cGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGctggcgaa                                                           >  1:6083410/1‑144 (MQ=11)
                                                                                            gacctggtcaAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGaaa                                                          >  1:1467320/8‑151 (MQ=11)
                                                                                              tatCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAAcc                                                        >  1:6591550/1‑151 (MQ=11)
                                                                                                tggtcaAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAg                                                      >  2:2411080/4‑151 (MQ=14)
                                                                                                    caAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAgtgcg                                                  >  1:2930579/1‑151 (MQ=11)
                                                                                                     aaGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAgtgcgc                                                 >  1:5604685/1‑151 (MQ=255)
                                                                                                     aaGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAgtgcgc                                                 >  1:5731772/1‑151 (MQ=255)
                                                                                                          cACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGTGCGCCGTac                                            >  1:361904/1‑151 (MQ=11)
                                                                                                           aCCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTaca                                           <  1:2500200/151‑1 (MQ=255)
                                                                                                            ccATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGTGCGCCGTacac                                          >  1:991708/1‑151 (MQ=255)
                                                                                                               tCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCa                                       <  2:6387461/151‑1 (MQ=11)
                                                                                                                        aCCAGACTTCGGTGGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGTGCGCCGTACACCCACGACGcgca                              <  1:6988852/151‑1 (MQ=11)
                                                                                                                            gACTTCGGTGGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGTGCGCCGTACACCCACGACGCGCAAGGc                          <  1:1935609/151‑1 (MQ=11)
                                                                                                                              cttCGGTGGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGTGCGCCGTACACCCACGACGCGCAAGGCGa                        >  2:1657527/3‑151 (MQ=11)
                                                                                                                                tCGGTGGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATg                      <  1:1731523/151‑1 (MQ=11)
                                                                                                                                   gTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCACAGGGCGATGACg                   >  1:645182/1‑151 (MQ=11)
                                                                                                                                    tGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCACAGGGCGATGACGt                  >  2:3664172/1‑151 (MQ=11)
                                                                                                                                      gccgagaACCTCAAGCCAACCTTCACCGTGCGCGTCATCCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGCGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCt                <  2:5681209/148‑1 (MQ=11)
                                                                                                                                          agaACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCACAGGGCGATGACGTCTACCa            >  1:4924539/1‑151 (MQ=11)
                                                                                                                                                    |                                                                                                                                                      
CAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCACAGGGCGATGACGTCTACCAGGACGCAGGC  >  CP009974/1897518‑1897816

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: