Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,385,964 T→G 14.0% Q49H (CAA→CAC ycjY ← S9 homolog non‑peptidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,385,9640TG14.0% 795.9 / 52.6 378Q49H (CAA→CACycjYS9 homolog non‑peptidase family protein
Reads supporting (aligned to +/- strand):  ref base T (91/226);  new base G (59/2);  total (150/228)
Fisher's exact test for biased strand distribution p-value = 4.03e-25
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 388 total aligned reads displayed.

TAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTC  >  CP009273/1385815‑1386112
                                                                                                                                                     |                                                                                                                                                    
tAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGtt                                                                                                                                                      >  2:3463719/1‑150 (MQ=255)
 aatcaatcACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGttt                                                                                                                                                     <  1:5321705/150‑1 (MQ=255)
  atcaatcaCGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGtgtg                                                                                                                                                    >  1:2682065/1‑150 (MQ=255)
  atcaatcaCGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTg                                                                                                                                                    <  1:4289698/150‑1 (MQ=255)
     aatcaCGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCCATAGTGACAAAACCCTGTTCCGCCAGTATTTTGGCATAGGTTCCGGCCGtgtgggc                                                                                                                                                 >  1:2401240/1‑147 (MQ=255)
       tcaCGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCtt                                                                                                                                               <  1:2137753/150‑1 (MQ=255)
        caCGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCttt                                                                                                                                              <  1:448382/150‑1 (MQ=255)
         aCGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTa                                                                                                                                             <  2:682287/150‑1 (MQ=255)
            gCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACg                                                                                                                                          <  1:2984796/150‑1 (MQ=255)
            gCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTGGTTCTTTaaag                                                                                                                                          >  2:2325151/1‑148 (MQ=255)
            gCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGGTGttttaag                                                                                                                                          >  1:2294101/1‑146 (MQ=255)
              aCTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTTTTCTTTAAAGcc                                                                                                                                        >  1:145550/1‑150 (MQ=255)
              aCTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGcc                                                                                                                                        <  2:1279625/150‑1 (MQ=255)
                tGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGcccc                                                                                                                                      <  2:63942/150‑1 (MQ=255)
                 gATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAAAGccccc                                                                                                                                     >  1:4625918/1‑150 (MQ=255)
                  aTGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGTTAATTAAAGCCCCCa                                                                                                                                    >  2:1869867/1‑150 (MQ=255)
                    gTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGGTTGTTCTTTAACGCCCCCAcc                                                                                                                                  >  2:2533015/1‑150 (MQ=255)
                    gTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCATTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTATTTAAAGCCCCCAcc                                                                                                                                  >  2:4033449/1‑150 (MQ=255)
                     tCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCt                                                                                                                                 <  1:4817792/150‑1 (MQ=255)
                      cTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTAGCCGCCACTTTCAACTTGATAAGATTCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGTTTCCGGCCGTTTTTTGTTTAAAACCCCCAccag                                                                                                                                >  1:3440427/1‑148 (MQ=255)
                         cAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGAGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCAGGGt                                                                                                                             >  2:513530/1‑150 (MQ=255)
                            tACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTgag                                                                                                                          <  1:2215810/150‑1 (MQ=255)
                            tACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTgag                                                                                                                          <  2:1617849/150‑1 (MQ=255)
                              cGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGTCATAGGTTCCGGCCGTTTTTTGTTTAAAGCCCCCAACAAGGTGatat                                                                                                                        >  1:1174964/1‑150 (MQ=255)
                               gAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATa                                                                                                                       >  2:752966/1‑150 (MQ=255)
                               gAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATa                                                                                                                       <  1:3278327/150‑1 (MQ=255)
                               gAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAAAGACCCCACCAGGGTGAGATa                                                                                                                       >  2:3689977/1‑150 (MQ=255)
                               gAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTATTTTAAGCCCCCACCTGGGTGTGATa                                                                                                                       >  1:197618/1‑150 (MQ=255)
                                   ttATGGGTTTTCTAACTGACGCTGCTCGCAGCCACTTTCACCTTGATAAGTTGCGTCATAAGCAATAGTGACAAATCACTTTTACTCGATTTTTTTGGCATTGGTTACGTCCGTTTGTTCTTTAACGCCCCCACCTGGGTGTGATAacac                                                                                                                   <  2:2031799/149‑1 (MQ=255)
                                   gTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTTTTATTTAAAGCCCCCACCAGGGGGGGATaaaaa                                                                                                                   >  1:3171594/1‑147 (MQ=255)
                                   gTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAacac                                                                                                                   <  1:3408744/150‑1 (MQ=255)
                                     aTGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAAAGCCCCCACCCGGGGGGGTTAACACaa                                                                                                                 >  1:620519/1‑150 (MQ=255)
                                       gggTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAAtt                                                                                                               <  2:4520507/150‑1 (MQ=255)
                                       gggTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAAtt                                                                                                               <  2:4346589/150‑1 (MQ=255)
                                        ggTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGTTGGTTAAAGCCCCCACCTGGGTGGGATTAAACAAttt                                                                                                              >  1:3827609/1‑149 (MQ=255)
                                          ttttCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGTTGTTTATCGCCCCCACCCGGGGGGTATAACACAAtatac                                                                                                            >  1:4511592/1‑146 (MQ=255)
                                               tAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTtata                                                                                                       <  2:4480971/150‑1 (MQ=255)
                                                aaCTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTATTTAACGCCCCCACCGGGGGGTGATAACACAAATGCCATATAc                                                                                                      >  1:1199755/1‑150 (MQ=255)
                                                 aCTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACt                                                                                                     <  2:4245484/150‑1 (MQ=255)
                                                    gACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGAc                                                                                                  <  2:2691797/150‑1 (MQ=255)
                                                    gACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGAc                                                                                                  <  1:2323965/150‑1 (MQ=255)
                                                      cgcgGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTGTTTAACGCCCCCACCGGGGTGAGATAACACAAATGCCGTATAATGACgg                                                                                                >  1:2601557/1‑150 (MQ=255)
                                                       gcgGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACAGACggg                                                                                               >  2:4873232/1‑150 (MQ=255)
                                                            tCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTtca                                                                                          <  1:3107775/150‑1 (MQ=255)
                                                              gccgccACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTtcatc                                                                                        <  1:1258273/150‑1 (MQ=255)
                                                              gccgccACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTtcatc                                                                                        <  1:1709357/150‑1 (MQ=255)
                                                              gccgccACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTtcatc                                                                                        <  2:1158218/150‑1 (MQ=255)
                                                              gccgccACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAAAGCCCCCACCAGGGTGAGATAACACAAATTCCTTATAATGACGGGTTtcata                                                                                        >  1:3059899/1‑149 (MQ=255)
                                                               ccgccACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTtcatca                                                                                       <  1:2026133/150‑1 (MQ=255)
                                                                    aCTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGTGATAACACAATTGCCTTATAATGACGGGTTTCATCAAAttt                                                                                  >  1:820973/1‑150 (MQ=255)
                                                                     cTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTTACGCCCCCACCTGGGTGTGATAACACAATTGCCTTATACTGTCTGGTTTAATCAAAtttt                                                                                 >  2:4187879/1‑150 (MQ=255)
                                                                     cTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAAtttt                                                                                 <  1:5173290/150‑1 (MQ=255)
                                                                         cACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAg                                                                             <  2:4513042/150‑1 (MQ=255)
                                                                         cACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGTTCGTTAACGCCCCCACCGGGGTGAGATAACACAATTGCCGTATACTGACGGGTTTCATCAAAATTTGGAg                                                                             >  2:1196250/1‑150 (MQ=255)
                                                                            cTTGATAAGATGCTTAATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGTTTCCGGCCGTTTTTTCTTTAACGCCCCCACCTGGTTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGa                                                                          <  2:209900/150‑1 (MQ=255)
                                                                               gATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAAt                                                                       <  2:4402785/150‑1 (MQ=255)
                                                                                aTAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGTTTCCGGCCGTTTGTTCTTTAAAGCCCCCACCAGGGGGAGTTAACACAATTGCCTTATAATGACGGGTTTCATCAAATTTTTGAGGGAAATa                                                                      >  1:291193/1‑150 (MQ=255)
                                                                                aTAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATa                                                                      <  1:3950091/150‑1 (MQ=255)
                                                                                   aGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAAGTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGAt                                                                   <  2:952024/150‑1 (MQ=255)
                                                                                    gATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATc                                                                  <  2:1544789/150‑1 (MQ=255)
                                                                                    gATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATc                                                                  <  1:2526585/150‑1 (MQ=255)
                                                                                    gATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATc                                                                  <  2:5207987/150‑1 (MQ=255)
                                                                                    gATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACGGACGGGTTTCATCAAATTTTGGAGGGAAATAGATc                                                                  <  2:378431/150‑1 (MQ=255)
                                                                                     aTGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTTTCGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGAAGGGTTTCATAAAATTTTGGAGGGAAATAGAtaa                                                                 >  2:2777105/1‑148 (MQ=255)
                                                                                       gCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACt                                                               <  2:4934183/150‑1 (MQ=255)
                                                                                         gTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGGTGTTTAACGCCCCCCCCAGGGGGTGATAAAACAAATTCCGTATACAGTAGGGGTTGATCAAATTTTTGTGGGGAATAGATAACAGc                                                             >  1:1045400/1‑150 (MQ=255)
                                                                                            aTAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGa                                                          <  2:1318130/150‑1 (MQ=255)
                                                                                              aaGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCGTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATAAAATTTTTGAGGGAAATAGATAACAGCAGACa                                                        >  2:90939/1‑150 (MQ=255)
                                                                                               aGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATAACTGCAGACaa                                                       >  2:3923579/1‑150 (MQ=255)
                                                                                                   aTAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAagag                                                   <  1:3526900/150‑1 (MQ=255)
                                                                                                     aGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGAt                                                 <  1:4316089/150‑1 (MQ=255)
                                                                                                      gTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATg                                                <  1:2792210/150‑1 (MQ=255)
                                                                                                       tGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATgg                                               <  2:1476437/150‑1 (MQ=255)
                                                                                                            aaTCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGa                                          <  2:595268/150‑1 (MQ=255)
                                                                                                            aaTCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAAAACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACAGCAGACAAAGAGATGGTTGGa                                          >  2:3066906/1‑150 (MQ=255)
                                                                                                              tCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGatt                                        <  1:1869867/150‑1 (MQ=255)
                                                                                                                ccTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGattat                                      <  2:82611/150‑1 (MQ=255)
                                                                                                                ccTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGTTCTTTTAAGCCCCCACCCGGGTGAGATAACACAAATTCCTTATAATGACGGGGTTAATCAAATTTTGGAGGGAAATAGATAACAGCAGAAAAACAGATGGTTGGattat                                      >  1:4228973/1‑150 (MQ=255)
                                                                                                                  ttttCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTg                                    <  1:1802575/150‑1 (MQ=255)
                                                                                                                  ttttCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCAGGGTGAGATAACACAATTGCCTTATAATGACGGGTTTCATCAAAATTTGGAGGGAAAAAGATAACAGCAGACAAAGAGATGGGTGGATTATTg                                    >  1:407754/1‑150 (MQ=255)
                                                                                                                    ttCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCt                                  <  2:3236449/150‑1 (MQ=255)
                                                                                                                    ttCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCt                                  <  2:520096/150‑1 (MQ=255)
                                                                                                                    ttCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCt                                  <  2:4109419/150‑1 (MQ=255)
                                                                                                                       cGCCAGTTTTTTGGCATAGGTTCCGGCCGTGTGTTCTTTAAAGCCCCCACCGGGGTGAGATAACACAATTGCCATATAATGACGGGTTTCATCAAATTTTTGAGGGAAAAAGATAACTGCCGACAAAGAGATGGTTGGTTTATTGCTGtt                               >  2:1607024/1‑150 (MQ=255)
                                                                                                                         ccAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTTCTGTTAg                             >  2:3528133/1‑150 (MQ=255)
                                                                                                                         ccAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCAGGGTGAGATAAAACAAAAGCCTTATACAGACGGGTTTAATCAAATTTTGGTGGGAAATAGATAACTGCAGACAAAGAGATGGTTGGATTATTTCTGTTAg                             >  2:1050987/1‑150 (MQ=255)
                                                                                                                              tttttGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAAGTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAg                        <  2:3246744/150‑1 (MQ=255)
                                                                                                                                tttGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCt                      <  2:5056572/150‑1 (MQ=255)
                                                                                                                                    gCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAAGTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACt                  <  2:2252890/150‑1 (MQ=255)
                                                                                                                                      aTAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACttt                <  1:1609206/150‑1 (MQ=255)
                                                                                                                                        aGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTttat              <  2:3641433/150‑1 (MQ=255)
                                                                                                                                        aGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTttat              <  1:3733956/150‑1 (MQ=255)
                                                                                                                                         ggTTCCGGCCGATTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTttatt             <  2:2378389/150‑1 (MQ=255)
                                                                                                                                          gTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTttatta            <  2:4695883/150‑1 (MQ=255)
                                                                                                                                            tCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTattatt          <  2:2167149/150‑1 (MQ=255)
                                                                                                                                              cGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtca        <  2:1698499/150‑1 (MQ=255)
                                                                                                                                                gCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATtcatc      >  2:1081511/1‑150 (MQ=255)
                                                                                                                                                  cGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTcatcat    <  1:1505183/150‑1 (MQ=255)
                                                                                                                                                    tttGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTc  >  1:1075359/1‑150 (MQ=255)
                                                                                                                                                    tttGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTTATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGAATTTATTATTCATCATTc  <  1:1657246/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAATCAATCACGGCACTGATGTCTTCAGTACGAATGTATGGGTTTTCTAACTGACGCGGCTCGCCGCCACTTTCACCTTGATAAGATGCGTCATAAGCAATAGTGACAAATCCCTTTTCCGCCAGTTTTTTGGCATAGGTTCCGGCCGTTTGTTCTTTAACGCCCCCACCTGGGTGAGATAACACAATTGCCTGATACTGACGGGTTTCATCAAATTTTGGAGGGAAATAGATCACTGCAGACAAAGAGATGGTTGGATTATTGCTGTTAGTGAAGCTGACTTTATTATTCATCATTC  >  CP009273/1385815‑1386112

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