Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,034,714 C→T 12.0% intergenic (+115/‑142) yedZ → / → zinT inner membrane heme subunit for periplasmic YedYZ reductase/zinc and cadmium binding protein, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,034,7140CT12.0% 860.6 / 85.9 286intergenic (+115/‑142)yedZ/zinTinner membrane heme subunit for periplasmic YedYZ reductase/zinc and cadmium binding protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (113/138);  new base T (17/17);  total (130/156)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 300 total aligned reads displayed.

TACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGAT  >  CP009273/2034565‑2034863
                                                                                                                                                     |                                                                                                                                                     
tACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCAcc                                                                                                                                                       <  2:1236934/150‑1 (MQ=255)
   gCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACTaaa                                                                                                                                                    <  1:1500590/150‑1 (MQ=255)
   gCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCaaa                                                                                                                                                    <  1:1454340/150‑1 (MQ=255)
     taGCAGGTTCAAAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAaat                                                                                                                                                  <  1:4092714/149‑1 (MQ=255)
     aaGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAaat                                                                                                                                                  <  1:1840452/150‑1 (MQ=255)
      aGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAaata                                                                                                                                                 <  2:2989457/150‑1 (MQ=255)
       gCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAaataa                                                                                                                                                >  1:2136359/1‑150 (MQ=255)
        cAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAaataat                                                                                                                                               <  2:3339124/150‑1 (MQ=255)
          ggTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAattt                                                                                                                                             <  2:129141/150‑1 (MQ=255)
          ggTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAattt                                                                                                                                             >  1:2592056/1‑150 (MQ=255)
              cacaACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAatttattt                                                                                                                                         <  2:2415454/150‑1 (MQ=255)
              cacaACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAatttattt                                                                                                                                         >  2:2603928/1‑150 (MQ=255)
               acaacaAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAAtttatttt                                                                                                                                        >  1:3984805/1‑150 (MQ=255)
                caacatATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAAtttatttta                                                                                                                                       <  2:2664959/150‑1 (MQ=255)
                   caaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTtattttattt                                                                                                                                    <  2:3412606/150‑1 (MQ=255)
                    aaaTTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACTAAAATAATTTATTTTATTTc                                                                                                                                   <  2:351943/150‑1 (MQ=255)
                    aaaTTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTc                                                                                                                                   >  1:3994110/1‑150 (MQ=255)
                        tATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATa                                                                                                                               <  2:2764969/150‑1 (MQ=255)
                         aTCTGTATAACTTCATAGTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATaa                                                                                                                              <  2:3994683/150‑1 (MQ=255)
                          tatGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAag                                                                                                                             <  1:4156231/148‑1 (MQ=255)
                          tCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAag                                                                                                                             <  1:3334503/150‑1 (MQ=255)
                          tCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAag                                                                                                                             <  1:732099/150‑1 (MQ=255)
                          tCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAag                                                                                                                             <  2:607882/150‑1 (MQ=255)
                           cTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAaga                                                                                                                            <  2:4164792/150‑1 (MQ=255)
                           cTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAaga                                                                                                                            <  2:3832101/150‑1 (MQ=255)
                           cTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAaga                                                                                                                            >  2:962698/1‑150 (MQ=255)
                            tGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagat                                                                                                                           <  2:1712425/150‑1 (MQ=255)
                             gTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagata                                                                                                                          >  2:2736490/1‑150 (MQ=255)
                              tataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagatag                                                                                                                         <  1:1545762/150‑1 (MQ=255)
                                 aaCTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAgatagata                                                                                                                      <  1:4211955/150‑1 (MQ=255)
                                    ttCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGt                                                                                                                   >  2:3537370/1‑150 (MQ=255)
                                     tCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTa                                                                                                                  <  2:2229992/150‑1 (MQ=255)
                                     tCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTa                                                                                                                  >  2:2365897/1‑150 (MQ=255)
                                     tCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTa                                                                                                                  >  2:3812664/1‑150 (MQ=255)
                                     tCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTAAGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTa                                                                                                                  >  2:1542138/1‑150 (MQ=255)
                                      cATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAg                                                                                                                 >  1:2689191/1‑150 (MQ=255)
                                       atatTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACTAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGa                                                                                                                >  1:602213/1‑150 (MQ=255)
                                       atatTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGa                                                                                                                <  2:888919/150‑1 (MQ=255)
                                        tatTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGaa                                                                                                               >  2:1203158/1‑150 (MQ=255)
                                         atTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACTAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAAc                                                                                                              >  2:1921453/1‑150 (MQ=255)
                                          ttGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACTAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACt                                                                                                             >  1:2109503/1‑150 (MQ=255)
                                          ttGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACt                                                                                                             <  2:91596/150‑1 (MQ=255)
                                          ttGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACt                                                                                                             <  1:924190/150‑1 (MQ=255)
                                           tGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTg                                                                                                            >  2:2046401/1‑150 (MQ=255)
                                                 cAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAg                                                                                                      <  1:2330095/150‑1 (MQ=255)
                                                  aaaGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGc                                                                                                     >  1:3828916/1‑150 (MQ=255)
                                                     gTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCAt                                                                                                  <  1:982627/150‑1 (MQ=255)
                                                     gTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCAt                                                                                                  <  1:3309077/150‑1 (MQ=255)
                                                       tcgctTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCTCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGc                                                                                                <  1:2816150/149‑1 (MQ=255)
                                                       aCGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGc                                                                                                <  1:2816353/150‑1 (MQ=255)
                                                        cgcgTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGct                                                                                               >  1:4334319/1‑150 (MQ=255)
                                                        cgcgTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGct                                                                                               <  1:4012177/150‑1 (MQ=255)
                                                        cgcgTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGct                                                                                               >  1:3434948/1‑150 (MQ=255)
                                                        cgcgTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGct                                                                                               >  2:1738363/1‑150 (MQ=255)
                                                           gTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCtctc                                                                                            <  1:2331162/150‑1 (MQ=255)
                                                           gTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCtctc                                                                                            <  2:214344/150‑1 (MQ=255)
                                                           gTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGAAGAACAGAGAAAGCCATGctcac                                                                                            >  1:4499196/1‑148 (MQ=255)
                                                              aTCCCTGCCGGATGTGGCTCGAGCGCCCTATCCGTCATACACAGATCGCATCAGTTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGtt                                                                                         <  1:2914869/150‑1 (MQ=255)
                                                               tCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGttt                                                                                        <  2:2943666/150‑1 (MQ=255)
                                                                 ccTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTcc                                                                                      <  1:1018884/150‑1 (MQ=255)
                                                                        gATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAgt                                                                               <  1:2477305/150‑1 (MQ=255)
                                                                        gATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAgt                                                                               <  1:3760020/150‑1 (MQ=255)
                                                                        gATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAgt                                                                               <  1:4406837/150‑1 (MQ=255)
                                                                                  gAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCAttt                                                                     <  1:2256163/150‑1 (MQ=255)
                                                                                   aGCCCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCAtttt                                                                    >  2:4810/1‑150 (MQ=255)
                                                                                    gcgcCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTg                                                                   <  1:1847939/150‑1 (MQ=255)
                                                                                    gcgcCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTg                                                                   <  2:1482663/150‑1 (MQ=255)
                                                                                    gcgcCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTg                                                                   >  1:4363960/1‑150 (MQ=255)
                                                                                      gcCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCt                                                                 <  1:1585004/150‑1 (MQ=255)
                                                                                             ccGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACTAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGtt                                                          <  1:2226733/150‑1 (MQ=255)
                                                                                              cGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTa                                                         <  1:3849136/150‑1 (MQ=255)
                                                                                               ggCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTAc                                                        <  1:962698/150‑1 (MQ=255)
                                                                                                 cATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACaa                                                      <  1:1964713/150‑1 (MQ=255)
                                                                                                         gATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACAt                                              <  2:1980734/150‑1 (MQ=255)
                                                                                                         gATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACAt                                              <  2:1642587/150‑1 (MQ=255)
                                                                                                         gATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACAt                                              <  2:1642224/150‑1 (MQ=255)
                                                                                                           tCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACTAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTa                                            <  1:1753743/150‑1 (MQ=255)
                                                                                                           tCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTa                                            >  1:3052829/1‑150 (MQ=255)
                                                                                                            cGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTac                                           >  2:4099684/1‑150 (MQ=255)
                                                                                                               aTCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTAcaca                                        >  1:2890268/1‑150 (MQ=255)
                                                                                                                   ggTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACAcatcat                                    <  2:3153888/150‑1 (MQ=255)
                                                                                                                   ggTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACAcatcat                                    <  2:1780668/150‑1 (MQ=255)
                                                                                                                         gTTTCGTTTGTCATCAACCTCAAATCACTAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACAt                              <  1:1790549/150‑1 (MQ=255)
                                                                                                                         gTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACAt                              >  1:804515/1‑150 (MQ=255)
                                                                                                                                ttGTCATCAACCTCAAATCACTAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAActct                       >  2:3287850/1‑150 (MQ=255)
                                                                                                                                   tcatcaACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTgga                    >  2:3715661/1‑150 (MQ=255)
                                                                                                                                      tcaACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTggagga                 <  2:2119093/150‑1 (MQ=255)
                                                                                                                                         aCCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAAc              <  2:3277291/150‑1 (MQ=255)
                                                                                                                                          ccTCAAATCACTAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACt             <  2:2885252/150‑1 (MQ=255)
                                                                                                                                           cTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTg            <  2:3768664/150‑1 (MQ=255)
                                                                                                                                             cAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGtt          >  1:3676017/1‑150 (MQ=255)
                                                                                                                                              aaaTCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGttt         >  2:3344754/1‑150 (MQ=255)
                                                                                                                                               aaTCACTAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACATATCATATACATTAACTCTGGAGGAAACTGtttt        >  1:1132617/1‑150 (MQ=255)
                                                                                                                                                aTCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTg       <  1:1842170/150‑1 (MQ=255)
                                                                                                                                                  cACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGc     <  2:3994007/150‑1 (MQ=255)
                                                                                                                                                  cACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGc     >  1:3914535/1‑150 (MQ=255)
                                                                                                                                                  cACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGc     >  2:4320107/1‑150 (MQ=255)
                                                                                                                                                    ccAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGa   <  1:2559533/150‑1 (MQ=255)
                                                                                                                                                     cAAAATACTTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGAt  <  2:1564548/150‑1 (MQ=255)
                                                                                                                                                     cAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGAt  >  2:132477/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGAT  >  CP009273/2034565‑2034863

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: