Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,364,366 G→A 57.4% 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,3660GA57.4% 150.8 / inf 336intergenic (+106/‑107)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base G (55/86);  new base A (67/124);  total (125/211)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 345 total aligned reads displayed.

CCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGT  >  CP009273/1364217‑1364514
                                                                                                                                                       |                                                                                                                                                    
ccGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                                                                                                                                                      <  1:3023889/150‑1 (MQ=255)
 cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                                                                                                                                                     <  1:4786329/150‑1 (MQ=255)
 cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                                                                                                                                                     <  1:592632/150‑1 (MQ=255)
   gTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTATTTCACACCATATTTAAAAACCCATAGATTATTTTCGGAAGCGCACAAATAAGCATTTTACCAGGTa                                                                                                                                                   >  1:3345205/1‑150 (MQ=255)
       ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCAAACCTTATTTAAATCCCCATAGATTATTTGCGGCAGCTCACAAAAACGCTTTTTACCGGGTaaaaa                                                                                                                                               >  2:1106413/1‑150 (MQ=255)
         tGAAAGACATCGCAATGCCGATGGTCAACGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATg                                                                                                                                            <  2:3288115/150‑1 (MQ=255)
         tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt                                                                                                                                             <  2:4406419/150‑1 (MQ=255)
         tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTGTTTTTTTATGTGACACCATATTTACAAACCCAAAGATTATTTGCGTCAGCGCACCAATAAGCATTTTGCCGAGTaaaaaaa                                                                                                                                             >  1:1193508/1‑149 (MQ=255)
           aaaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTGTGTCACACCATATTTACAACCCCAAAGATTATTTTCGTCAGCTCACAAATAAGCATTTTACCAAGTAAAAAATGa                                                                                                                                           >  1:293325/1‑150 (MQ=255)
           aaaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGa                                                                                                                                           <  1:1063887/150‑1 (MQ=255)
              gACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGAttt                                                                                                                                        <  2:3382749/150‑1 (MQ=255)
               aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTATGTGACACCATATTTAAATACCCAGCGATTATTTGCGTCACCGCACAAATACCCATTTTCCCGCGTAAAAAATGAtttg                                                                                                                                       >  2:1956163/1‑149 (MQ=255)
                  tCGCAATGCCGAATGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATTCGCTTTTTCCCTGGTAAAAAATGATTTCCTg                                                                                                                                    <  1:2877607/150‑1 (MQ=255)
                  tCGAAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTg                                                                                                                                    <  1:345167/150‑1 (MQ=255)
                    gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCg                                                                                                                                  <  1:5251192/150‑1 (MQ=255)
                     cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGt                                                                                                                                 <  1:3232234/150‑1 (MQ=255)
                            gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCATATTTACATACCCATAGATTATTTGCGTCAGCTCACAAATACCCATTTTACCTAGTAAAAAATGATTTCCAGCGTGACTaaa                                                                                                                          >  2:1223072/1‑150 (MQ=255)
                             aaGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCGTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTTAAAAATGTTTTCCAGCGTGAATaaaa                                                                                                                         >  2:845429/1‑150 (MQ=255)
                              aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTATGTGACACCGTTTTTAAAACCCCCAAGAATATTTTAGTCACCTCACAAATAAGCATTTTACCGAGTAAAAAATTTTTTACAGCGGGACTaaaaa                                                                                                                        >  1:4392567/1‑149 (MQ=255)
                               ggTCAAAGGTTAAAAGAATTCACCATGCGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAAc                                                                                                                        <  1:1538240/150‑1 (MQ=255)
                                 tCAAAGGTTAACAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCAGAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCt                                                                                                                     <  2:3604784/150‑1 (MQ=255)
                                   aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTGTGAAACCATATTTAAAAACCCCAAGATTATTTTTGTAAGATAACAAAAAACCATTTTTACGAGTAAAAAATGTTTTACTGCGTGAATaaaaaacaag                                                                                                                   >  1:5175706/1‑144 (MQ=255)
                                     aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTTTCACACCCAATTTAAAAACCCCAAGATTAATTGAGTCAGCGAACAAAAAAACATTTTACCATGTTAAAAAATATTTACATCGTGAAAAAAAcacagttg                                                                                                                 >  1:1673961/1‑143 (MQ=255)
                                     aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTATTTAAAACCATAATTAAAAACCCCAAGAATATTTTCGTCAGCGCACAAAAAAGCATTTTACCGAGTAAAAAAAAATTTAAAGAGTGACAAAAAaccaagtg                                                                                                                 >  1:4334431/1‑145 (MQ=255)
                                       gTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCt                                                                                                               <  2:320411/150‑1 (MQ=255)
                                       gTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCt                                                                                                               <  2:2988621/150‑1 (MQ=255)
                                         tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTGATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAATATGATTTCCTGCGTGACTAAAACCCTTGTGCTCa                                                                                                             <  1:3110086/150‑1 (MQ=255)
                                          aaaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCaa                                                                                                            <  2:1505093/150‑1 (MQ=255)
                                           aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                                                                                                           <  1:3551490/150‑1 (MQ=255)
                                            aaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAtt                                                                                                          <  2:424750/150‑1 (MQ=255)
                                            aaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAtt                                                                                                          <  2:7267/150‑1 (MQ=255)
                                               aaTTCACCATGAGCGGCGCTTATGCGCCGATTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc                                                                                                       <  1:5049447/150‑1 (MQ=255)
                                                aTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACa                                                                                                      <  1:3858033/150‑1 (MQ=255)
                                                 ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAg                                                                                                     <  1:853715/150‑1 (MQ=255)
                                                      cATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAt                                                                                                <  2:1004280/150‑1 (MQ=255)
                                                       aTGAGCGGCGCTTATGCGCCGTTTTTTTTTATTTGACACCATAATTAAAAACCCAAAGAATAATTTCGTCAGCTCACAAAAACGCATTTTACCAAGTAAAAAAATATTTACAGCGTGTCAAAAAACCAATTGCTAAAATTACAGTTTAtt                                                                                               >  2:1435928/1‑150 (MQ=255)
                                                       aTGAGCGGCGCTTATGCGCCGTTTTTTTTTATGTCACACCCTATTTAAAAACCCCAAGAATATTTTCGTCAGCTCACAAAAAAGCATTTTACCCGGTAAAAAATTAATTACAGAGGGTCTAAAAAACAAGTGGTCAAATGAGAGTTTAtt                                                                                               >  2:5136561/1‑150 (MQ=255)
                                                        tGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAttt                                                                                              <  2:1151952/150‑1 (MQ=255)
                                                        tGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAttt                                                                                              <  2:4143411/150‑1 (MQ=255)
                                                         gAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCATATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtttt                                                                                             <  2:184602/150‑1 (MQ=255)
                                                          aGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtttt                                                                                             <  2:2677854/150‑1 (MQ=255)
                                                           gcggcgCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCt                                                                                           <  2:1790847/150‑1 (MQ=255)
                                                                 cTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGag                                                                                     <  1:15030/150‑1 (MQ=255)
                                                                  ttATGCGCCGTTTTTTTTTATGTAACACCATATTTAAATACCCAAATAATATTTGAGTCATCTCACAAAAAAGCATTTTAACTGGTAAAAAAATATTTAAATCGTGACAAAAAACCATTTGGGCAAAATACAGTTTATTTTATGGGGagt                                                                                    >  2:2430577/1‑150 (MQ=255)
                                                                   tATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagta                                                                                   <  1:1850256/150‑1 (MQ=255)
                                                                     tGCGCCGATTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagtagt                                                                                 <  2:2217403/150‑1 (MQ=255)
                                                                      gcgcCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagtagt                                                                                 <  1:393429/150‑1 (MQ=255)
                                                                        gcCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGtctc                                                                              <  2:4582574/150‑1 (MQ=255)
                                                                          cGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGtctctc                                                                            >  1:1940624/1‑150 (MQ=255)
                                                                           gTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCg                                                                           >  2:1377344/1‑150 (MQ=255)
                                                                           cgattttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGtctctc                                                                            <  2:1671125/147‑1 (MQ=255)
                                                                              ttttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGtt                                                                         <  1:3310999/150‑1 (MQ=255)
                                                                                ttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTc                                                                       <  2:2539999/150‑1 (MQ=255)
                                                                                ttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTc                                                                       <  2:5081946/150‑1 (MQ=255)
                                                                                     tCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATggg                                                                  <  1:2581275/150‑1 (MQ=255)
                                                                                     tCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATggg                                                                  <  2:2331338/150‑1 (MQ=255)
                                                                                       tGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGAc                                                                <  1:2795437/150‑1 (MQ=255)
                                                                                            cacCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTa                                                           >  2:4714329/1‑150 (MQ=255)
                                                                                             acCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTAc                                                          <  1:4831494/150‑1 (MQ=255)
                                                                                              ccTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTAcc                                                         >  2:4925721/1‑150 (MQ=255)
                                                                                              ccTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTAcc                                                         >  2:885498/1‑150 (MQ=255)
                                                                                              ccTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTAcc                                                         >  2:1584240/1‑150 (MQ=255)
                                                                                                ttATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCAc                                                       <  1:2094874/150‑1 (MQ=255)
                                                                                                   tttACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGa                                                    <  1:1631912/150‑1 (MQ=255)
                                                                                                    ttACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGaa                                                   <  1:3988387/150‑1 (MQ=255)
                                                                                                    ttACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGaa                                                   <  2:614997/150‑1 (MQ=255)
                                                                                                    ttACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGaa                                                   <  2:128314/150‑1 (MQ=255)
                                                                                                      aCATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAg                                                 <  2:2794066/150‑1 (MQ=255)
                                                                                                       cATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAgg                                                <  2:4334431/150‑1 (MQ=255)
                                                                                                         tCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCa                                              >  1:4136760/1‑150 (MQ=255)
                                                                                                           cccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAAc                                            >  1:1985757/1‑150 (MQ=255)
                                                                                                            ccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACa                                           <  1:1223072/150‑1 (MQ=255)
                                                                                                             cATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACAt                                          <  1:158551/150‑1 (MQ=255)
                                                                                                             cATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACAt                                          >  2:2004582/1‑150 (MQ=255)
                                                                                                             cATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACAt                                          >  2:3348192/1‑150 (MQ=255)
                                                                                                              aTAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAACGGCAACATg                                         >  2:4051490/1‑150 (MQ=255)
                                                                                                                aGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGaa                                       <  2:5018005/150‑1 (MQ=255)
                                                                                                                aGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGaa                                       <  2:5221856/150‑1 (MQ=255)
                                                                                                                aGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGaa                                       >  1:1715805/1‑150 (MQ=255)
                                                                                                                aGATAATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGAATTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATATTCTGCGGAGTAGACTCACGTTACATGTGACCGCTACCACGGAAAGGCAACATGaa                                       >  2:4893060/1‑150 (MQ=255)
                                                                                                                     attTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGa                                  <  1:905240/150‑1 (MQ=255)
                                                                                                                       ttGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaaa                                >  2:2107781/1‑150 (MQ=255)
                                                                                                                        tGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaaaa                               >  1:1121946/1‑150 (MQ=255)
                                                                                                                        tGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaaaa                               <  1:1055135/150‑1 (MQ=255)
                                                                                                                           gTCAGCTCACAACTACGCTTTTTCCCTAGTAACAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTa                            <  2:3514262/150‑1 (MQ=255)
                                                                                                                           gTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTa                            >  1:2648510/1‑150 (MQ=255)
                                                                                                                           gTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTa                            <  2:3001237/150‑1 (MQ=255)
                                                                                                                              aGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACgg                         <  2:293325/150‑1 (MQ=255)
                                                                                                                                     aaaTACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCt                  <  1:4969972/150‑1 (MQ=255)
                                                                                                                                      aaTACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCACCATGAAACAGAAAATTACGGATTACCTg                 >  2:3442212/1‑150 (MQ=255)
                                                                                                                                       aTACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTgg                <  1:4233697/150‑1 (MQ=255)
                                                                                                                                        tACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGa               <  1:2290983/150‑1 (MQ=255)
                                                                                                                                         aCGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGAc              >  2:3825197/1‑150 (MQ=255)
                                                                                                                                          cGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACg             >  2:2346657/1‑150 (MQ=255)
                                                                                                                                          cGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACg             <  2:528030/150‑1 (MQ=255)
                                                                                                                                            cTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGaa           >  1:556682/1‑150 (MQ=255)
                                                                                                                                               ttttCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAAt         <  1:2430577/150‑1 (MQ=255)
                                                                                                                                               ttttCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAAt         <  2:1940573/150‑1 (MQ=255)
                                                                                                                                                tttCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATc        >  2:1050695/1‑150 (MQ=255)
                                                                                                                                                 ttCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCt       >  1:4873551/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGT  >  CP009273/1364217‑1364514

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: