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.1% 418.2 / 0.0 137intergenic (+350/‑201)RPPX_17095/RPPX_17100proteophosphoglycan precursor/cupin
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (78/55);  total (80/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 137 total aligned reads displayed.

CTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGCCATTT  >  CP009974/3837704‑3837946
                                                                                               |                                                                                                                                                   
cTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTCTTCCCGGGTAAACCCGCTCCAGCACACGCATCggtggt                                                                                              >  1:1326978/1‑151 (MQ=11)
   cacacGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGtggttgg                                                                                           >  1:5047004/1‑151 (MQ=9)
    acacGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGtggttggt                                                                                          <  1:1540499/151‑1 (MQ=255)
    acacGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGtggttggt                                                                                          >  2:4655247/1‑151 (MQ=255)
    acacGCACCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGtggttggt                                                                                          <  1:1541744/151‑1 (MQ=255)
     cacGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCTCACGCAGCGGGGGTTGGTc                                                                                         >  1:2670419/1‑151 (MQ=11)
       cGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACAGTTGCCGCCGCTACCGTATTCGCGGGTGAACCAGCTCCTGGACACGCATCGGTGGCTGGTcca                                                                                       >  2:3136998/1‑149 (MQ=11)
        gCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTcaca                                                                                      >  2:612151/1‑151 (MQ=255)
        gCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTcaca                                                                                      <  1:2664202/151‑1 (MQ=255)
         cATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAg                                                                                     >  2:1404092/1‑151 (MQ=255)
          aTCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCCTCGGGGGTTGGTCACAgg                                                                                    >  1:4154472/1‑151 (MQ=11)
             gAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcg                                                                                 <  1:2746024/151‑1 (MQ=255)
              aGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgt                                                                                >  2:4352168/1‑151 (MQ=255)
               gggCGGCGCAGGCCGTGGTAGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTAGTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgtc                                                                               <  2:4955671/151‑1 (MQ=14)
                  cggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTa                                                                            >  1:1106779/1‑151 (MQ=255)
                  cggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTa                                                                            <  1:4766392/151‑1 (MQ=255)
                   ggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCCCCCGCTCACGCCCCGGTGGCTGGGCACAGGCCGTCGCAg                                                                           >  2:195752/1‑151 (MQ=11)
                    gcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAgg                                                                          >  2:842508/1‑151 (MQ=255)
                     cacaGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGa                                                                         >  1:3362563/1‑151 (MQ=255)
                      acaGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGCAg                                                                        >  2:1561242/1‑151 (MQ=255)
                        aGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCg                                                                      >  1:2060636/1‑151 (MQ=255)
                         ggCCGTTGTAGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCgg                                                                     <  2:1981728/151‑1 (MQ=14)
                         ggCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCgg                                                                     >  2:415579/1‑151 (MQ=255)
                         ggCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCgg                                                                     <  2:3095580/151‑1 (MQ=255)
                            cGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGtt                                                                  >  2:1662761/1‑151 (MQ=255)
                                ttAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCAcc                                                              <  1:4482402/150‑1 (MQ=255)
                                gTAGGCGCGGGTTGACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCAcc                                                              <  2:1052279/151‑1 (MQ=255)
                                 tAGGAGCTGGTGCACCCGCGAAAGGGCCGGGCCTGCGTGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCGACCGTATTCGCGGGTAATCCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCAccc                                                             <  2:590134/151‑1 (MQ=14)
                                  aGGAGCGGGTACACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcg                                                            <  2:332421/151‑1 (MQ=255)
                                   ggAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                                           >  1:1933904/1‑151 (MQ=255)
                                     aGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGa                                                         >  2:4539898/1‑151 (MQ=255)
                                     aGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGa                                                         <  2:3737219/151‑1 (MQ=255)
                                      gCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGCGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCCaa                                                        >  2:2661177/1‑151 (MQ=14)
                                        gggTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGaaaa                                                      >  1:1727971/1‑151 (MQ=255)
                                         ggtaaaCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAg                                                     <  2:1462053/146‑1 (MQ=14)
                                         ggTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAg                                                     >  1:1159365/1‑151 (MQ=255)
                                          gtggagccGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGACGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAgg                                                    <  2:2263540/145‑1 (MQ=14)
                                           ttCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGc                                                   >  1:1860971/1‑151 (MQ=255)
                                               cccGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGt                                               >  1:3900569/1‑151 (MQ=255)
                                                 cgcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGCGCGGGTTCACCCGCGAACAGGCCGGTAc                                             >  2:4709308/1‑151 (MQ=14)
                                                 cgcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTAc                                             <  1:1175574/151‑1 (MQ=255)
                                                  gcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACa                                            >  1:782524/1‑151 (MQ=255)
                                                    gAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAgg                                          <  1:3086456/151‑1 (MQ=255)
                                                    gAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                                           >  1:515672/1‑134 (MQ=255)
                                                    gAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                                           <  2:515672/134‑1 (MQ=255)
                                                     caaGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGc                                         <  1:3507664/150‑1 (MQ=255)
                                                     aaaGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGc                                         <  1:4418896/151‑1 (MQ=255)
                                                         ggccgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGCAAAGGCCGGTACAGGCGaaa                                     >  2:4558870/1‑151 (MQ=255)
                                                          gccgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGaaaa                                    >  1:5293/1‑151 (MQ=255)
                                                           ccgggccTGCTGTCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGTCAACCCCGCCCCTGCACACGCATCGCTCGTCGGTCACAGGCCGTCGTAGGACCTGGTTCACCCGCGAAACGGCCCGTCCCAGCGAAAAt                                   >  2:3755987/1‑151 (MQ=14)
                                                           ccgggccTGCTGGGCAATCTCGACTCGTCGTTCCACGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAAt                                   <  2:3986931/151‑1 (MQ=18)
                                                           ccgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGACGAAAAt                                   >  2:3218224/1‑151 (MQ=255)
                                                            cgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGCAAATc                                  >  2:3334103/1‑151 (MQ=255)
                                                            cgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATc                                  <  1:4349891/151‑1 (MQ=255)
                                                            cgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATc                                  >  1:2082343/1‑151 (MQ=255)
                                                             gggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCg                                 <  1:2226269/151‑1 (MQ=255)
                                                             gcgccTGCTGGCCAATTTCCACCCTCGGTTCACCGTCGCCGCCGCTACCGTATTCGCGGGTGAACCCGCGCCTGCACACGCATCGGTGGTGGGTCACAGGCCGTGGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCg                                 <  1:2226318/149‑1 (MQ=17)
                                                              ggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCgg                                >  2:3653819/1‑151 (MQ=255)
                                                               gccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGc                               >  1:2711800/1‑151 (MQ=255)
                                                               gccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGc                               <  1:2827459/151‑1 (MQ=255)
                                                               gccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGCTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGc                               >  2:4013095/1‑151 (MQ=255)
                                                                ccTGCTGGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGcc                              <  1:319941/151‑1 (MQ=255)
                                                                ccTGCGGGCCAATTTCCCCCCGCGGCTCACCGTTGCGGCGGCGACCGTGTTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTGGGGCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGcc                              <  2:1396955/151‑1 (MQ=14)
                                                                 ctgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccc                             >  1:892767/1‑151 (MQ=255)
                                                                  tgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGcccc                            >  1:4865977/1‑151 (MQ=255)
                                                                   gctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCCTCGGTGGTCGGTCACAGGCCGTCGTAGGAGCCGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAACCGGccccc                           >  2:4365754/1‑151 (MQ=17)
                                                                   gctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccccc                           <  1:1869288/151‑1 (MQ=255)
                                                                    ctgGCCACTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCt                          <  2:1518846/151‑1 (MQ=255)
                                                                    cgggCCAATTTCCCACCGACGTTCACCGTTGCCTCCGCTATCGTATTCGCGGGTAAACCCGCTCCTGCACACGCTTCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCt                          <  2:2316928/149‑1 (MQ=18)
                                                                     tgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTg                         >  2:4757691/1‑151 (MQ=255)
                                                                     tgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTg                         <  1:821614/151‑1 (MQ=255)
                                                                     tgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTg                         >  2:3407229/1‑151 (MQ=255)
                                                                      ggCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGGTACCGTATTGGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGGGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGc                        <  1:1949510/151‑1 (MQ=18)
                                                                      ggCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGc                        <  1:4872419/151‑1 (MQ=255)
                                                                         cAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCttt                     <  2:4875736/151‑1 (MQ=255)
                                                                          aaTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGTTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCtttt                    <  1:2899447/151‑1 (MQ=255)
                                                                          aaTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCtttt                    <  2:724072/151‑1 (MQ=255)
                                                                          aaTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCtttt                    >  1:3090902/1‑151 (MQ=255)
                                                                          aaTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCCCTTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGACAATCGGCCCCCTGCtttt                    >  1:1639820/1‑151 (MQ=18)
                                                                           atttcCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTc                   <  2:2729628/148‑1 (MQ=255)
                                                                            tttcCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTcc                  <  1:2367061/149‑1 (MQ=255)
                                                                            tttcCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTcc                  >  1:4143456/3‑151 (MQ=255)
                                                                                 ccccGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCgg             <  1:3588564/150‑1 (MQ=255)
                                                                                 aCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCgg             >  1:321053/1‑151 (MQ=255)
                                                                                 aCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCgg             >  2:1331891/1‑151 (MQ=255)
                                                                                   ccgccgTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGGACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAt           >  1:903283/1‑151 (MQ=255)
                                                                                     gccgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTg         >  2:231242/3‑151 (MQ=255)
                                                                                     gccgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTg         <  2:2947029/149‑1 (MQ=255)
                                                                                     gccgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTg         <  2:2943959/149‑1 (MQ=255)
                                                                                     gccgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTg         >  2:231690/3‑151 (MQ=255)
                                                                                      ccgtGCACCGTTGCGGCCGGTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACCGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTgg        <  1:3924857/150‑1 (MQ=21)
                                                                                       cgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGc       >  2:1139239/1‑151 (MQ=255)
                                                                                        gtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGcc      >  1:1685794/1‑151 (MQ=255)
                                                                                        gtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGcc      <  1:385451/151‑1 (MQ=255)
                                                                                        gtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGACTGGcc      >  2:2890385/1‑151 (MQ=255)
                                                                                         ttCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGCCa     >  2:2274871/1‑151 (MQ=255)
                                                                                         ttCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCCTAGGAGCGGGCTCACCCGCGAAAAGGCCGGTGCAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGccc     >  2:1129489/1‑150 (MQ=255)
                                                                                           cACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGCCAtt   >  1:2791781/1‑151 (MQ=255)
                                                                                           cACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTc                                                                 >  1:2407200/1‑89 (MQ=255)
                                                                                            aCCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGCTTGGCCCttt  >  2:638329/1‑151 (MQ=255)
                                                                                               |                                                                                                                                                   
CTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGCCATTT  >  CP009974/3837704‑3837946

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