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.4% 690.8 / 5.9 227intergenic (+350/‑201)RPPX_17095/RPPX_17100proteophosphoglycan precursor/cupin
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base G (109/111);  total (114/113)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 227 total aligned reads displayed.

CTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGCCATTT  >  CP009974/3837704‑3837946
                                                                                               |                                                                                                                                                   
cTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCCACTGCACACGCATCggtggt                                                                                              >  1:4923199/1‑151 (MQ=11)
 tACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGtggtt                                                                                             <  1:3403249/151‑1 (MQ=9)
      acGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTca                                                                                        >  1:189984/1‑151 (MQ=255)
        gCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTcaca                                                                                      <  1:6128858/151‑1 (MQ=255)
         cATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAg                                                                                     >  1:5564670/1‑151 (MQ=255)
          aTCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAgg                                                                                    >  2:899719/1‑151 (MQ=255)
          aTCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAgg                                                                                    <  1:1809876/151‑1 (MQ=255)
           tCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGc                                                                                   <  1:419778/151‑1 (MQ=255)
           tCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGCTGGTCACAgtc                                                                                   >  2:1744445/1‑149 (MQ=11)
            ggAGGGCGCCGCAGGCCGTAGTAGTAGCCGTCCCATCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGcc                                                                                  <  1:2330484/150‑1 (MQ=14)
            cGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGcc                                                                                  <  1:6636268/151‑1 (MQ=255)
              ggggCGGCACAGGCCGTTGTAGGACCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgt                                                                                <  1:6525588/150‑1 (MQ=255)
              cgggCGGCGCAGGCCGTTGTAGGGGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgt                                                                                <  1:1188139/150‑1 (MQ=14)
              agtgcgacAGAGGCCGGTGTAGTCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgt                                                                                <  1:3116938/148‑1 (MQ=14)
              aGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACCCATCGGCGGTTGGTCACAGGCcgc                                                                                >  1:2800095/1‑150 (MQ=11)
               tggcggcACATGCCGTGCCCGCCGCGCGCTCCCCTGCGACCGGGCCCGGCCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgtc                                                                               <  2:4326962/150‑1 (MQ=14)
               gggCGGCAGAGGCCGTTGTAGGAGCGGGTGCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgtc                                                                               <  2:1916931/151‑1 (MQ=14)
                ggcggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgtcg                                                                              >  2:7018823/1‑151 (MQ=255)
                ggcggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcggcg                                                                              >  2:7019126/1‑151 (MQ=255)
                ggcggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGAGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgtcg                                                                              >  2:6949252/1‑151 (MQ=255)
                  cggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTTACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTa                                                                            <  1:701954/151‑1 (MQ=14)
                   ggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCCTCGGCGGCTGGTCACAGGCcctcggac                                                                           >  2:6718605/1‑148 (MQ=11)
                       caGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGc                                                                       >  2:1855154/1‑151 (MQ=255)
                        aGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCg                                                                      >  1:898553/1‑151 (MQ=255)
                        aGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCg                                                                      <  1:5530906/151‑1 (MQ=255)
                          gCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGGCGTCGTAGGAGCggg                                                                    >  2:96765/1‑151 (MQ=255)
                          gCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCggg                                                                    <  1:4419519/151‑1 (MQ=255)
                          gCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCggg                                                                    <  1:3359100/151‑1 (MQ=255)
                           ccGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGt                                                                   <  1:4945289/151‑1 (MQ=255)
                            cggtGTAGGAGCGGTTTAACCCGCGAGAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGtt                                                                  <  2:5928878/148‑1 (MQ=14)
                               tGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCAc                                                               >  1:2610921/1‑151 (MQ=255)
                               tGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGCTTCAc                                                               >  2:6827571/1‑151 (MQ=255)
                               gggcggAGCGGGGGCACCCGCGAAAGGGGGGGGGCTGCGGGCCAAATTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCGCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCAc                                                               <  2:6430653/147‑1 (MQ=14)
                                gTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCAcc                                                              <  1:4568455/151‑1 (MQ=255)
                                  aGGGGCGGGTTCACCCGCGAAAGGGACGGGCCTGCTGGCCAATTTCCACCCGGCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcg                                                            <  2:1050510/151‑1 (MQ=14)
                                   ggggCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                                           <  2:2241318/148‑1 (MQ=14)
                                     aGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGa                                                         >  1:5753974/1‑151 (MQ=255)
                                      gCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGaa                                                        <  1:2960605/151‑1 (MQ=255)
                                       cGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGaaa                                                       <  2:2048370/151‑1 (MQ=255)
                                       cGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGCTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGCCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGaaa                                                       >  2:274867/1‑151 (MQ=14)
                                        gggTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGaaaa                                                      <  2:1241559/151‑1 (MQ=255)
                                         ggTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAg                                                     >  1:3523222/1‑151 (MQ=255)
                                         ggTACACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAg                                                     <  1:1375659/151‑1 (MQ=255)
                                           ttCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGc                                                   >  1:4750248/1‑151 (MQ=255)
                                           ttCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGCGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGc                                                   >  2:6679651/1‑151 (MQ=255)
                                             cACCAGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCg                                                 <  1:1165776/151‑1 (MQ=255)
                                              gcccGCGAAAGGGCCGGGCCTGCGGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCgg                                                <  2:4844428/150‑1 (MQ=255)
                                              aCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCgg                                                >  1:673431/1‑151 (MQ=255)
                                                ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTa                                              >  1:108191/1‑151 (MQ=255)
                                                 cgcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTAc                                             >  1:1355765/1‑151 (MQ=255)
                                                 cgcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTAc                                             <  2:3040818/151‑1 (MQ=255)
                                                   cgAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGca                                                                                                            >  2:6309836/1‑86 (MQ=33)
                                                   cgAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAg                                           <  1:109657/151‑1 (MQ=255)
                                                    gAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGCAAAGGCCGGTACAgg                                          >  2:1546063/1‑151 (MQ=255)
                                                    gAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAgg                                          <  1:3703072/151‑1 (MQ=255)
                                                      aaGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCg                                        >  2:3474977/1‑151 (MQ=255)
                                                       aGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGa                                       <  2:1121023/151‑1 (MQ=255)
                                                         ggccgggccgTCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGaaa                                     <  1:5219630/151‑1 (MQ=17)
                                                         ggccgggccTGCTGGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGaaa                                     <  1:4227453/151‑1 (MQ=255)
                                                         ggccgggccTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGaaa                                     >  2:797001/1‑151 (MQ=255)
                                                           ccgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAAt                                   >  2:2122055/1‑151 (MQ=255)
                                                            cgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATc                                  >  2:6027491/1‑151 (MQ=255)
                                                            cgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATc                                  >  2:5394309/1‑151 (MQ=255)
                                                            cgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATc                                  >  1:4132339/1‑151 (MQ=255)
                                                            cgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGGAAACCCGCTCCTGCACACGCAGCGGTGGTTGGTCACAGGCCGTCGCAGGAGCGGGTTCACGCGCGAAAAGGCCGGTACAGGCGgaaaga                                  >  1:2237234/1‑149 (MQ=17)
                                                                ccTGGTGGCCAATCTCCCCCCGCCGTTAACCGTTGCCGCCCCTACCGTAGTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGcc                              <  2:1984140/151‑1 (MQ=17)
                                                                 ctgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccc                             >  2:5371323/1‑151 (MQ=255)
                                                                 ctgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGCACAGGCGAACATCGGccc                             >  2:3934263/1‑151 (MQ=17)
                                                                 ctgctgGCCAATCTCCACCCGCCGTCCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGTATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccc                             <  2:2903844/151‑1 (MQ=255)
                                                                   gctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTGGGTCACAGGCCCTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGTGAAAATCGGccccc                           >  1:3842669/1‑151 (MQ=17)
                                                                    ctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCt                          >  1:1688198/1‑151 (MQ=255)
                                                                    ctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGGTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCt                          >  1:2169393/1‑151 (MQ=255)
                                                                     tgGCCAATTTCCACCCGTCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTg                         <  2:2075207/151‑1 (MQ=255)
                                                                     tgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTg                         <  1:2680190/151‑1 (MQ=255)
                                                                      ggCCAATTTCCACCCGCCGTCCACCGTTGCCGCCGCTAACGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGc                        <  2:1685591/151‑1 (MQ=255)
                                                                       gCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCt                       >  2:4135145/1‑151 (MQ=255)
                                                                       gCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCt                       >  1:3414305/1‑151 (MQ=255)
                                                                        ccAATCTCCGCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCtt                      <  1:4440450/151‑1 (MQ=255)
                                                                         cAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCttt                     <  1:5903861/151‑1 (MQ=255)
                                                                         cAATTTCCACCCGCCGTTCCCCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCttt                     <  2:458106/151‑1 (MQ=255)
                                                                         cAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCGGCTCCTGCACACGCGTCGGTGGGTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAATCCCGGGACAGGCGAAAATCGGCCCCCTGCgtt                     >  2:2670291/1‑151 (MQ=17)
                                                                         cAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCttt                     <  2:7118188/151‑1 (MQ=255)
                                                                          aaTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCtttt                    >  2:2446869/1‑151 (MQ=255)
                                                                           atttcCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTc                   >  1:5969721/4‑151 (MQ=255)
                                                                              tcCACCCGCCGTCCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTcccc                <  1:4016987/151‑1 (MQ=255)
                                                                                cACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCg              <  2:469093/151‑1 (MQ=255)
                                                                                 aCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCgg             >  1:5479926/1‑151 (MQ=255)
                                                                                  cccGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGa            >  1:3439063/1‑151 (MQ=255)
                                                                                     gccgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGCttc         >  2:1784159/3‑150 (MQ=255)
                                                                                     gccgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGCTTg         >  1:3266323/3‑151 (MQ=255)
                                                                                     gccgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGCTTg         >  2:5656195/3‑151 (MQ=255)
                                                                                     gccgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTg         >  2:1514131/3‑151 (MQ=255)
                                                                                       cgtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGc       >  2:4181472/1‑151 (MQ=255)
                                                                                        gtgcACCGTTCCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGcc      <  1:3071521/148‑1 (MQ=255)
                                                                                        gtTCGCCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGcc      >  2:5972313/1‑151 (MQ=255)
                                                                                        gtTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGCTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCAGCCCCCTGCTTTCCCCCCGGATTGGcc      >  2:2545027/1‑151 (MQ=255)
                                                                                         ttCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGCCa     >  1:761790/1‑151 (MQ=255)
                                                                                         ttCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGGACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGTTTGGCCa     >  1:5313188/1‑151 (MQ=255)
                                                                                            aCCGTTGCCGCCTCTCCAGTATTCGCGGGTAAGCCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGCCAttt  <  2:3219066/151‑1 (MQ=255)
                                                                                            aCCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGATTGGCCAttt  <  2:3514366/151‑1 (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: