Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 3,837,799 A→G intergenic (+350/‑201) RPPX_17095 → / → RPPX_17100 proteophosphoglycan precursor/cupin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099743,837,7990AG97.9% 546.0 / ‑4.4 189intergenic (+350/‑201)RPPX_17095/RPPX_17100proteophosphoglycan precursor/cupin
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (87/96);  total (88/101)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 189 total aligned reads displayed.

CTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGCCATTT  >  CP009974/3837704‑3837946
                                                                                               |                                                                                                                                                   
cTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCggtggt                                                                                              >  1:1296588/1‑151 (MQ=2)
 tACGCACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGtggtt                                                                                             <  1:443423/151‑1 (MQ=255)
   cacacGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGtggttgg                                                                                           >  1:4181419/1‑151 (MQ=9)
     cacGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTc                                                                                         >  2:5681839/1‑151 (MQ=255)
       cGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGCTGCCGCCGCTCCCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTcac                                                                                       >  2:856077/1‑151 (MQ=11)
         cATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAg                                                                                     >  2:318653/1‑151 (MQ=255)
            cGAGTGTGGCACAGGCCGTTGTAGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGcc                                                                                  <  1:6055884/151‑1 (MQ=14)
            cGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGcc                                                                                  <  1:4853792/151‑1 (MQ=255)
            cGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGcc                                                                                  <  1:5235674/151‑1 (MQ=255)
            cGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGcc                                                                                  >  1:2326002/1‑151 (MQ=255)
             gAGGGCGGTACAGGCCGTTGTAGGCGCGGTTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcg                                                                                 <  1:2319521/151‑1 (MQ=14)
              aGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgt                                                                                >  1:3224119/1‑151 (MQ=255)
              aGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgt                                                                                >  1:3223922/1‑151 (MQ=255)
              aGGGCAGCCCCGGCCGTTGCAGCAGCGGGTCCACCCGCGAAAGGNCCGGGCCTGCTGGCCAATATCNACCCGCCGTTCACCATTGACGACGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTAACAGACcnt                                                                                <  2:3166791/151‑1 (MQ=14)
               gggCCGCACAGGCCGTTCTAGGAGCTGCTCCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgtc                                                                               <  1:6533456/151‑1 (MQ=14)
                ggcggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCCCACGCATCGGTGGTTGGTCACCGGCcgtcg                                                                              >  1:5276952/1‑151 (MQ=14)
                ggcggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgtcg                                                                              >  1:1785015/1‑151 (MQ=255)
                 gcggcACAGGCCGTTGTAGGACCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgtcgt                                                                             <  1:6426518/151‑1 (MQ=255)
                  gggcgtaGGCCGTGGTAGGGCCGGGGCAACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTa                                                                            <  2:3067230/145‑1 (MQ=14)
                  cgacacaGGCCGTCCTCGACGCGGGTTCACCCGCGAAAGGGCCGGGCCGGCTGGCCAATCTCCACCCGCCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTa                                                                            <  2:2183577/148‑1 (MQ=14)
                   ggcAGAGGCGGGTGGAGGGGCGGGTGCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAg                                                                           <  2:3357692/151‑1 (MQ=14)
                   ggcACAGGCCGTTGTAGGAGCGGGTTCCCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAg                                                                           <  1:569463/151‑1 (MQ=255)
                   ggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAg                                                                           <  1:942825/151‑1 (MQ=255)
                   ggcACAGGCCGTGGAGGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAg                                                                           <  2:6321615/151‑1 (MQ=14)
                     cacaGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGa                                                                         <  1:1349423/151‑1 (MQ=255)
                      acaGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACCCGCATCGGTGGTTGGCCACAGCCCGTCGTAgcac                                                                        >  1:6070119/1‑148 (MQ=11)
                      acaGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAg                                                                        >  2:4070480/1‑151 (MQ=255)
                       caGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGc                                                                       >  1:109654/1‑151 (MQ=255)
                       caGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGc                                                                       <  1:6441413/151‑1 (MQ=255)
                        gggCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCg                                                                      <  2:466786/150‑1 (MQ=255)
                         ggCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCgg                                                                     >  1:4044593/1‑151 (MQ=255)
                         ggCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCgg                                                                     >  1:3457410/1‑151 (MQ=255)
                          gCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCggg                                                                    >  1:1944493/1‑151 (MQ=255)
                           ccGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGt                                                                   >  1:844766/1‑151 (MQ=255)
                              ttGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCa                                                                >  1:3612585/1‑151 (MQ=255)
                                gTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCAcc                                                              >  1:6072588/1‑151 (MQ=255)
                                   ggggggggTTCCCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                                           <  1:5202090/146‑1 (MQ=14)
                                   ggAGCGGGTGCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                                           <  1:4063440/151‑1 (MQ=255)
                                     aGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGa                                                         <  1:5771594/151‑1 (MQ=255)
                                     aGCGGGTTAACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGa                                                         <  2:2832942/151‑1 (MQ=255)
                                       ccggtaaaCCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTGGGTCACAGGCCGTCGTAGGAGCGGGCCCCCGCGCGaaa                                                       >  2:6215275/8‑151 (MQ=14)
                                       ccggTTCCCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGaaa                                                       <  2:1331707/149‑1 (MQ=255)
                                         ggTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCCCCCGCGAAAAg                                                     >  1:6161332/1‑151 (MQ=255)
                                         ggTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAg                                                     <  2:4736649/151‑1 (MQ=255)
                                          gTTTACCCGCGAAGGGGCGGGGCCTGCGGGGGAGTTTCCTCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTGCGCCCGCTCCTGCGCACGCGGCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAgg                                                    <  2:1951760/151‑1 (MQ=14)
                                          gTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAgg                                                    <  1:4479517/151‑1 (MQ=255)
                                           ttCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTGACAGGCCGCCGTAGGAGCGGGCTCACCCGCCAAAAGGc                                                   >  2:2001165/1‑151 (MQ=14)
                                           ttCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGc                                                   <  2:3799721/151‑1 (MQ=255)
                                           ttCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGc                                                   >  1:2401865/1‑151 (MQ=255)
                                             cACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCg                                                 >  1:388468/1‑151 (MQ=255)
                                                ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTa                                              <  1:616227/151‑1 (MQ=255)
                                                 ggggggggggcggggcgTGCTGGCCCGTTTCGACGCGCCGTTCAGCGTTGCCGCCGCTGCCGTGTGCGCGGGGAGACCGGCGCCTGCGCACGCGTCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCGCCCGCGAAAAGGCCGGTAc                                             <  2:5968768/144‑1 (MQ=9)
                                                 cgcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTAc                                             <  1:1008616/151‑1 (MQ=255)
                                                  gcgAAAGGGCCGGGCCTGCTGGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACa                                            <  1:260151/151‑1 (MQ=255)
                                                  gcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACa                                            <  1:516493/151‑1 (MQ=255)
                                                    gAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAgg                                          <  1:3790693/151‑1 (MQ=255)
                                                     aaaGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGc                                         <  1:5910525/151‑1 (MQ=255)
                                                      aaGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCg                                        <  1:3572442/151‑1 (MQ=255)
                                                      aaGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCg                                        >  1:5564744/1‑151 (MQ=255)
                                                       gggggcgggccTGCTCGCCACTCTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCGCCTGCACACGCATCGGGGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGa                                       <  2:5209766/150‑1 (MQ=14)
                                                          gccgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGaaaa                                    >  1:1871321/1‑151 (MQ=255)
                                                           ccgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAAt                                   <  2:5410206/151‑1 (MQ=255)
                                                             gggccTGCTGGCCAATGTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCg                                 <  2:1771426/151‑1 (MQ=255)
                                                              ggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCgg                                <  2:4979476/151‑1 (MQ=255)
                                                              ggccGGCTGGCCAAGCTCCCTCCGCCCTTCACCGGTGCCGTCGCTTCCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCgg                                <  2:4992878/151‑1 (MQ=17)
                                                               gccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGc                               >  2:1149627/1‑151 (MQ=255)
                                                                ccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGcc                              >  1:755685/1‑151 (MQ=255)
                                                                 ctgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCTGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccc                             >  1:1804495/1‑151 (MQ=255)
                                                                 ctgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccc                             >  1:6267790/1‑151 (MQ=255)
                                                                 ctgctgGCCAATTTCCACCCGCCGTCCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccc                             <  1:4284869/151‑1 (MQ=255)
                                                                 cggctcGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccc                             <  2:900852/149‑1 (MQ=255)
                                                                  tgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACCCATCGGTGGTTGGGCACAGGCCGTCGGAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACCGGCGAAAATCGGcccc                            >  2:4000327/1‑151 (MQ=14)
                                                                   gctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTCCCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCCAACCTCGGccccc                           >  2:3641270/1‑151 (MQ=14)
                                                                    ctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCt                          <  1:2850739/151‑1 (MQ=255)
                                                                     tgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTg                         <  1:1888209/151‑1 (MQ=255)
                                                                     tgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTg                         >  2:1485434/1‑151 (MQ=255)
                                                                      ggCCAATCTCCACCCGTCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGGGGCTGGTCACAGGCCGTCGGAGGCGCGGGCTCCCCCGCGAAGAGCCCGGCACAGGCGAAAATCGGCCCCCTGc                        >  1:3971196/1‑151 (MQ=18)
                                                                       gCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAg                                                                           <  1:604282/99‑1 (MQ=14)
                                                                         cAATTTCCGCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCttt                     <  2:4767074/151‑1 (MQ=255)
                                                                         cAATTTCCACCCGCCGTTCACCGTTGCCGCGGCTACCGTATTCGCGCGTAAACCCGCTCCTGCACACGCATCGGGGGCTGGTCGCAGGCCGTCGTAGGCGCGGGGTCCCCCGCGAAAAGGCCGGTACAGGCGCACATCGGCCCCCTGCgtt                     >  2:3905430/1‑151 (MQ=17)
                                                                          aaTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCtttt                    <  1:1567266/151‑1 (MQ=255)
                                                                          aaTCGCCACCGGCAGTTCACCTTTGCCGCCGCGAGCGTATGCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCtttt                    <  2:4551650/151‑1 (MQ=18)
                                                                           ctttcCACCCGCCGTTCACACTTGCGCCCGCTACCGAATTCGCGGGTAAACCCGCTCCTGCACCCGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTc                   <  2:1042267/148‑1 (MQ=21)
                                                                           atttcCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTc                   <  2:785217/148‑1 (MQ=255)
                                                                            tttcCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTcc                  >  1:5357221/3‑151 (MQ=255)
                                                                                cACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCGCACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCACCCTGCTTTTCCCCCg              >  2:2490648/1‑151 (MQ=255)
                                                                                   ccgccgTTCCCCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAt           <  2:4308224/151‑1 (MQ=255)
                                                                                    cgccgtTCTCCGTTGCCGCCGCGGCCGTATTCGCGGGTAACCCCGCTCCTGCACACGCAGCGGTGGTGGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAtt          <  2:6070119/148‑1 (MQ=21)
                                                                                    cgccgtTCACCGTTGCCGCCGCTACCGTGTTCGCGGGTAGACCCGCTCCTGCACCCGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAtt          <  1:1198386/148‑1 (MQ=255)
                                                                                    cgccgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAGACCCGCCCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAtt          <  1:6369452/148‑1 (MQ=255)
                                                                                    cgccgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAtt          <  1:2790192/148‑1 (MQ=255)
                                                                                    cgccgtTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAtt          <  1:6215275/148‑1 (MQ=255)
                                                                                     gccttgcACCGTTGCCGCCGCTACCGTATTCGCCGGTAAACCCGGTCCTGCGCACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGCGCGGGTTCCCCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTg         <  2:4038987/145‑1 (MQ=25)
                                                                                     gccgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTg         <  1:1013241/149‑1 (MQ=255)
                                                                                       cgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGc       >  2:1724700/1‑151 (MQ=255)
                                                                                        gtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGcc      <  2:2335469/151‑1 (MQ=255)
                                                                                        gtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGcc      <  2:684293/151‑1 (MQ=255)
                                                                                        gtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGcc      <  1:4754998/151‑1 (MQ=255)
                                                                                         ttCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGCCa     <  2:556231/151‑1 (MQ=255)
                                                                                           cACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGCTTGGCCAtt   >  1:5373472/1‑151 (MQ=255)
                                                                                               |                                                                                                                                                   
CTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGCCATTT  >  CP009974/3837704‑3837946

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