Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,555,325 T→C 18.8% intergenic (‑142/‑48) glpE ← / → glpD thiosulfate:cyanide sulfurtransferase (rhodanese)/sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,555,3250TC18.8% 776.7 / 172.3 322intergenic (‑142/‑48)glpE/glpDthiosulfate:cyanide sulfurtransferase (rhodanese)/sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (144/116);  new base C (36/24);  total (181/141)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 327 total aligned reads displayed.

TGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAG  >  CP009273/3555176‑3555471
                                                                                                                                                     |                                                                                                                                                  
 gATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGtt                                                                                                                                                   >  1:689970/1‑150 (MQ=255)
 gATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGtt                                                                                                                                                   >  2:3474380/1‑150 (MQ=255)
  aTCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTc                                                                                                                                                  <  1:4392505/150‑1 (MQ=255)
      aTGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATa                                                                                                                                              <  2:2902215/150‑1 (MQ=255)
       tGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATaa                                                                                                                                             >  2:4949258/1‑150 (MQ=255)
       tGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATaa                                                                                                                                             >  1:3074101/1‑150 (MQ=255)
       tGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATaa                                                                                                                                             >  1:1779616/1‑150 (MQ=255)
        gCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAAc                                                                                                                                            <  1:2242561/150‑1 (MQ=255)
        gCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAAc                                                                                                                                            <  2:2264033/150‑1 (MQ=255)
         cGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACg                                                                                                                                           >  1:797347/1‑150 (MQ=255)
          gTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGCTCGATAACGa                                                                                                                                          >  1:1737873/1‑150 (MQ=255)
           tctctcTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                                                                                                                                         <  2:3010005/150‑1 (MQ=255)
              ctctTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAt                                                                                                                                      >  2:1569185/1‑150 (MQ=255)
               tctttcttTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAtt                                                                                                                                     <  1:4431639/150‑1 (MQ=255)
               tctttcttTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGCTCGATAACGAACAtt                                                                                                                                     <  1:936227/150‑1 (MQ=255)
                ctttctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt                                                                                                                                    <  1:2915427/150‑1 (MQ=255)
                 tttctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa                                                                                                                                   >  1:502126/1‑150 (MQ=255)
                   tctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATg                                                                                                                                 <  1:4675915/150‑1 (MQ=255)
                    ctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGa                                                                                                                                <  1:4474368/150‑1 (MQ=255)
                    ctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGa                                                                                                                                <  2:475037/150‑1 (MQ=255)
                     tttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAg                                                                                                                               >  2:5078934/1‑150 (MQ=255)
                      ttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGTTAATATGCTCGATAACGAACATTTATGAGc                                                                                                                              <  2:3316125/150‑1 (MQ=255)
                        acaaacaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCtt                                                                                                                            >  1:3386149/1‑150 (MQ=255)
                         caaacaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCttt                                                                                                                           <  1:2345780/150‑1 (MQ=255)
                          aaacaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTa                                                                                                                          >  1:1344846/1‑150 (MQ=255)
                          aaacaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTa                                                                                                                          >  1:321900/1‑150 (MQ=255)
                          aaacaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATCCTTTGCTAATATGCTCGATAACGAACATTTATGAGCTTTa                                                                                                                          <  1:3707559/150‑1 (MQ=255)
                           aacaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTaa                                                                                                                         >  2:1831005/1‑150 (MQ=255)
                            acaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAAc                                                                                                                        <  1:521917/150‑1 (MQ=255)
                             caaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACg                                                                                                                       >  1:375909/1‑150 (MQ=255)
                                 tGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAg                                                                                                                   <  2:3777809/150‑1 (MQ=255)
                                     cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGaa                                                                                                               <  2:1985069/150‑1 (MQ=255)
                                      aaaTTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAAt                                                                                                              <  1:3265914/150‑1 (MQ=255)
                                       aaTTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATg                                                                                                             <  1:4824198/150‑1 (MQ=255)
                                       aaTTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGCTCGATAACGAACATTTATTAGCTTTAAAGAAAGTGAATg                                                                                                             >  1:252934/1‑150 (MQ=255)
                                        aTTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGa                                                                                                            >  2:2958111/1‑150 (MQ=255)
                                             ccGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTCATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGgca                                                                                                       >  2:3602308/1‑150 (MQ=255)
                                                cacaGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGgcagca                                                                                                    >  2:126307/1‑150 (MQ=255)
                                                cacaGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATCCTTTGCTAATATGCTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGgcagca                                                                                                    <  2:4236061/150‑1 (MQ=255)
                                                 acaGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCAt                                                                                                   <  1:3365045/150‑1 (MQ=255)
                                                      ttACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGaaa                                                                                              <  2:913106/150‑1 (MQ=255)
                                                      ttACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGaaa                                                                                              <  2:1751917/150‑1 (MQ=255)
                                                       tACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAAc                                                                                             >  2:616548/1‑150 (MQ=255)
                                                       tACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAAc                                                                                             >  1:2153475/1‑150 (MQ=255)
                                                       tACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAAc                                                                                             <  1:3527854/150‑1 (MQ=255)
                                                         cGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCa                                                                                           <  1:622328/150‑1 (MQ=255)
                                                             gAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGa                                                                                       <  1:4289808/150‑1 (MQ=255)
                                                              aaGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCACAGAt                                                                                      >  2:4662741/1‑150 (MQ=255)
                                                              aaGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGAt                                                                                      >  1:1296868/1‑150 (MQ=255)
                                                               aGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCTATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATc                                                                                     >  1:2884827/1‑150 (MQ=255)
                                                               aGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATc                                                                                     >  1:1064282/1‑150 (MQ=255)
                                                                 cGGCAGATAAACGCCATAATGTTATACATATCACTCTTAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                   <  2:4034936/150‑1 (MQ=255)
                                                                 cGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                   <  2:1423458/150‑1 (MQ=255)
                                                                 cGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                   <  1:3254504/150‑1 (MQ=255)
                                                                 cGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                   <  1:3884/150‑1 (MQ=255)
                                                                      gATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATtgt                                                                              <  1:1579468/150‑1 (MQ=255)
                                                                        tAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGa                                                                            >  1:2325028/1‑150 (MQ=255)
                                                                        tAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGCTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAACTTCTGATGGTGa                                                                            >  2:4831712/1‑150 (MQ=255)
                                                                         aaaCGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGAt                                                                           <  2:4925626/150‑1 (MQ=255)
                                                                         aaaCGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGCTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGAt                                                                           <  1:1340406/150‑1 (MQ=255)
                                                                         aaaCGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGCTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGAt                                                                           >  1:1963742/1‑150 (MQ=255)
                                                                           aCGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAg                                                                         <  2:3637031/150‑1 (MQ=255)
                                                                              ccATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAgggg                                                                      >  1:4098160/1‑150 (MQ=255)
                                                                              ccATAATGTTATAAATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAgggg                                                                      <  2:2509287/150‑1 (MQ=255)
                                                                                aTAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGggc                                                                    >  2:664387/1‑150 (MQ=255)
                                                                                    tGTTATACATATCACTCTTAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCa                                                                <  2:2330238/150‑1 (MQ=255)
                                                                                       tataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCa                                                             >  1:4186677/1‑150 (MQ=255)
                                                                                         taaatatCACTCTAAAAAGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATTTTAGAAAACTAAAAATTATTACATTAAACGAAAGTGAATTAGCGCATAATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAAt                                                           >  2:1778951/4‑150 (MQ=255)
                                                                                           cATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGCTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATgg                                                         >  2:399297/1‑150 (MQ=255)
                                                                                            atatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGCTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGt                                                        >  1:4007609/1‑150 (MQ=255)
                                                                                             tatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTg                                                       >  2:2127204/1‑150 (MQ=255)
                                                                                             tatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTg                                                       >  1:3183500/1‑150 (MQ=255)
                                                                                             tatCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTg                                                       <  1:1426488/150‑1 (MQ=255)
                                                                                              atCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGc                                                      >  2:247934/1‑150 (MQ=255)
                                                                                              atCAATCTAAAATGTATTTTCAATGTTACCTAAAGCGAGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAAAGAAAGTGAATGAGGGCAGCAAGGAAAACAAAGATCTGAATGTGATAGTGGGCGTCATCAATGGAGc                                                      >  2:1487668/1‑150 (MQ=255)
                                                                                                cACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGAGTGTGATAGGGGGCGGCATCAATGGTGCTg                                                    >  2:45649/1‑150 (MQ=255)
                                                                                                    ctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATTTTCGATAACGAACAAATATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACAAAAGATCTGAGTGTGATATGGGGCGGCATAAATGTTGCTGGTaa                                                >  1:2752886/1‑149 (MQ=255)
                                                                                                               tttCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAcg                                     <  1:2485185/150‑1 (MQ=255)
                                                                                                               tttCAATGTTACCTAAAGCGCGATCCTTTGCTAATATGCTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAcg                                     >  1:4696946/1‑150 (MQ=255)
                                                                                                                 tCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAcgcc                                   <  2:639872/150‑1 (MQ=255)
                                                                                                                 tCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAcgcc                                   <  2:5141314/150‑1 (MQ=255)
                                                                                                                  cAATGTTACCTAAAGCGCGATTCTTTGCTAATATGCTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAcgccg                                  <  1:1507205/150‑1 (MQ=255)
                                                                                                                    aTGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAACGATCTGATGGTGAGAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCt                                >  2:2151596/1‑150 (MQ=255)
                                                                                                                    aTGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCt                                >  2:2647030/1‑150 (MQ=255)
                                                                                                                      gTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTgg                              <  2:799455/150‑1 (MQ=255)
                                                                                                                              aaGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGtt                      >  1:587211/1‑150 (MQ=255)
                                                                                                                              aaGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGtt                      >  1:2963105/1‑150 (MQ=255)
                                                                                                                              aaGCGCGATTCTTTGCTAATATGCTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGtt                      >  1:725239/1‑150 (MQ=255)
                                                                                                                               aGCGCGATTCTTTTCTAATATGCTCGATAACTAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGttt                     <  2:2735955/150‑1 (MQ=255)
                                                                                                                               aGCGCGATTCTTTGCTAATATGATCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGttt                     <  2:2663180/150‑1 (MQ=255)
                                                                                                                                 cgcgATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTAt                   <  2:4318913/150‑1 (MQ=255)
                                                                                                                                   ataTTCTTTGCTAATATGTTAGATAAAGAACATTTTTGAGCTTTAACGATAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGTTAGGGGGCGGCATAAATGGTGCTGGTATCGCGGCAGACGCCGCTGGTCGCGGTTTATcc                 <  2:4692207/148‑1 (MQ=255)
                                                                                                                                    gATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCg                >  2:931172/1‑150 (MQ=255)
                                                                                                                                    gATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCg                >  2:4430305/1‑150 (MQ=255)
                                                                                                                                    gATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCg                >  2:3111352/1‑150 (MQ=255)
                                                                                                                                     aTTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTGTATCCGt               >  2:3856490/1‑150 (MQ=255)
                                                                                                                                        cTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCt            >  1:2750708/1‑150 (MQ=255)
                                                                                                                                            gCTAATATGCTCGAGAACGAACATTTATGAGCTTTACCGAAAGTGAATGAGGGCAGCATGGAATCCAAAGATCTGAGTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGTTg        <  1:3139580/150‑1 (MQ=255)
                                                                                                                                             cTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGc       <  1:1226903/150‑1 (MQ=255)
                                                                                                                                                atatGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGTTAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCTGACGCCGCTGGACGCGGTTTATCCGTGCTGATGATgg    >  2:2819857/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGATCCATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAG  >  CP009273/3555176‑3555471

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: