Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 3,837,678 C→T intergenic (+229/‑322) 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,6780CT85.4% 150.8 / 21.0 89intergenic (+229/‑322)RPPX_17095/RPPX_17100proteophosphoglycan precursor/cupin
Reads supporting (aligned to +/- strand):  ref base C (9/4);  new base T (31/45);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 7.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

TTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTA  >  CP009974/3837529‑3837816
                                                                                                                                                     |                                                                                                                                          
ttCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCACTTTTCCCCCGCCGCTCACCGTTGCTGCCGCCAccc                                                                                                                                        >  2:2202478/1‑150 (MQ=2)
 tCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTAc                                                                                                                                          >  1:2443675/1‑151 (MQ=255)
 ccACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTAc                                                                                                                                          <  2:709736/150‑1 (MQ=255)
  cACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTAcc                                                                                                                                         <  1:992327/151‑1 (MQ=255)
  cACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCAGCTGGCCAATTTCCACCCGCCCTTCACCCTTGCCGCCGCTAcc                                                                                                                                         >  2:633630/1‑151 (MQ=255)
   aCCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGGTCACCCGCGCAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTAccc                                                                                                                                        >  1:217464/1‑150 (MQ=255)
    ccATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGt                                                                                                                                       <  1:70630/151‑1 (MQ=255)
    ccATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTCGCCAATTACCACCCGCCGTTCACCGTTGCCGCCGCTACCGt                                                                                                                                       >  1:802435/1‑151 (MQ=255)
    ccATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTGGCAGGAGCGGGGTCACCCGCGAGAGAGCCGGGCCAGCTGGCCAACTTCCCCCCGCCCTTCACCGTCGCGGCCGCCACTGt                                                                                                                                       >  2:2229877/1‑151 (MQ=14)
     cATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTCCCGtt                                                                                                                                      >  1:2421116/1‑150 (MQ=255)
     cATTGCCGCCGCTACCGTAGTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTa                                                                                                                                      <  1:1134483/151‑1 (MQ=255)
      aTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGCGGGCGGCACAGGCCGCTGTAGGCGCGGGGTCACCCGCGCAAGGGCCGGGCCTGCTGGCCAATTTCCAACCGCCGTTCCCCGTTGCCGCCGCTACCGTAt                                                                                                                                     >  2:4766797/1‑151 (MQ=11)
       ttGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTAtt                                                                                                                                    <  1:2013621/151‑1 (MQ=38)
       ttGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTAtt                                                                                                                                    <  1:2013048/151‑1 (MQ=38)
        gggcgccgcTACCGTAGGCGCGGGTAGCCCCGCTCCGGCACACGCAGCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTc                                                                                                                                   <  2:2407442/148‑1 (MQ=11)
         gccgccGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTcg                                                                                                                                  <  1:799036/151‑1 (MQ=38)
           cgccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTcgcg                                                                                                                                >  1:2949723/1‑151 (MQ=37)
            gccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCgg                                                                                                                               <  1:5104766/151‑1 (MQ=38)
            gccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCgg                                                                                                                               >  2:4294160/1‑151 (MQ=38)
              cgcgagcGTATTCGCGGGTAAACCCGCTCCTACACACGCAGCGAGGGCGGCAGAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGt                                                                                                                             <  1:5269376/145‑1 (MQ=11)
              cgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGt                                                                                                                             <  1:4069824/151‑1 (MQ=37)
              cgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGt                                                                                                                             <  1:573860/151‑1 (MQ=38)
              cgcTACCGTATGCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGt                                                                                                                             <  1:4071443/151‑1 (MQ=255)
               gctagggTATTCGCGGGTAAACCCGCTCCTACGCACGCGGCGAGGGCGGCACAGGCCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTa                                                                                                                            <  2:2284682/145‑1 (MQ=11)
                cTACGGGATTCGCTGGTAAACCCGCTGCTACACGCGCAGCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTaa                                                                                                                           <  2:1353523/151‑1 (MQ=11)
                cTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTaa                                                                                                                           <  1:3941279/151‑1 (MQ=37)
                 tACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTaaa                                                                                                                          >  1:2836770/1‑151 (MQ=21)
                 tACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCCGTGCTGCCCCTTTCTTATTCGCGGGTaag                                                                                                                          >  2:4340562/1‑150 (MQ=11)
                  agggTAGTCGCGGGTAGACCGGGTCCAACACCGGCGTCGAGGGCGGCACAGGCCGTCGTATGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAc                                                                                                                         <  2:3287090/148‑1 (MQ=11)
                  aCCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAc                                                                                                                         >  2:943127/1‑151 (MQ=25)
                   gcggaTTCGCGGGTAAACCCGCCGCTACACACGCATCGAGGGCGGCCCAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAcc                                                                                                                        <  2:397137/147‑1 (MQ=11)
                   ccGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAcc                                                                                                                        >  2:3492075/1‑151 (MQ=18)
                     gTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCg                                                                                                                      >  1:1120971/1‑151 (MQ=25)
                      tATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGc                                                                                                                     <  1:785044/151‑1 (MQ=25)
                      tATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGc                                                                                                                     <  1:786730/151‑1 (MQ=25)
                       aTTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGCAAACCCGcc                                                                                                                    >  2:2885058/1‑150 (MQ=25)
                        ttCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTc                                                                                                                   >  1:132710/1‑151 (MQ=21)
                         tCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTcc                                                                                                                  >  1:899208/1‑151 (MQ=21)
                          tgcgGGTGAACCGGGTCCTACACATGCAGCGAGGGCGGCACAGGCCGTTGTAGGGGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCt                                                                                                                 <  2:2703079/150‑1 (MQ=11)
                          cgcgGGTAAGACCGCTCCTACACCCGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCt                                                                                                                 <  1:1055907/151‑1 (MQ=255)
                          cgcgGGTAAACCCGCTCCTACCCACGCAGCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCt                                                                                                                 <  1:2410448/151‑1 (MQ=255)
                          cgcgGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCt                                                                                                                 >  1:3332573/1‑151 (MQ=25)
                           gcgGGTAAACCCGCTCCTGCACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCtg                                                                                                                <  2:1552571/151‑2 (MQ=19)
                            cgGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTACACCCGCTCCtgc                                                                                                               >  2:2985084/1‑149 (MQ=12)
                            cgGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCtgc                                                                                                               >  1:493947/1‑149 (MQ=18)
                            cgGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCtgc                                                                                                               >  1:5261179/1‑149 (MQ=18)
                             gggTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTgca                                                                                                              <  1:1680035/151‑1 (MQ=18)
                              ggTAACCCGGCTCCGACACACGCATCGAGGGTGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTgcac                                                                                                             <  2:807970/151‑1 (MQ=11)
                              ggTAAACCCGCTCCTACACACGCATCGGGGGCGGCACAGGCCGTTTGAGGAGCGGGTTCACCCGCGAAAGGGCCGTGCCTGCTGGCCACTTTCTACCCGCAGTCCACCGTTGCGCCGGCTACTCTATTCGCTGGTACAACGGCTCCTtcac                                                                                                             >  2:3456076/1‑151 (MQ=11)
                              gctagaCCCGCCCCTACGCAGGCATCGAGGGCTGCAGAGGCCGTTGTAGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTgcac                                                                                                             <  2:4290510/146‑1 (MQ=11)
                               gTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTcccgcacc                                                                                                            >  1:1958521/1‑145 (MQ=16)
                                tAACCCCGCTCCTACACACGCCTCGTGGGCGGCGCTGCCCCTTGTACGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTgcacac                                                                                                           <  2:4356744/151‑1 (MQ=11)
                                 aaaCCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACg                                                                                                          >  2:3752476/1‑151 (MQ=18)
                                 aaaCCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACg                                                                                                          >  2:731777/1‑151 (MQ=18)
                                     ccGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCCCCCGCACACGCCTc                                                                                                      >  2:2268122/1‑151 (MQ=255)
                                       gCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCgg                                                                                                    >  2:5288804/1‑150 (MQ=2)
                                        cTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCggt                                                                                                   <  1:5078801/151‑3 (MQ=9)
                                        cTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCggt                                                                                                   >  1:162919/1‑149 (MQ=9)
                                                                                        cagccGCGAAAGGGGCGGGCCTGCTGGTCAATTTTCACCGGCCGATCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGGGTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCct                                                <  1:4181969/148‑1 (MQ=2)
                                                                                         cccagcgATAGGGCCGGGTCTGCTGGCCAATTCTCACCTGCCGTTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACTCTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctg                                               <  2:1787585/147‑1 (MQ=9)
                                                                                         aCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctg                                               <  1:4284316/151‑1 (MQ=2)
                                                                                         aCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctg                                               <  1:5286651/151‑1 (MQ=2)
                                                                                           ccGCGAAGGGGCCGGTCCTGCTGGCCAATTTTCACCAGCCGTTCTCGGTTGCTGCCGCTACCGTAGTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                                             <  2:2345951/151‑1 (MQ=2)
                                                                                           ccGCGAAAGGGCCGGGCGTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCGGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                                             <  1:3502551/151‑1 (MQ=2)
                                                                                           ccGCGAAAGGGCCGGGCCTGCTGGCGAATTTTCACCAGCCGTTCAACGTTGCTGCCGCTACCGTATTCGCGGGTAACCCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                                             <  1:4434930/151‑1 (MQ=2)
                                                                                           ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTCCCGTATTCGCGGGTAACCCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                                             <  2:31401/151‑1 (MQ=2)
                                                                                           ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                                             >  2:235082/1‑151 (MQ=11)
                                                                                           ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                                             <  1:170097/151‑1 (MQ=2)
                                                                                            cgcgAAAGGGCCGGGCCTGCTGGCCAATTTTCACCGGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgctg                                            <  2:150742/151‑1 (MQ=11)
                                                                                            cgcgAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgctg                                            <  2:3227558/151‑1 (MQ=2)
                                                                                             ccgAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTgg                                           <  2:2073831/150‑1 (MQ=2)
                                                                                              cgAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGACAGGGCCGGGCCTGCTGGc                                          >  2:4041131/1‑151 (MQ=2)
                                                                                               gAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCCACCGAATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGcc                                         <  2:5069675/151‑1 (MQ=9)
                                                                                                aaaGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCa                                        <  1:3455756/151‑1 (MQ=2)
                                                                                                 aaGGGCCGGGCCTGCTGGCCAGTTTTCACCAGCCGTTCACCGTTGCTGCCGCTCCCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCaa                                       <  2:1451941/151‑1 (MQ=2)
                                                                                                 aaGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCaa                                       <  2:3097971/151‑1 (MQ=2)
                                                                                                  aGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAt                                      >  1:5408570/1‑151 (MQ=2)
                                                                                                  aGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAt                                      <  1:618014/151‑1 (MQ=2)
                                                                                                  aGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAt                                      <  1:39988/151‑1 (MQ=2)
                                                                                                       cgggccTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAttttca                                 <  1:4166300/151‑6 (MQ=2)
                                                                                                       cgggccTGCTGGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAttttcc                                 >  1:3461420/1‑146 (MQ=2)
                                                                                                         ggccTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTACACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAcc                               <  1:615587/151‑1 (MQ=2)
                                                                                                           ccTGCTGGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccag                             >  1:3869106/1‑149 (MQ=2)
                                                                                                            ctgctgGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagc                            >  2:613653/1‑148 (MQ=2)
                                                                                                              gctgGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagccg                          >  1:2480953/1‑146 (MQ=2)
                                                                                                                tgGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagccgcc                        >  2:2099875/1‑144 (MQ=2)
                                                                                                                tgGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagccgct                        >  2:862381/1‑144 (MQ=2)
                                                                                                                tgGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagccgcc                        >  2:567042/1‑144 (MQ=2)
                                                                                                                                      cACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCCGGCCCATCTTCACCACCCGCTCACCCCTGCCCCCGCCACCGTa  >  2:1106872/1‑151 (MQ=2)
                                                                                                                                                     |                                                                                                                                          
TTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTA  >  CP009974/3837529‑3837816

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