Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,193,578 G→T 6.3% D111Y (GAT→TAT)  zraR → fused DNA‑binding response regulator in two‑component regulatory system with ZraS: response regulator/sigma54 interaction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,193,5780GT 6.3% 1206.3 / 34.9 314D111Y (GAT→TAT) zraRfused DNA‑binding response regulator in two‑component regulatory system with ZraS: response regulator/sigma54 interaction protein
Reads supporting (aligned to +/- strand):  ref base G (148/146);  new base T (5/15);  total (153/161)
Fisher's exact test for biased strand distribution p-value = 3.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 320 total aligned reads displayed.

AGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCTCGTCGCGCCA  >  CP009273/4193429‑4193727
                                                                                                                                                     |                                                                                                                                                     
agaTGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCg                                                                                                                                                       <  2:4576715/150‑1 (MQ=255)
agaTGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCg                                                                                                                                                       <  2:1961541/150‑1 (MQ=255)
     gACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAAc                                                                                                                                                  >  1:196264/1‑150 (MQ=255)
       cGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCt                                                                                                                                                <  2:4513037/150‑1 (MQ=255)
        ggCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTa                                                                                                                                               <  2:4331548/150‑1 (MQ=255)
        ggCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTa                                                                                                                                               >  1:277197/1‑150 (MQ=255)
           aTCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAg                                                                                                                                            >  1:2637132/1‑150 (MQ=255)
             cGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGc                                                                                                                                          <  2:2131232/150‑1 (MQ=255)
              gCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCg                                                                                                                                         >  2:2910456/1‑150 (MQ=255)
              gCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCg                                                                                                                                         <  2:3297404/150‑1 (MQ=255)
              gCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCg                                                                                                                                         <  2:2222308/150‑1 (MQ=255)
                cACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGAc                                                                                                                                       >  1:3804112/1‑150 (MQ=255)
                 aCGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACg                                                                                                                                      >  2:3401849/1‑150 (MQ=255)
                 aCGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACg                                                                                                                                      >  1:2914091/1‑150 (MQ=255)
                      gAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGa                                                                                                                                 >  1:4576411/1‑150 (MQ=255)
                      gAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGa                                                                                                                                 <  2:2713209/150‑1 (MQ=255)
                       aaaGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGaa                                                                                                                                >  1:1875948/1‑150 (MQ=255)
                         agagATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGaaaa                                                                                                                              <  1:2930811/150‑1 (MQ=255)
                          gagaTCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGaaaaa                                                                                                                             >  2:1557142/1‑150 (MQ=255)
                             aTCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAgcg                                                                                                                          <  2:4295204/150‑1 (MQ=255)
                              tCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAgcgc                                                                                                                         >  1:1403510/1‑150 (MQ=255)
                               cAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCt                                                                                                                        >  2:2463593/1‑150 (MQ=255)
                                aaaGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTc                                                                                                                       <  2:4015937/150‑1 (MQ=255)
                                aaaGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTc                                                                                                                       >  1:4085401/1‑150 (MQ=255)
                                 aaGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCGGGAAAAAGCGCTcg                                                                                                                      >  2:3577571/1‑150 (MQ=255)
                                  aGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTcgc                                                                                                                     >  2:3957554/1‑150 (MQ=255)
                                      ttAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCAt                                                                                                                 <  1:2273218/150‑1 (MQ=255)
                                      ttAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCAt                                                                                                                 >  2:3126795/1‑150 (MQ=255)
                                      ttAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCAt                                                                                                                 <  2:3768145/150‑1 (MQ=255)
                                      ttAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCAt                                                                                                                 <  2:1428852/150‑1 (MQ=255)
                                      ctcaccccGGCCATTCCGGTCCTGATTATGACTGCGTACTCCAGCGTCGAGCCGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTTCAGGCGACGCTGGAAAAAGCGCTCGCGCAt                                                                                                                 <  1:3368677/145‑1 (MQ=255)
                                         aaCCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACg                                                                                                              <  2:4053504/150‑1 (MQ=255)
                                              ggCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAg                                                                                                         <  1:1262339/150‑1 (MQ=255)
                                              ggCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAg                                                                                                         <  1:4293790/150‑1 (MQ=255)
                                              ggCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAg                                                                                                         >  1:1717771/1‑150 (MQ=255)
                                               gCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGt                                                                                                        <  1:211475/150‑1 (MQ=255)
                                               gCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGt                                                                                                        >  2:2589779/1‑150 (MQ=255)
                                                 aaTTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTAt                                                                                                      <  1:3961401/150‑1 (MQ=255)
                                                  aTTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTAtt                                                                                                     <  1:4009273/150‑1 (MQ=255)
                                                  aTTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTCtt                                                                                                     >  1:1376452/1‑150 (MQ=255)
                                                       ggTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGc                                                                                                <  1:1995485/150‑1 (MQ=255)
                                                          gCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGa                                                                                             <  2:756087/150‑1 (MQ=255)
                                                           cTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGaa                                                                                            >  1:820403/1‑150 (MQ=255)
                                                              aTTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAaca                                                                                         <  2:4089087/150‑1 (MQ=255)
                                                                      tGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGt                                                                                 >  2:2313157/1‑150 (MQ=255)
                                                                      tGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGt                                                                                 >  2:4610431/1‑150 (MQ=255)
                                                                        cGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGa                                                                               >  1:71247/1‑150 (MQ=255)
                                                                        cGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGa                                                                               >  1:576279/1‑150 (MQ=255)
                                                                         gTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGAc                                                                              <  1:962676/150‑1 (MQ=255)
                                                                          tACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACt                                                                             >  2:1638854/1‑150 (MQ=255)
                                                                           aCTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTg                                                                            >  2:3196507/1‑150 (MQ=255)
                                                                            cTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTgc                                                                           >  2:2969127/1‑150 (MQ=255)
                                                                                aGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTgccagc                                                                       >  2:652983/1‑150 (MQ=255)
                                                                                   gTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGccagccag                                                                    <  2:4229821/150‑1 (MQ=255)
                                                                                      gagaCGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTc                                                                 <  2:404348/150‑1 (MQ=255)
                                                                                       agaCGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCg                                                                >  1:1300912/1‑150 (MQ=255)
                                                                                          cggcggTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTa                                                             <  2:1666867/150‑1 (MQ=255)
                                                                                            gcggTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATg                                                           <  2:228189/150‑1 (MQ=255)
                                                                                              ggTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGt                                                         >  1:3848815/1‑150 (MQ=255)
                                                                                               gTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTc                                                        <  1:1645306/150‑1 (MQ=255)
                                                                                               gTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTc                                                        <  2:872762/150‑1 (MQ=255)
                                                                                                 agagGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCgg                                                      >  2:594876/1‑150 (MQ=255)
                                                                                                   agGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTa                                                    <  1:4434960/150‑1 (MQ=255)
                                                                                                         tGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGcc                                              >  2:1883853/1‑150 (MQ=255)
                                                                                                         tGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGcc                                              >  1:3250267/1‑150 (MQ=255)
                                                                                                            aaaCTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCgg                                           >  2:1114912/1‑150 (MQ=255)
                                                                                                              aCTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCg                                         <  1:366062/150‑1 (MQ=255)
                                                                                                                 ggggCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATg                                      <  2:4488488/150‑1 (MQ=255)
                                                                                                                 ggggCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATg                                      <  2:927135/150‑1 (MQ=255)
                                                                                                                 ggggCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATg                                      >  1:1573735/1‑150 (MQ=255)
                                                                                                                 ggggCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATg                                      <  1:1788994/150‑1 (MQ=255)
                                                                                                                  gggCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGc                                     >  2:82790/1‑150 (MQ=255)
                                                                                                                      gcTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAaca                                 <  1:1216305/150‑1 (MQ=255)
                                                                                                                         ggATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCt                              <  2:2023192/150‑1 (MQ=255)
                                                                                                                         ggATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCt                              <  1:355833/150‑1 (MQ=255)
                                                                                                                            ttATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCt                           <  2:1585839/150‑1 (MQ=255)
                                                                                                                              aTCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCa                         >  1:1012107/1‑150 (MQ=255)
                                                                                                                              aTCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCa                         <  2:1108373/150‑1 (MQ=255)
                                                                                                                               tctcATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAg                        >  1:192418/1‑150 (MQ=255)
                                                                                                                                ctcATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGt                       >  2:3850584/1‑150 (MQ=255)
                                                                                                                                  catcaAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGa                     >  2:3354814/1‑150 (MQ=255)
                                                                                                                                    tcaAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGaaa                   >  1:1936518/1‑150 (MQ=255)
                                                                                                                                    tcaAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGaaa                   >  1:4263399/1‑150 (MQ=255)
                                                                                                                                     caAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAAt                  >  2:3012722/1‑150 (MQ=255)
                                                                                                                                        gCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGc               <  1:1172336/150‑1 (MQ=255)
                                                                                                                                        gCCGCTCGATATCGAAAACCTACAGGCGACGCTGTAAAAATCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGAGATGCAACACCTGCTCAGTGAAATCGc               >  2:480179/1‑150 (MQ=255)
                                                                                                                                         ccGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGcc              >  1:2257706/1‑150 (MQ=255)
                                                                                                                                         ccGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGcc              >  1:503393/1‑150 (MQ=255)
                                                                                                                                           gCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCt            <  2:3678049/150‑1 (MQ=255)
                                                                                                                                            cTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtc           >  1:4630666/1‑150 (MQ=255)
                                                                                                                                            cTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtc           >  2:3894319/1‑150 (MQ=255)
                                                                                                                                             tGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATTCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtcg          >  2:3099981/1‑150 (MQ=255)
                                                                                                                                             tGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtcg          <  2:1253798/150‑1 (MQ=255)
                                                                                                                                              ggATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtcgt         <  2:430706/150‑1 (MQ=255)
                                                                                                                                                aTTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtcgtcg       >  1:243366/1‑150 (MQ=255)
                                                                                                                                                aTTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtcgtcg       >  1:1507476/1‑150 (MQ=255)
                                                                                                                                                  ttCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCTCGTcgcg     >  1:427368/1‑150 (MQ=255)
                                                                                                                                                  ttCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCTCGTcgcg     <  2:2888834/150‑1 (MQ=255)
                                                                                                                                                    cGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCTCGTCGCGcc   >  1:3279155/1‑150 (MQ=255)
                                                                                                                                                     gATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCTCGTCGCGCCa  >  1:4180500/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGATGGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCTCGTCGCGCCA  >  CP009273/4193429‑4193727

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: