Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,647,590 C→A 8.2% E235* (GAA→TAA)  sseB ← rhodanase‑like enzyme, sulfur transfer from thiosulfate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,647,5900CA 8.2% 1066.1 / 47.4 287E235* (GAA→TAA) sseBrhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base C (137/126);  new base A (15/9);  total (152/135)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 293 total aligned reads displayed.

GGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAAT  >  CP009273/2647441‑2647739
                                                                                                                                                     |                                                                                                                                                     
ggCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTc                                                                                                                                                       <  1:5055115/150‑1 (MQ=255)
     cAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaa                                                                                                                                                  >  1:3493433/1‑150 (MQ=255)
     cAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTAGGtaa                                                                                                                                                  <  2:4475023/150‑1 (MQ=255)
        tttACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaataa                                                                                                                                               >  1:3593931/1‑150 (MQ=255)
        tttACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaataa                                                                                                                                               <  1:5113680/150‑1 (MQ=255)
          tACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAg                                                                                                                                             >  1:4780330/1‑150 (MQ=255)
            cTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTg                                                                                                                                           <  1:54595/150‑1 (MQ=255)
             tggtggTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTgg                                                                                                                                          <  1:43985/150‑1 (MQ=255)
               gtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCACTATGTTCGGTAATAAAGTGGCt                                                                                                                                        >  2:3014654/1‑150 (MQ=255)
               gtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCt                                                                                                                                        <  1:1487149/150‑1 (MQ=255)
                   ttAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAAtt                                                                                                                                    <  2:1179703/150‑1 (MQ=255)
                    tAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTc                                                                                                                                   <  2:420612/150‑1 (MQ=255)
                     aacaacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTcc                                                                                                                                  >  1:2058689/1‑150 (MQ=255)
                      acaacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCt                                                                                                                                 <  2:4716172/150‑1 (MQ=255)
                       caacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCtt                                                                                                                                <  1:182314/150‑1 (MQ=255)
                       caacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCtt                                                                                                                                <  1:2946321/150‑1 (MQ=255)
                        aacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTAACGCAAAAAAAAACCCCAGCGAAGTTCATAGAATGGCGCAATATGTTAGGTAATAAAGTTGCTAATTCCttt                                                                                                                               >  1:2100644/1‑150 (MQ=255)
                          ccTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCttttt                                                                                                                             <  2:2227660/150‑1 (MQ=255)
                                gagCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtt                                                                                                                       <  2:904851/150‑1 (MQ=255)
                                  gcagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtttt                                                                                                                     <  2:4581357/150‑1 (MQ=255)
                                  gcagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtttt                                                                                                                     <  1:3557779/150‑1 (MQ=255)
                                   cagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCttttt                                                                                                                    <  1:3636902/150‑1 (MQ=255)
                                    agcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtttttt                                                                                                                   >  1:2672263/1‑150 (MQ=255)
                                     gcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTCTCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTCATAAAGTGGCTAATTCCTGTTTCCCCTTTTTTc                                                                                                                  >  2:3439436/1‑150 (MQ=255)
                                      caAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCa                                                                                                                 <  2:1075494/150‑1 (MQ=255)
                                       aaGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCAc                                                                                                                >  2:4024222/1‑150 (MQ=255)
                                        aGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCAcc                                                                                                               <  1:2850451/150‑1 (MQ=255)
                                        aGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCAcc                                                                                                               <  1:1582015/150‑1 (MQ=255)
                                         gTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCt                                                                                                              >  2:1533155/1‑150 (MQ=255)
                                              cccGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTACTAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACaa                                                                                                         >  2:2273379/1‑150 (MQ=255)
                                               ccGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACaaa                                                                                                        >  1:1004394/1‑150 (MQ=255)
                                                cGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAat                                                                                                       <  1:4602204/150‑1 (MQ=255)
                                                cGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTAGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAat                                                                                                       <  2:4672978/150‑1 (MQ=255)
                                                  aaTGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAatat                                                                                                     >  2:4789038/1‑150 (MQ=255)
                                                     gccgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGa                                                                                                  >  1:1036682/1‑150 (MQ=255)
                                                     gccgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGa                                                                                                  >  1:3651832/1‑150 (MQ=255)
                                                      ccgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGAt                                                                                                 <  2:4905202/150‑1 (MQ=255)
                                                       cgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGAtt                                                                                                <  2:2026908/150‑1 (MQ=255)
                                                          ccTTTGGGCGGCATATTTTAGATTATCCGCTTCTGTTTAAAGTCACGCAATAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAATTATCGATTGGt                                                                                             <  2:635224/150‑1 (MQ=255)
                                                          ccTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGt                                                                                             >  2:3320341/1‑150 (MQ=255)
                                                            tttGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTc                                                                                           >  1:132786/1‑150 (MQ=255)
                                                               gggCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATc                                                                                        <  2:1296373/150‑1 (MQ=255)
                                                                  cggcATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATGGGTTCATCGcc                                                                                     >  1:3118909/1‑150 (MQ=255)
                                                                   ggcATATTTTAGATTCTCCTATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTAGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCa                                                                                    <  2:3880725/150‑1 (MQ=255)
                                                                   ggcATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCa                                                                                    >  2:1718190/1‑150 (MQ=255)
                                                                    gcATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATATAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTAAACTGACAAATATCGATTGGTTCATCGCCAg                                                                                   >  1:142503/1‑150 (MQ=255)
                                                                     cATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAgg                                                                                  <  1:1257652/150‑1 (MQ=255)
                                                                     cATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAgg                                                                                  >  1:794301/1‑150 (MQ=255)
                                                                       tatTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTa                                                                                >  2:3750117/1‑150 (MQ=255)
                                                                       tatTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTa                                                                                <  2:2591187/150‑1 (MQ=255)
                                                                       tatTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTa                                                                                >  2:3373088/1‑150 (MQ=255)
                                                                        atTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTaa                                                                               >  2:1033194/1‑150 (MQ=255)
                                                                        atTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTaa                                                                               >  1:3764520/1‑150 (MQ=255)
                                                                            tAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTa                                                                           >  2:1806755/1‑150 (MQ=255)
                                                                            tAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTAGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTa                                                                           >  1:2743441/1‑150 (MQ=255)
                                                                              gATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATc                                                                         <  2:1630673/150‑1 (MQ=255)
                                                                               aTTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCg                                                                        >  2:2904145/1‑150 (MQ=255)
                                                                                ttATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCgg                                                                       <  2:1404716/150‑1 (MQ=255)
                                                                                ttATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTAGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCgg                                                                       <  1:1658921/150‑1 (MQ=255)
                                                                                     cGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTACTGTATCGGTCGCt                                                                  >  2:76/1‑150 (MQ=255)
                                                                                      gATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTa                                                                 >  1:3408536/1‑150 (MQ=255)
                                                                                        ttCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTaca                                                               <  1:3307826/150‑1 (MQ=255)
                                                                                         tCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTacac                                                              >  1:2976511/1‑150 (MQ=255)
                                                                                           tGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACtt                                                            >  2:1455099/1‑150 (MQ=255)
                                                                                           tGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACtt                                                            >  1:4286831/1‑150 (MQ=255)
                                                                                             tttAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTcc                                                          <  2:3636642/150‑1 (MQ=255)
                                                                                               tAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTACTGTATCGGTCGCTACACTTCCCg                                                        >  2:1106749/1‑150 (MQ=255)
                                                                                               tAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCg                                                        >  2:2703171/1‑150 (MQ=255)
                                                                                               tAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCg                                                        <  2:2374835/150‑1 (MQ=255)
                                                                                                aaaGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGt                                                       >  1:4024252/1‑150 (MQ=255)
                                                                                                 aaGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTc                                                      <  1:1957132/150‑1 (MQ=255)
                                                                                                 aaGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTc                                                      >  1:3322519/1‑150 (MQ=255)
                                                                                                     cACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCt                                                  <  2:1993003/150‑1 (MQ=255)
                                                                                                       cGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGa                                                >  1:3992909/1‑150 (MQ=255)
                                                                                                          aaaaaaCCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGaata                                             <  1:4171903/150‑1 (MQ=255)
                                                                                                          aaaaaaCCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGaata                                             >  1:4387981/1‑150 (MQ=255)
                                                                                                           aaaaaCCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGaataa                                            >  2:4343489/1‑150 (MQ=255)
                                                                                                           aaaaaCCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGaataa                                            >  2:4549394/1‑150 (MQ=255)
                                                                                                            aaaaCCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGaataat                                           <  2:3812062/150‑1 (MQ=255)
                                                                                                            aaaaCCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGaataat                                           >  1:2488011/1‑150 (MQ=255)
                                                                                                                 caccCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATttctt                                      >  1:730505/1‑150 (MQ=255)
                                                                                                                 caccCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATttctt                                      >  2:1325632/1‑150 (MQ=255)
                                                                                                                  accCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATttcttc                                     <  2:3117149/150‑1 (MQ=255)
                                                                                                                   ccccAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCg                                    <  2:2673501/150‑1 (MQ=255)
                                                                                                                      cAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGata                                 >  2:5024254/1‑150 (MQ=255)
                                                                                                                        gCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATc                               <  1:2935056/150‑1 (MQ=255)
                                                                                                                          gACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCAc                             >  2:3570471/1‑150 (MQ=255)
                                                                                                                          gACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCAc                             >  2:2124627/1‑150 (MQ=255)
                                                                                                                            cGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCa                           <  2:1036682/150‑1 (MQ=255)
                                                                                                                             gTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCAt                          >  1:2629515/1‑150 (MQ=255)
                                                                                                                              ttCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATc                         >  1:835035/1‑150 (MQ=255)
                                                                                                                                cATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCgg                       >  2:3989958/1‑150 (MQ=255)
                                                                                                                                 aTAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGc                      >  2:390538/1‑150 (MQ=255)
                                                                                                                                  tAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCt                     >  1:2789621/1‑150 (MQ=255)
                                                                                                                                      aTGGCGCAATATGTTAGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCaa                 >  1:3549965/1‑150 (MQ=255)
                                                                                                                                           gcAATATGTTCGGTAATAAAGTGGCTAATTCCTTCTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCa            >  1:1122308/1‑150 (MQ=255)
                                                                                                                                                atGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAg       <  2:411103/150‑1 (MQ=255)
                                                                                                                                                atGTTAGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAg       >  1:3181938/1‑150 (MQ=255)
                                                                                                                                                 tGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGt      <  1:989068/150‑1 (MQ=255)
                                                                                                                                                     cGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAAt  >  2:2984300/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCATCGCCAGGTAATGTATCGGTCGCTACACTTCCCGTCGCCTGAATAATTTCTTCGATATCACCATCGGCTTCAATGCCAATAAGTAAAT  >  CP009273/2647441‑2647739

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