Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 3,566,025 T→G intergenic (‑154/+160) RPPX_15935 ← / ← RPPX_15940 molecular chaperone/recombinase RdgC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099743,566,0250TG91.9% 838.8 / 48.2 273intergenic (‑154/+160)RPPX_15935/RPPX_15940molecular chaperone/recombinase RdgC
Reads supporting (aligned to +/- strand):  ref base T (15/7);  new base G (113/137);  total (128/145)
Fisher's exact test for biased strand distribution p-value = 4.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 273 total aligned reads displayed.

GTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACACACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAAGCCAGCACCGC  >  CP009974/3565878‑3566172
                                                                                                                                                   |                                                                                                                                                   
gtCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTgc                                                                                                                                                  <  1:1108068/151‑1 (MQ=255)
 tctcCATCGTTTGAATGAGGGGGCCGTCCCCGGAGTGGTCACGCGGGGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTgcc                                                                                                                                                 <  2:5266481/151‑1 (MQ=255)
    ccATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTgccggc                                                                                                                                              <  1:1650095/151‑3 (MQ=255)
      aTCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTgccggcgc                                                                                                                                            >  2:474016/1‑151 (MQ=255)
      aTCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTgccggcgc                                                                                                                                            <  2:1707050/151‑1 (MQ=255)
       tCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTgccggcgcc                                                                                                                                           <  2:2355528/151‑1 (MQ=255)
          tttGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCAcc                                                                                                                                        <  2:410761/151‑1 (MQ=255)
          tttGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCAcc                                                                                                                                        >  1:5312767/1‑151 (MQ=255)
                 gatgGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTcg                                                                                                                                 >  2:1431206/1‑151 (MQ=255)
                   tgGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTcgcg                                                                                                                               >  1:1025865/1‑151 (MQ=255)
                     gggCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTCTTcgcgcg                                                                                                                             >  2:4143409/1‑149 (MQ=255)
                       gCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTa                                                                                                                           >  2:2498428/1‑151 (MQ=255)
                        ccGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTaa                                                                                                                          >  2:1155379/1‑151 (MQ=255)
                         cGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTaaa                                                                                                                         >  1:4552369/1‑151 (MQ=255)
                         cGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTaaa                                                                                                                         <  1:3324753/151‑1 (MQ=255)
                           tCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAAcc                                                                                                                       >  1:2786575/1‑151 (MQ=255)
                                  aGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCt                                                                                                                <  1:499137/151‑1 (MQ=255)
                                  aGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCt                                                                                                                >  2:2866565/1‑151 (MQ=255)
                                     ggTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGt                                                                                                                            >  1:5017059/1‑136 (MQ=255)
                                     ggTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGt                                                                                                                            <  2:5017059/136‑1 (MQ=255)
                                      gTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTacac                                                                                                            >  2:2937528/1‑151 (MQ=255)
                                       tCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTAcaca                                                                                                           <  1:1378487/151‑1 (MQ=255)
                                        cACGCGGTGCGTCGGGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTGTTCGCGCGTAAACCCGCACCTACACAg                                                                                                          <  2:4051448/151‑1 (MQ=37)
                                          cgcgGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGac                                                                                                        >  2:4448554/1‑150 (MQ=255)
                                          cgcgGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGAt                                                                                                        <  1:1409400/151‑1 (MQ=255)
                                           gcgGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGAtt                                                                                                       <  1:3974378/151‑1 (MQ=255)
                                            cgGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATtg                                                                                                      >  1:4733421/1‑151 (MQ=255)
                                            cgGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGGAAACCCGCTCCTACACAGATtg                                                                                                      >  2:2175861/1‑151 (MQ=255)
                                             ggTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCAGCGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGCGTACACCGGCACCTACACAGATtgg                                                                                                     <  2:3393394/151‑1 (MQ=25)
                                              gTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATtggt                                                                                                    <  1:1001789/151‑1 (MQ=255)
                                              gTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATtggt                                                                                                    <  2:1213314/151‑1 (MQ=255)
                                                    cgtGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTGAACCCGCGCCTACACAGATTGGCCGTGGa                                                                                              >  2:1387467/1‑151 (MQ=25)
                                                    cgtGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGa                                                                                              <  1:5172379/151‑1 (MQ=255)
                                                     gtGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGac                                                                                             <  2:1985278/151‑1 (MQ=255)
                                                      tGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGaca                                                                                            <  1:2937528/151‑1 (MQ=255)
                                                      tGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTGCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGaca                                                                                            <  2:4980198/151‑1 (MQ=255)
                                                       gAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGacat                                                                                           <  2:3234211/151‑2 (MQ=255)
                                                        aGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGacata                                                                                          <  1:4386195/151‑3 (MQ=255)
                                                        aGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGacata                                                                                          <  2:4641491/151‑3 (MQ=255)
                                                           tcgtcCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGacatacat                                                                                       <  1:3799412/151‑1 (MQ=255)
                                                             gtcCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCg                                                                                     <  1:423239/151‑1 (MQ=255)
                                                             gtcCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCg                                                                                     >  1:1475342/1‑151 (MQ=255)
                                                              tcCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGCAAACCCGCCCCTACACAGATTGGTGGTGGACATACATCgg                                                                                    >  2:5097220/1‑151 (MQ=21)
                                                               ccTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGa                                                                                   <  1:1709721/151‑1 (MQ=255)
                                                                cTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGa                                                                                              <  2:2064659/139‑1 (MQ=255)
                                                                cTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAg                                                                                  <  2:1262500/151‑1 (MQ=255)
                                                                cTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAg                                                                                  >  1:695673/1‑151 (MQ=255)
                                                                 tCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAgg                                                                                 >  2:2555976/1‑151 (MQ=255)
                                                                    caGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTgg                                                                              <  1:5189453/151‑1 (MQ=255)
                                                                      gTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGa                                                                            >  1:371190/1‑151 (MQ=255)
                                                                       tGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACCGATCGGTGGTGGACATACATCGGACGTGGGac                                                                           >  2:4570757/1‑151 (MQ=18)
                                                                         ggCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGacac                                                                         <  2:576926/151‑1 (MQ=255)
                                                                            aCCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCCGGTAAACCCGCTCCTACACAGATTGGTGCTGGACATACATCGTCGGTGGGACACAgg                                                                      >  1:2564964/1‑151 (MQ=18)
                                                                             ccATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGCCTGGTGGTGGACACACATCGGCGGTGGGACCCACGc                                                                     >  2:1033009/1‑151 (MQ=18)
                                                                               aTGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCcg                                                                   >  1:2053666/1‑151 (MQ=255)
                                                                                    gCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTa                                                              <  1:3444114/151‑1 (MQ=255)
                                                                                     ccTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAg                                                             <  2:2767369/151‑1 (MQ=255)
                                                                                        ggCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTAGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGCTTGGTGGTGGACGTACATCGGAGGTGGGACGCAGGCCGCCGTCGGAg                                                          >  2:3733064/1‑151 (MQ=17)
                                                                                            gCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCggg                                                      >  1:1913261/1‑151 (MQ=255)
                                                                                              ggCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGtt                                                    >  1:55284/1‑151 (MQ=255)
                                                                                              ggCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGGCGTAGGAGCGGGtt                                                    >  2:1055918/1‑151 (MQ=255)
                                                                                              ggCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGCCGTAGGAGCGGGtt                                                    >  1:4673130/1‑151 (MQ=255)
                                                                                                   tCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTAccc                                               >  1:2702468/1‑151 (MQ=255)
                                                                                                     ggggCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCcgc                                             <  2:1110283/150‑1 (MQ=255)
                                                                                                     cGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGGCGCAGGCGCGGGGTTACCcgc                                             >  2:2418995/1‑151 (MQ=17)
                                                                                                           aGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAagag                                       <  1:4783633/151‑1 (MQ=255)
                                                                                                              ccggcccTGCTGGCCAATCCCCACCTGCCGTTCACCATTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCGCATACACAGATTGGTGGTGGACACGTATCGAGGGTGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCg                                   >  1:2468368/6‑151 (MQ=11)
                                                                                                                ggcccTGCTGGCCAATCCCCACCTGCCGTTCACCATTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCGCATACACAGATTGGTGGTGGACACGTATCGAGGGTGGACACAGGCCGTCGTAGGAGCGGGGTTACCCGCGAAGAGGCCggc                                 >  1:3082924/4‑150 (MQ=11)
                                                                                                                 gcccTGCTGGCCAATCCCCACCTGCCGTGCACCATTGCCGCCGCCGCCGTATTCGCGGGTAAACCCGCGCATACACAGATTGGTGGTGGACACGTATCGAGGGTGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCggcc                                <  2:190600/149‑3 (MQ=14)
                                                                                                                  gcCTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTg                                <  2:4237883/151‑1 (MQ=255)
                                                                                                                  cccTGCTGGCCAATCCCCACCTGCCGTTCACCATTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCGCATACACAGATTGGTGGTGGACACGTATCGAGGGTGGACACCGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCgtccc                               >  2:1420103/2‑147 (MQ=11)
                                                                                                                   ccTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGCGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCCCGAAGAGCCCGGTGc                               >  2:2622956/1‑151 (MQ=17)
                                                                                                                   ccTGCTGGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCCCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGGCGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGTGc                               >  1:4298086/1‑151 (MQ=14)
                                                                                                                   ccTGCTGGCCAATCCCCACCTGCCGTTCACCATTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCGCATACACAGATTGGTGGTGGACACGTATCGAGGGTGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCggccct                              >  1:400340/1‑147 (MQ=14)
                                                                                                                   ccTGCTGGCCAATCCCCACCTGCCGTTCACCATTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCGCATACACAGATTGGTGGTGGACACGTATCGAGGGTGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCggccct                              >  1:227724/1‑147 (MQ=11)
                                                                                                                   ccTGCTGGCCAATCCCCACCTGCCGTTCACCACTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCGCATACACAGATTGGTGGTGGACACGTATCGAGGGTGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCggccct                              <  1:5089664/151‑5 (MQ=14)
                                                                                                                    ctgctgGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCt                              <  2:3024624/151‑1 (MQ=255)
                                                                                                                    ctgctgGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCt                              <  2:1551419/151‑1 (MQ=255)
                                                                                                                    ctgctgGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCt                              <  2:4538832/151‑1 (MQ=255)
                                                                                                                    ctgctgGCCAATCCCCACCTGCCGTTCACCATTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCGCATACACAGATTGGTGGTGGACACGTATCGAGGGTGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCggccctg                             <  2:1309115/151‑6 (MQ=14)
                                                                                                                     tgctgGCTAATTTCCACACGCCGTTCACCGGCGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTg                             <  1:1550055/151‑1 (MQ=17)
                                                                                                                     tgctgGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACACACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTg                             >  1:1803193/1‑151 (MQ=255)
                                                                                                                      gctgGCCAATGGCCACCCGCCGTTCACCGGTGCGGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTgg                            <  2:3666853/151‑1 (MQ=17)
                                                                                                                      gctgGCCAATCCCCACCTGCCGTTCACCATTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCGCATACACAGATTGGTGGTGGACACGTATCGAGGGTGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCggccctgct                           <  1:4425363/151‑8 (MQ=14)
                                                                                                                      gctgGCCAAGCCCCACCTGCCGTTCACCATTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCGCATACACAGATTGGTGGTGGACACGTATCGAGGGTGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCggccctgct                           <  1:3261600/151‑8 (MQ=14)
                                                                                                                       ctgGCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGc                           >  1:5129433/1‑151 (MQ=255)
                                                                                                                          gCCAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGa                        <  1:1979371/151‑1 (MQ=255)
                                                                                                                           ccAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGaa                       <  1:4154037/151‑1 (MQ=255)
                                                                                                                            cAATTTCCCCCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGaaa                      <  1:2121253/151‑1 (MQ=255)
                                                                                                                            cAATTTCCACCCGCCGTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGaaa                      <  1:640452/151‑1 (MQ=255)
                                                                                                                                 tCCACCCGCCGTGGACGGGTGCGGGCGCCACCGGATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGAt                 <  1:3640229/151‑1 (MQ=18)
                                                                                                                                       cgccgTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAAGc           <  2:1773073/151‑1 (MQ=255)
                                                                                                                                         ccgTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAAGCCa         <  1:804900/151‑1 (MQ=255)
                                                                                                                                          cgTTCACCGGTGCCGGCGCCAGCGGATTCGCGGGTAAACCCGCTGCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAAGCCAg        <  2:5347587/151‑1 (MQ=21)
                                                                                                                                          cgTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACCTACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAAGCCAg        >  2:789473/1‑151 (MQ=255)
                                                                                                                                          cgTTCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTcaagcaag        >  2:635599/1‑151 (MQ=255)
                                                                                                                                             tCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGACGAGGCCGGTGCTGGCAGCAACAGATTCAAGCCAGCAc     >  2:4807289/1‑151 (MQ=18)
                                                                                                                                             tCACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAAGCCAGCAc     <  1:2415408/151‑1 (MQ=255)
                                                                                                                                             acACCGGTGCCGGCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACATACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAAGCCAGCAc     <  2:2675266/150‑1 (MQ=255)
                                                                                                                                               aCCGGTGCCGGCGACACCGTATTCGCGGGTAAACCCGCTCCTACACAAGTTGGTGGTGGAAATGCATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTGTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAAGCCAGCACCg   <  1:4451250/151‑1 (MQ=21)
                                                                                                                                                   |                                                                                                                                                   
GTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACAGTGGGCACCATGGTGCCTGGCTGCGGCCTTCCGGGCAAGTCCGCGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACAGATTGGTGGTGGACACACATCGGAGGTGGGACACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTGCTGGCAGAAAACGATTCAAGCCAGCACCGC  >  CP009974/3565878‑3566172

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: