Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 426,157 G→T 92.9% E33* (GAG→TAG)  yajD → HNH nuclease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273426,1570GT92.9% 564.7 / 16.0 236E33* (GAG→TAG) yajDHNH nuclease family protein
Reads supporting (aligned to +/- strand):  ref base G (20/2);  new base T (105/109);  total (125/111)
Fisher's exact test for biased strand distribution p-value = 1.72e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 238 total aligned reads displayed.

GCATCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCA  >  CP009273/426011‑426302
                                                                                                                                                  |                                                                                                                                                 
gCATCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAg                                                                                                                                                 <  2:3539209/149‑1 (MQ=255)
   tCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTCTCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAgttt                                                                                                                                              >  1:2655665/1‑149 (MQ=255)
   tCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAgttt                                                                                                                                              <  1:2959825/149‑1 (MQ=255)
         gAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTAt                                                                                                                                        >  2:2428336/1‑149 (MQ=255)
             cTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCa                                                                                                                                    <  2:3187407/149‑1 (MQ=255)
                   gTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGc                                                                                                                              >  1:315962/1‑149 (MQ=255)
                   gTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAACTATCCGTGGGTCTGAGGTCGCTGTTCCCGCGAGTTTGTTTATTACAACCTGc                                                                                                                              >  1:2602773/1‑149 (MQ=255)
                    tAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCg                                                                                                                             <  1:3789991/149‑1 (MQ=255)
                     aaCCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTGCCGCTAGTTTGTTTTTTCCAACCAGCGt                                                                                                                            >  1:799595/1‑149 (MQ=255)
                     aaCCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCAGCGt                                                                                                                            >  1:1140775/1‑149 (MQ=255)
                       cccTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGa                                                                                                                          >  1:622764/1‑149 (MQ=255)
                       cccTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGa                                                                                                                          >  1:3367999/1‑149 (MQ=255)
                         cTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAAc                                                                                                                        <  2:1660346/149‑1 (MQ=255)
                         cTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAAc                                                                                                                        >  2:1957051/1‑149 (MQ=255)
                          tACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACt                                                                                                                       <  1:1238971/149‑1 (MQ=255)
                          tACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACt                                                                                                                       <  1:241414/149‑1 (MQ=255)
                          tACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACt                                                                                                                       >  1:325853/1‑149 (MQ=255)
                            cATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTa                                                                                                                     <  2:1789918/149‑1 (MQ=255)
                              tttCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGGTTATTCCAACCTGCGTGAACTTAcc                                                                                                                   >  2:3942412/1‑149 (MQ=255)
                                  aTCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTACCCCTAGTTTGTTTATTAAAAACTGCGTGGAATTACCCtta                                                                                                               >  1:1753604/1‑148 (MQ=255)
                                    ccAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTcac                                                                                                             <  2:4109428/149‑1 (MQ=255)
                                            cTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGGTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGa                                                                                                     >  1:3606990/1‑149 (MQ=255)
                                             ttttCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGGCGCTGTTCACGCTAGTTTGTTTATTCCAACCTGCGTGAACTTTCCGTTCACCACATTGAt                                                                                                    >  2:2477612/1‑149 (MQ=255)
                                                tCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCAc                                                                                                 >  2:3004318/1‑149 (MQ=255)
                                                tCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCAc                                                                                                 >  2:3615724/1‑149 (MQ=255)
                                                  aTGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGa                                                                                               <  2:3744199/149‑1 (MQ=255)
                                                   tGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGAc                                                                                              >  2:1985220/1‑149 (MQ=255)
                                                   tGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGAc                                                                                              >  2:1830414/1‑149 (MQ=255)
                                                     gCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTTATCACGACCa                                                                                            >  1:3693222/1‑149 (MQ=255)
                                                     gCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCa                                                                                            <  1:3326596/149‑1 (MQ=255)
                                                     gCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCa                                                                                            <  1:3725462/149‑1 (MQ=255)
                                                     gCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGCGTTTGTTTATTGCAACCTGCGTGAACTTACCGTTCACCCCAATGATCACGACCa                                                                                            >  1:1553693/1‑149 (MQ=255)
                                                         tcatcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATAcc                                                                                        >  2:3428549/1‑149 (MQ=255)
                                                            tcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAAt                                                                                     >  2:870161/1‑149 (MQ=255)
                                                            tcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAAt                                                                                     >  2:3840641/1‑149 (MQ=255)
                                                              ccAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATaa                                                                                   <  2:754675/149‑1 (MQ=255)
                                                                    aaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTACCGCTAGGTTGTTTATTCCAACCAGCGTGAACTTACCGTTCACCACATAGATAAAGACCATAACAATAAcactga                                                                             >  1:4088837/1‑144 (MQ=255)
                                                                     aCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCCCCACATTGATCACGACCATACCAATAACCCGGaa                                                                            >  2:3869226/1‑149 (MQ=255)
                                                                       tATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTACAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGa                                                                          >  1:3982855/1‑149 (MQ=255)
                                                                        aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGAt                                                                         <  2:2229531/149‑1 (MQ=255)
                                                                         tGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTGCCGCGNGTTTTTTTATTCCAACCTGCGTGAACTTACCGTTCACCACAATGATCACGACCCTACCAATAAACCCGAAGATg                                                                        >  1:2081776/1‑149 (MQ=255)
                                                                         tGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATg                                                                        <  2:3232012/149‑1 (MQ=255)
                                                                         tGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAAAAACCCGGAAGATg                                                                        >  1:3960920/1‑149 (MQ=255)
                                                                           cgcgGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGc                                                                      <  1:1784838/149‑1 (MQ=255)
                                                                                ttAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTaa                                                                 <  2:583103/149‑1 (MQ=255)
                                                                                ttAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTAACGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTaa                                                                 >  2:821728/1‑149 (MQ=255)
                                                                                ttAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGGTTATTCCAACCTGCGTGAACTTACCGTTCCCCACATTGATCACGGCCATACCACTAACCCGGCAGATGGCAGTac                                                                 >  2:2627191/1‑148 (MQ=255)
                                                                                  aGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACt                                                               >  2:2354692/1‑149 (MQ=255)
                                                                                   gAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTg                                                              <  2:894691/149‑1 (MQ=255)
                                                                                     aaGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTggg                                                            >  2:3033726/1‑149 (MQ=255)
                                                                                     aaGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTggg                                                            <  1:2815389/149‑1 (MQ=255)
                                                                                      aGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCCGTAACTGGGa                                                           >  2:2325633/1‑149 (MQ=255)
                                                                                      aGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGa                                                           >  2:3872080/1‑149 (MQ=255)
                                                                                      aGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGa                                                           >  2:2283702/1‑149 (MQ=255)
                                                                                      aGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGa                                                           >  2:2064114/1‑149 (MQ=255)
                                                                                      aGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGa                                                           >  2:1610828/1‑149 (MQ=255)
                                                                                       gCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGaa                                                          <  1:531592/149‑1 (MQ=255)
                                                                                         ggCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAAtt                                                        <  2:636153/149‑1 (MQ=255)
                                                                                          gCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAAAAACCCGGAAGATGGCAGTAACTGGGAAttg                                                       >  2:140873/1‑149 (MQ=255)
                                                                                             aTCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCGAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAAttgttg                                                    <  2:1762878/149‑1 (MQ=255)
                                                                                                  gAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGtct                                               <  1:3195493/149‑1 (MQ=255)
                                                                                                  gAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGtct                                               >  2:2463604/1‑149 (MQ=255)
                                                                                                   aaaaaGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTGCCGCTAGGTTGTTTATTTCAAACAGCGTGAACATACCGTTAACCACAATGATAACCACCAAACCAAAAACCCCGAAGATGGCAGTAACTGGGGATTGTTttgtgtg                                              >  1:700018/1‑146 (MQ=255)
                                                                                                    aaaaGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTctct                                             >  2:1193510/1‑149 (MQ=255)
                                                                                                      aaGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTAt                                           >  2:3341270/1‑149 (MQ=255)
                                                                                                      aaGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTAt                                           >  2:3343126/1‑149 (MQ=255)
                                                                                                       aGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACCTTGATCACGACCATACCAATAACCCCGAAGATGGCAGTAACTGGGAATTTTTGTGTCTCTAtt                                          >  2:2782216/1‑149 (MQ=255)
                                                                                                       aGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTAtt                                          >  2:1429207/1‑149 (MQ=255)
                                                                                                       aGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTACAACCTGCGTGACCTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGAGGGCAGTAACTGGGAATTGTTGTGTCTCTAtt                                          >  2:1104535/1‑149 (MQ=255)
                                                                                                       aGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGGTTGTTTATTCCAACCAGCGTGAACATACCGTTCACCACATTTATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTAtt                                          >  1:688096/1‑149 (MQ=255)
                                                                                                          aTTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGcc                                       >  2:737591/1‑149 (MQ=255)
                                                                                                              aaaaTCTCTCCGTGGGTCTGCGGTCGCTGGTCCCGCTAGTTTGATTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGa                                   <  1:2481656/149‑1 (MQ=255)
                                                                                                              aaaaTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGa                                   >  2:2546531/1‑149 (MQ=255)
                                                                                                               aaaTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGAt                                  <  1:2763651/149‑1 (MQ=255)
                                                                                                                aaTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATc                                 <  2:2766340/149‑1 (MQ=255)
                                                                                                                aaTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATc                                 <  2:2655665/149‑1 (MQ=255)
                                                                                                                 aTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCa                                <  1:1359498/149‑1 (MQ=255)
                                                                                                                 aTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCa                                >  2:2721451/1‑149 (MQ=255)
                                                                                                                    tATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGa                             <  1:245992/149‑1 (MQ=255)
                                                                                                                    tATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCAGCGTGAACTTACCGTTCACCACCTTGATCACGACCCTACCAATAACCCGGCAGATGGCAGTAACGGGGAAATGTTGTGTCTCTATTGCCACGATCATGa                             >  2:1053503/1‑149 (MQ=255)
                                                                                                                      tCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGc                           >  2:3086618/1‑149 (MQ=255)
                                                                                                                      tCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGc                           <  1:4108009/149‑1 (MQ=255)
                                                                                                                       ccGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCa                          <  1:2772816/149‑1 (MQ=255)
                                                                                                                       ccGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCa                          <  2:21866/149‑1 (MQ=255)
                                                                                                                       ccGTGGGTCTGCGGTCGCTGTTACCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTAACGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGACTTGTTGTGTCTCTATTTCCACGATCATGAGCa                          >  2:2083497/1‑149 (MQ=255)
                                                                                                                        cGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCAt                         <  1:1990393/149‑1 (MQ=255)
                                                                                                                        cGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCAt                         <  1:3049973/149‑1 (MQ=255)
                                                                                                                         gTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCAtt                        <  2:37529/149‑1 (MQ=255)
                                                                                                                            ggTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCAAACCAATAACCCGGAAGATGGCAGTAACTGGGCATTGTTGTGTCTCTATTGCCACGCTCATGAGCATTCGa                     >  2:3528987/1‑149 (MQ=255)
                                                                                                                              tCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGaaa                   <  1:1438778/149‑1 (MQ=255)
                                                                                                                              tCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGaaa                   <  1:3805773/149‑1 (MQ=255)
                                                                                                                                  cGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAAtata               <  2:2860012/149‑1 (MQ=255)
                                                                                                                                   ggTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATAc              >  2:1264801/1‑149 (MQ=255)
                                                                                                                                   ggTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATAc              <  1:870161/149‑1 (MQ=255)
                                                                                                                                   ggTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCCTTCACCACATTTATAACGACCATACCAATAAACCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATAc              >  1:106541/1‑149 (MQ=255)
                                                                                                                                    gTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATAcc             >  1:1921351/1‑149 (MQ=255)
                                                                                                                                      cGCTGGTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGa           <  1:940295/149‑1 (MQ=255)
                                                                                                                                        cTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAg         >  2:3879227/1‑149 (MQ=255)
                                                                                                                                            tCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGa     >  1:2852629/1‑149 (MQ=255)
                                                                                                                                              ccGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATc   >  2:2955209/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCATCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCA  >  CP009273/426011‑426302

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: ‑