Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,030,927 G→A intergenic (+69/‑715) RPPX_17905 → / → RPPX_17915 23S rRNA methyltransferase/sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,030,9270GA99.6% 877.1 / NA 249intergenic (+69/‑715)RPPX_17905/RPPX_1791523S rRNA methyltransferase/sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (100/148);  total (100/149)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 249 total aligned reads displayed.

CGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGC  >  CP009974/4030779‑4031057
                                                                                                                                                    |                                                                                                                                  
cGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATa                                                                                                                                  <  1:6434011/151‑1 (MQ=255)
cGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATa                                                                                                                                  <  1:261293/151‑1 (MQ=255)
 gATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAg                                                                                                                                 >  2:6027584/1‑151 (MQ=255)
   ttGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGtt                                                                                                                               <  1:4481618/151‑1 (MQ=255)
     gCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCa                                                                                                                             <  2:7165906/151‑1 (MQ=255)
      ccTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAg                                                                                                                            <  2:2502091/151‑1 (MQ=255)
      ccTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAg                                                                                                                            <  1:5615829/151‑1 (MQ=255)
      ccTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAg                                                                                                                            >  1:5571689/1‑151 (MQ=255)
       cTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAgg                                                                                                                           <  2:4352737/151‑1 (MQ=255)
        tGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGc                                                                                                                          <  1:4718504/151‑1 (MQ=255)
         gCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCg                                                                                                                         >  1:1009773/1‑151 (MQ=255)
         gCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCg                                                                                                                         <  2:3147927/151‑1 (MQ=255)
           aaGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGAc                                                                                                                       <  2:2334130/151‑1 (MQ=255)
            agcagcTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACg                                                                                                                      <  2:3740941/151‑1 (MQ=255)
              cagcTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCAGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCa                                                                                                                    <  1:4585590/151‑1 (MQ=255)
                 cTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGAt                                                                                                                 >  2:563960/1‑151 (MQ=255)
                 cTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGAt                                                                                                                 >  2:4959880/1‑151 (MQ=255)
                  tGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGCCGCAGAtt                                                                                                                >  1:6396682/1‑151 (MQ=255)
                  tGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGAtt                                                                                                                <  2:4798150/151‑1 (MQ=255)
                      ccACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGa                                                                                                            <  1:470872/151‑1 (MQ=255)
                        aCGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCCGATTGCGACt                                                                                                          >  1:656261/1‑151 (MQ=255)
                        aCGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACt                                                                                                          <  2:2553837/151‑1 (MQ=255)
                         cGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACtt                                                                                                         >  1:3992768/1‑151 (MQ=255)
                         cGACCGCGAATAATTCACCTGCCGCCTGCGGCGGCTCGACTTGAGGCCACGCTGAGGCCAGCCCGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGGTTCTAGGCGCGCGATAGTCCAGGCGGCGCAGATTGCGACtt                                                                                                         <  2:1323983/151‑1 (MQ=21)
                          gACCTTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTc                                                                                                        >  2:3449248/1‑151 (MQ=255)
                           aCCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGAAGATTGCGACTTCa                                                                                                       >  1:3352613/1‑151 (MQ=255)
                            ccGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCag                                                                                                      <  2:6462039/151‑1 (MQ=255)
                             cGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCagg                                                                                                     <  1:1357914/151‑1 (MQ=255)
                              gTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCagga                                                                                                    <  1:4288981/151‑1 (MQ=255)
                                gaagaaGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCaggagg                                                                                                  <  2:7014535/151‑1 (MQ=255)
                                 aagaagTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGTGACTTCAGGAGGc                                                                                                 >  1:1251575/1‑151 (MQ=255)
                                  agaagTAACCTGCCGCCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGcc                                                                                                <  1:6319292/151‑1 (MQ=255)
                                  agaagTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGcc                                                                                                <  1:3758435/151‑1 (MQ=255)
                                  agaagTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGcc                                                                                                <  1:2529757/151‑1 (MQ=255)
                                    gagTGACCTGCCGCCTGCGGCGGCTCGGCTTGAGGCCACGCTGGGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGa                                                                                              <  2:5861221/150‑1 (MQ=18)
                                    aagTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGa                                                                                              <  1:1628933/151‑1 (MQ=255)
                                       tAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGCCTTCAGGAGGCCGAGTg                                                                                           >  2:6865920/1‑151 (MQ=255)
                                       tAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCagg                                                                                                     <  2:1430472/141‑1 (MQ=255)
                                       tAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCagg                                                                                                     >  1:1430472/1‑141 (MQ=255)
                                        aaCCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGc                                                                                          <  1:906873/151‑1 (MQ=255)
                                         aCCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCa                                                                                         <  1:5713116/151‑1 (MQ=255)
                                          ccTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAg                                                                                        >  1:3391875/1‑151 (MQ=255)
                                            tGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAggg                                                                                      >  2:2299155/1‑151 (MQ=255)
                                             gCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGc                                                                                     >  2:279417/1‑151 (MQ=255)
                                               cACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCtt                                                                                   <  2:6712189/151‑1 (MQ=255)
                                               cACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCtt                                                                                   >  2:801042/1‑151 (MQ=255)
                                               cACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCtt                                                                                   >  1:5190959/1‑151 (MQ=255)
                                                  cTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTTCg                                                                                >  2:5110952/1‑151 (MQ=255)
                                                  cTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCg                                                                                <  1:4812322/151‑1 (MQ=255)
                                                  cTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCg                                                                                <  2:3079306/151‑1 (MQ=255)
                                                   tGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGa                                                                               >  2:4077403/1‑151 (MQ=255)
                                                    gcggcgGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGaa                                                                              >  1:3913542/1‑151 (MQ=255)
                                                    gcggcgGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGaa                                                                              >  2:1835260/1‑151 (MQ=255)
                                                    gcggcgGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGaa                                                                              <  1:1898974/151‑1 (MQ=255)
                                                    gcggcgGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGaa                                                                              <  1:760301/151‑1 (MQ=255)
                                                      ggcggcTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAgg                                                                            >  1:2788626/1‑151 (MQ=255)
                                                      ggcggcTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAgg                                                                            >  2:4782809/1‑151 (MQ=255)
                                                       gcggcTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAggg                                                                           <  1:2419470/151‑1 (MQ=255)
                                                        cggcTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGa                                                                          <  2:1481869/151‑1 (MQ=255)
                                                        cggcTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGa                                                                          <  1:6791977/151‑1 (MQ=255)
                                                         ggcTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAg                                                                         <  1:4385023/151‑1 (MQ=255)
                                                            tCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTg                                                                      >  1:311063/1‑151 (MQ=255)
                                                              gACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGAc                                                                    >  2:2373078/1‑151 (MQ=255)
                                                               aCTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACg                                                                   >  2:422376/1‑151 (MQ=255)
                                                                cTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACgg                                                                  <  1:3668895/151‑1 (MQ=255)
                                                                  tGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGc                                                                >  1:790700/1‑151 (MQ=255)
                                                                   gAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCa                                                               <  2:587707/151‑1 (MQ=255)
                                                                     aGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATg                                                             <  2:1705061/151‑1 (MQ=255)
                                                                          cGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGc                                                        >  1:2151335/1‑151 (MQ=255)
                                                                          cGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGc                                                        <  1:1851200/151‑1 (MQ=255)
                                                                          cGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCCTGCGAAGGGAGGTGACGGGCATGGATGc                                                        >  2:5508304/1‑151 (MQ=255)
                                                                             tGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCg                                                     <  2:1009775/151‑1 (MQ=255)
                                                                              gAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGt                                                    <  2:6614308/151‑1 (MQ=255)
                                                                              gAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGt                                                    <  1:689175/151‑1 (MQ=255)
                                                                               aGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTc                                                   <  1:3584596/151‑1 (MQ=255)
                                                                                  ccAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGTGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATCCATGCCCGTCacc                                                <  2:1251575/151‑3 (MQ=14)
                                                                                  ccAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAg                                                <  2:5759747/151‑1 (MQ=255)
                                                                                    aGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAgcg                                              >  1:3220638/1‑151 (MQ=255)
                                                                                     gCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAgcgc                                             >  2:6736893/1‑151 (MQ=255)
                                                                                       aCGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTg                                           <  1:4884561/151‑1 (MQ=255)
                                                                                       aCGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTg                                           <  1:4637587/151‑1 (MQ=255)
                                                                                       aCGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATCCATGCCCGTCacctgcct                                           >  1:4470289/1‑144 (MQ=14)
                                                                                        cgtcgtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTgg                                          <  1:3260293/151‑2 (MQ=255)
                                                                                          tcgtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTgggc                                        >  1:3166324/1‑148 (MQ=255)
                                                                                           cgtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTgggcc                                       >  1:6725597/1‑147 (MQ=21)
                                                                                            gtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCTCTG‑GGccc                                     >  1:1282355/1‑151 (MQ=255)
                                                                                              tctgcttctgcttctgcttctgcttctgcttctgctTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCa                                   >  1:4050422/1‑151 (MQ=255)
                                                                                                tgcttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAgg                                 <  2:1965994/151‑1 (MQ=16)
                                                                                                tgcttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAgg                                 >  2:1032245/1‑151 (MQ=16)
                                                                                                   ttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATg                              <  1:5335912/151‑1 (MQ=16)
                                                                                                    tctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATgg                             >  2:1386478/1‑151 (MQ=16)
                                                                                                     ctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATggg                            <  1:1993146/151‑1 (MQ=255)
                                                                                                      tgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATgggg                           >  2:7068625/1‑151 (MQ=19)
                                                                                                       gcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGc                          <  2:1509729/151‑1 (MQ=16)
                                                                                                          tctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGCCtg                       <  2:2066250/151‑2 (MQ=16)
                                                                                                          tctgcttctgcttccgcctctgctgTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGCCtg                       <  2:3791108/151‑2 (MQ=11)
                                                                                                           ctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGCCtgc                      <  2:5374913/151‑3 (MQ=12)
                                                                                                            ttcttcttcttttgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGCCtgca                     >  2:4569028/3‑148 (MQ=11)
                                                                                                              cttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGCCtgcagc                   <  1:1934493/151‑6 (MQ=9)
                                                                                                               ttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGCCtgcagcg                  >  2:3509149/1‑145 (MQ=255)
                                                                                                                                                    |                                                                                                                                  
CGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGC  >  CP009974/4030779‑4031057

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