Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,364,366 G→A 95.0% intergenic (+106/‑107) pspE → / → ycjM thiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,364,3660GA95.0% 548.1 / 16.3 178intergenic (+106/‑107)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base A (62/107);  total (67/111)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 178 total aligned reads displayed.

GCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCT  >  CP009273/1364223‑1364509
                                                                                                                                               |                                                                                                                                               
gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTGTGTCACACCTTATTTAAATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTACCTAGTaa                                                                                                                                             >  2:3136595/1‑148 (MQ=255)
 ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTaaa                                                                                                                                            <  1:1904731/148‑1 (MQ=255)
    gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTGTGTGAAACCGTATTTAAATAACCAAAGATTTTTTGGGTCAGCGCACAAAAAAGCATTTTGACTAGTaaaaaa                                                                                                                                         >  1:556375/1‑148 (MQ=255)
       aGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGa                                                                                                                                      <  1:1169299/148‑1 (MQ=255)
          cATCGCAATGCCGAAGGTCAACGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGAttt                                                                                                                                   <  2:2675347/148‑1 (MQ=255)
          cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGAttt                                                                                                                                   <  2:2431905/148‑1 (MQ=255)
          cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGAttt                                                                                                                                   <  1:378866/148‑1 (MQ=255)
           aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTc                                                                                                                                  <  1:1654028/148‑1 (MQ=255)
            tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTcc                                                                                                                                 <  2:628568/148‑1 (MQ=255)
             cGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCATTGACACCATTTTTAAATACCCATAATTTTTTTTCGTAAACTCAAAAAAAAGCTTTTTACCAACTAAAAAATTAtttaca                                                                                                                                >  1:2919221/1‑145 (MQ=255)
               cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGc                                                                                                                              <  1:2314362/148‑1 (MQ=255)
                aaTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCg                                                                                                                             <  1:419884/148‑1 (MQ=255)
                  tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTGTGACACCTTTTTTAAATACCCAAAGATTATTTGCGTCAGCTCACAAATAAGCTTTTTACCTAGTAAAAAAATATTTACTGCGTg                                                                                                                           >  2:1020524/1‑148 (MQ=255)
                  tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTg                                                                                                                           <  1:1991034/148‑1 (MQ=255)
                    ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGAc                                                                                                                         <  1:794840/148‑1 (MQ=255)
                    ccGAAGGTCAAAGGTAAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGAc                                                                                                                         <  1:3008763/148‑1 (MQ=255)
                        aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTaaa                                                                                                                     <  2:1640924/148‑1 (MQ=255)
                         ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTaaaa                                                                                                                    <  2:738871/148‑1 (MQ=255)
                              acggTAAAAAGAATACACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                                                                                               <  2:3358/146‑1 (MQ=255)
                              aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                                                                                               <  2:1656691/148‑1 (MQ=255)
                               aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCGTATTTACATCCCCCAAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCGAGTAAAAAATTATTTACTGCGTGACTAAAACCCTtg                                                                                                              >  2:108567/1‑148 (MQ=255)
                                  ttAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGc                                                                                                           <  1:1647200/148‑1 (MQ=255)
                                     aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCa                                                                                                        <  1:2057976/148‑1 (MQ=255)
                                      aaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCaa                                                                                                       <  1:811949/148‑1 (MQ=255)
                                      aaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCaa                                                                                                       <  1:811951/148‑1 (MQ=255)
                                       aGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                                                                                                      <  1:2646926/148‑1 (MQ=255)
                                       aGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                                                                                                      <  1:1978769/148‑1 (MQ=255)
                                         aaTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                                                                                                    <  1:1585035/148‑1 (MQ=255)
                                           ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc                                                                                                  <  2:1981084/148‑1 (MQ=255)
                                            tCACCATGAGCGGCGCTTATGCGCCGTATTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACa                                                                                                 <  2:52193/148‑1 (MQ=255)
                                               ccATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGtt                                                                                              <  2:2141152/148‑1 (MQ=255)
                                               ccATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGtt                                                                                              <  2:3107893/148‑1 (MQ=255)
                                                 aTGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTa                                                                                            <  1:697463/148‑1 (MQ=255)
                                                  tGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTAAATCCCCATAGATTATTTTCGTCAGCTCACAAAAACGCTTTTTCCCAAGTAAAAAATTTTTTCATGCGTGACTAAAAACCGTGTGCTAAATTGACAGTTTAt                                                                                           >  2:2641307/1‑148 (MQ=255)
                                                  tGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCCTATTTAAATACCCAAAGATTATTTTAGTCATCTCACAAAAAAGCTTTTTACCTAGTAAAAAAATATTTACATGGTGACAAAAACCCATGTGGTCAAATTACAGTTTAt                                                                                           >  2:1373914/1‑148 (MQ=255)
                                                     gcggcgCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtttt                                                                                        <  2:1278442/148‑1 (MQ=255)
                                                     gcggcgCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtttt                                                                                        <  1:2024023/148‑1 (MQ=255)
                                                      cggcgCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTc                                                                                       <  1:3357390/148‑1 (MQ=255)
                                                        gcgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTg                                                                                     <  2:2706605/148‑1 (MQ=255)
                                                         cgcTTATGCGCCGTTTTTTTTTCTGTCACACCATATTTACAACCCCAAAGATTATTTTCGTCAGCTAACAAAAAAGCTTTTTACCAATTAAAAAAATATTTACTTCGTGAAAAAAAAACATTTTCTCAAATTAAACTTTAttttattc                                                                                    >  1:2349612/1‑144 (MQ=255)
                                                           cTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCgg                                                                                  <  1:295450/148‑1 (MQ=255)
                                                            ttATGCGCCGTTTTTTTTTATTTGAAAACATATTTTAAAACCCAAATATTATTTTAGTAAGCTAAAAAATAAGCATTTTAACTAATTAAAAATTATTTACATCGTTAATAAAACACATGTTCTGAATTTACAGTTTAttttattagga                                                                                 >  2:3273709/1‑141 (MQ=255)
                                                             tATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGag                                                                                <  1:3372556/148‑1 (MQ=255)
                                                               tGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagta                                                                              <  1:2734676/148‑1 (MQ=255)
                                                                    cGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGtctc                                                                         <  1:74476/148‑1 (MQ=255)
                                                                    cGTTTTTTTTTCTGTCACACCTTATTTAAATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAAAGATTTCCTGCGTGACAAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagtagtata                                                                         >  2:376050/1‑145 (MQ=255)
                                                                       ttttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCg                                                                      <  2:1551931/148‑1 (MQ=255)
                                                                        tttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGt                                                                     <  1:1737615/148‑1 (MQ=255)
                                                                        tttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGt                                                                     <  2:1831029/148‑1 (MQ=255)
                                                                           ttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTc                                                                  <  2:853354/148‑1 (MQ=255)
                                                                                 gTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGa                                                            <  1:2610896/148‑1 (MQ=255)
                                                                                   cacaCCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGAcc                                                          >  1:2987568/1‑148 (MQ=255)
                                                                                    acacCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCg                                                         >  2:2133043/1‑148 (MQ=255)
                                                                                      acCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGttt                                                                                             >  1:2117144/1‑110 (MQ=255)
                                                                                       ccTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTa                                                      >  1:849210/1‑148 (MQ=255)
                                                                                         ttATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTAcc                                                    >  1:1049601/1‑148 (MQ=255)
                                                                                          tATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCa                                                   >  1:1710136/1‑148 (MQ=255)
                                                                                          tATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCGAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCa                                                   <  1:48724/148‑1 (MQ=255)
                                                                                            tttACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACg                                                 <  1:1605678/148‑1 (MQ=255)
                                                                                            tttACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACg                                                 <  1:1605713/148‑1 (MQ=255)
                                                                                            tttACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACg                                                 >  2:1998717/1‑148 (MQ=255)
                                                                                                  tCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAgg                                           >  1:675920/1‑148 (MQ=255)
                                                                                                    cccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCa                                         <  1:341796/148‑1 (MQ=255)
                                                                                                     ccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCaa                                        <  2:3252429/148‑1 (MQ=255)
                                                                                                     ccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCaa                                        <  1:1933319/148‑1 (MQ=255)
                                                                                                      cATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAAc                                       <  1:2647358/148‑1 (MQ=255)
                                                                                                      cATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAAc                                       <  1:2828731/148‑1 (MQ=255)
                                                                                                       aTAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACa                                      >  1:200692/1‑148 (MQ=255)
                                                                                                       aTAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACa                                      >  2:3033594/1‑148 (MQ=255)
                                                                                                        tAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACAt                                     >  1:262909/1‑148 (MQ=255)
                                                                                                         aGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATg                                    >  2:1712702/1‑148 (MQ=255)
                                                                                                         aGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATg                                    >  2:3453671/1‑148 (MQ=255)
                                                                                                           attattTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGaa                                  >  1:1351314/1‑148 (MQ=255)
                                                                                                             tattTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAAc                                >  1:355259/1‑148 (MQ=255)
                                                                                                             tattTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAAc                                >  1:2466939/1‑148 (MQ=255)
                                                                                                              attTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACa                               >  2:21364/1‑148 (MQ=255)
                                                                                                               tttGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAg                              >  1:1020310/1‑148 (MQ=255)
                                                                                                               tttGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAg                              >  2:2823352/1‑148 (MQ=255)
                                                                                                               tttGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAg                              <  1:1489735/148‑1 (MQ=255)
                                                                                                               tttGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAg                              >  2:2373255/1‑148 (MQ=255)
                                                                                                               tttGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAg                              >  1:2772986/1‑148 (MQ=255)
                                                                                                                   cGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaaaa                          <  2:2554375/148‑1 (MQ=255)
                                                                                                                   cGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaaaa                          <  2:667059/148‑1 (MQ=255)
                                                                                                                     tCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAAtt                        >  2:2590762/1‑148 (MQ=255)
                                                                                                                     tCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAAtt                        >  2:2397334/1‑148 (MQ=255)
                                                                                                                     tCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAAtt                        >  1:1302309/1‑148 (MQ=255)
                                                                                                                       aGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTAc                      <  1:1104692/148‑1 (MQ=255)
                                                                                                                         cTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACgg                    >  1:1964352/1‑148 (MQ=255)
                                                                                                                         cTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACgg                    >  1:156521/1‑148 (MQ=255)
                                                                                                                          tCACAAATACGCTTGTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGGTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGa                   <  2:2386139/148‑1 (MQ=255)
                                                                                                                           cacaAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGAt                  <  1:3136595/148‑1 (MQ=255)
                                                                                                                            acaAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGAtt                 >  1:1784224/1‑148 (MQ=255)
                                                                                                                              aaaTACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTAc               <  2:3232790/148‑1 (MQ=255)
                                                                                                                              aaaTACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTAc               <  1:3061599/148‑1 (MQ=255)
                                                                                                                              aaaTACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTAc               <  2:2165118/148‑1 (MQ=255)
                                                                                                                                 tACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTg            >  2:1221523/1‑148 (MQ=255)
                                                                                                                                     cTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACg        >  2:185990/1‑148 (MQ=255)
                                                                                                                                       ttttCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGaa      >  2:2065528/1‑148 (MQ=255)
                                                                                                                                        tttCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGaaa     >  2:2038244/1‑148 (MQ=255)
                                                                                                                                        tttCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGaaa     >  1:1763526/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCT  >  CP009273/1364223‑1364509

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