Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,193,578 G→T 93.6% 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,5780GT93.6% 442.4 / 20.6 155D111Y (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 (4/6);  new base T (77/68);  total (81/74)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 156 total aligned reads displayed.

GGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCTCGTCG  >  CP009273/4193433‑4193722
                                                                                                                                                 |                                                                                                                                                
ggACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTAt                                                                                                                                                <  1:3132738/148‑1 (MQ=255)
 gACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCtata                                                                                                                                               >  2:576632/1‑148 (MQ=255)
    gtcATAGCCACGCTGTAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGGGTACTCCAGCGTCGTGACGGCGGTAGAGGCACTGAAAACTGGGGCTCTGGATTATGTCTTCAATCCGCTGGATTTCTATAAcc                                                                                                                                            <  1:2968095/146‑1 (MQ=255)
     gCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCt                                                                                                                                           <  1:1286940/148‑1 (MQ=255)
      cATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTa                                                                                                                                          <  2:847780/148‑1 (MQ=255)
       aTCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTAc                                                                                                                                         <  2:3065916/148‑1 (MQ=255)
        tCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACa                                                                                                                                        <  1:381627/148‑1 (MQ=255)
                  gAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTg                                                                                                                              >  2:2812248/1‑148 (MQ=255)
                   aaaGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTgg                                                                                                                             >  1:1287253/1‑148 (MQ=255)
                    aaGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACATACAGGCGACGCTGGa                                                                                                                            >  2:920354/1‑148 (MQ=255)
                     agagATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGaa                                                                                                                           >  2:1653143/1‑148 (MQ=255)
                     agagATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGaa                                                                                                                           >  2:1653050/1‑148 (MQ=255)
                     agagATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGaa                                                                                                                           >  2:1128899/1‑148 (MQ=255)
                       agaTCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGaaaa                                                                                                                         <  1:898398/148‑1 (MQ=255)
                       agaTCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGaaaa                                                                                                                         <  1:898363/148‑1 (MQ=255)
                       agaTCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGaaaa                                                                                                                         <  1:710624/148‑1 (MQ=255)
                       agaTCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGaaaa                                                                                                                         <  1:2371059/148‑1 (MQ=255)
                          tCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAgc                                                                                                                      >  1:2375900/1‑148 (MQ=255)
                           cAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAgcg                                                                                                                     <  2:2613719/148‑1 (MQ=255)
                             aaGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCt                                                                                                                   >  2:2956064/1‑148 (MQ=255)
                             aaGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCt                                                                                                                   >  2:433353/1‑148 (MQ=255)
                                   tAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCa                                                                                                             >  1:3258986/1‑148 (MQ=255)
                                   tAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCa                                                                                                             <  1:102733/148‑1 (MQ=255)
                                   tAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCa                                                                                                             <  1:102671/148‑1 (MQ=255)
                                    aaaCCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCAt                                                                                                            <  1:377847/148‑1 (MQ=255)
                                     aaCCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATa                                                                                                           <  1:2265940/148‑1 (MQ=255)
                                       cccGGCAATTCCGGTGCTTATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACg                                                                                                         <  1:295088/148‑1 (MQ=255)
                                       cccGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACg                                                                                                         <  2:1997069/148‑1 (MQ=255)
                                       cccGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACg                                                                                                         <  2:1997063/148‑1 (MQ=255)
                                       cccGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACg                                                                                                         >  1:1556058/1‑148 (MQ=255)
                                         cGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGca                                                                                                       <  1:576632/148‑1 (MQ=255)
                                          ggCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGcac                                                                                                      >  1:2809083/1‑148 (MQ=255)
                                           gCAATTCCGGTGCTGATTCTGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGcaca                                                                                                     <  2:2366880/148‑1 (MQ=255)
                                            cAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAg                                                                                                    >  2:1890254/1‑148 (MQ=255)
                                             aaTTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGt                                                                                                   <  1:2399409/148‑1 (MQ=255)
                                             aaTTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCGGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGt                                                                                                   <  1:1399702/148‑1 (MQ=255)
                                               ttCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTAt                                                                                                 >  2:2308531/1‑148 (MQ=255)
                                               ttCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTAt                                                                                                 <  2:2184863/148‑1 (MQ=255)
                                               ttCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTAt                                                                                                 >  2:1735823/1‑148 (MQ=255)
                                                 ccGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTg                                                                                               <  2:158138/148‑1 (MQ=255)
                                                 ccGGTGCTGATTATCACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTg                                                                                               <  1:563507/148‑1 (MQ=255)
                                                    gTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATg                                                                                            <  2:2213961/148‑1 (MQ=255)
                                                    gTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATg                                                                                            <  2:2728400/148‑1 (MQ=255)
                                                     tGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGc                                                                                           <  2:3171329/148‑1 (MQ=255)
                                                     tGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGc                                                                                           >  1:1281711/1‑148 (MQ=255)
                                                      gCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCt                                                                                          >  1:2793522/1‑148 (MQ=255)
                                                          aTTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGaaa                                                                                      >  2:808756/1‑148 (MQ=255)
                                                            tATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAaca                                                                                    <  2:2113192/148‑1 (MQ=255)
                                                                  tGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCg                                                                              >  1:289699/1‑148 (MQ=255)
                                                                  tGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCg                                                                              >  1:289994/1‑148 (MQ=255)
                                                                      tACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTCTTGATGCTGAAACACCTGCGGTGa                                                                          >  2:40570/1‑148 (MQ=255)
                                                                      tACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGa                                                                          <  2:2014226/148‑1 (MQ=255)
                                                                              cGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTgccagc                                                                  >  2:1879449/1‑148 (MQ=255)
                                                                                  gagaCGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGt                                                              <  1:726376/148‑1 (MQ=255)
                                                                                     aCGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCg                                                           <  1:2521119/148‑1 (MQ=255)
                                                                                     aCGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCg                                                           >  2:822782/1‑148 (MQ=255)
                                                                                     aCGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCg                                                           >  2:822623/1‑148 (MQ=255)
                                                                                       ggcggTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGt                                                         >  2:3216847/1‑148 (MQ=255)
                                                                                            tAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGt                                                    <  1:2155280/148‑1 (MQ=255)
                                                                                            tAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGt                                                    >  1:3021345/1‑148 (MQ=255)
                                                                                            tAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGt                                                    <  1:2155006/148‑1 (MQ=255)
                                                                                            tAGAGGAACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAAACTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACAACTGCGGTGACTGCCAGCCAGTTCGTGATGGt                                                    >  1:1374159/1‑148 (MQ=255)
                                                                                             agagGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTCTGGTc                                                   >  2:340673/1‑148 (MQ=255)
                                                                                             agagGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTc                                                   >  2:340642/1‑148 (MQ=255)
                                                                                             agagGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTc                                                   >  1:2889743/1‑148 (MQ=255)
                                                                                             agagGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTc                                                   >  2:379446/1‑148 (MQ=255)
                                                                                              gagGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCg                                                  <  1:1860241/148‑1 (MQ=255)
                                                                                               agGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCgg                                                 >  2:3161301/1‑148 (MQ=255)
                                                                                                   aCTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTaaa                                             <  1:2004618/148‑1 (MQ=255)
                                                                                                      gAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGc                                          <  2:1286896/148‑1 (MQ=255)
                                                                                                      gAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGc                                          <  2:1287253/148‑1 (MQ=255)
                                                                                                        aaaCTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGccc                                        >  1:2006361/1‑148 (MQ=255)
                                                                                                         aaCTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCg                                       >  2:2050172/1‑148 (MQ=255)
                                                                                                         aaCTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCg                                       >  2:2049201/1‑148 (MQ=255)
                                                                                                          aCTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGcc                                         >  1:134220/1‑145 (MQ=255)
                                                                                                            tGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCg                                    <  1:2812248/148‑1 (MQ=255)
                                                                                                              gggCGCTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGAt                                  >  1:2518737/1‑148 (MQ=255)
                                                                                                                 cgcTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCa                               <  1:1128899/148‑1 (MQ=255)
                                                                                                                 cgcTGGATTATCTAATCTAGCCGCTGGATTTCTATAACCTTCAGGCGACGCTGGAAAAAGCGCTCGAGCATACGCACAGTATTGATGCTGAAACACATGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCa                               <  2:2521747/148‑1 (MQ=255)
                                                                                                                  gcTGGATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCaa                              <  2:1947523/148‑1 (MQ=255)
                                                                                                                    tGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAaca                            <  2:1260055/148‑1 (MQ=255)
                                                                                                                      gATTATCTCATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACAcc                          >  2:2197216/1‑148 (MQ=255)
                                                                                                                       aTTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCt                         >  1:3110002/1‑148 (MQ=255)
                                                                                                                           tctcATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTc                     >  2:3290655/1‑148 (MQ=255)
                                                                                                                            ctcATCAAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCa                    >  2:1971349/1‑148 (MQ=255)
                                                                                                                             tcatcaAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAg                   >  2:294986/1‑148 (MQ=255)
                                                                                                                              catcaAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGt                  >  2:1558365/1‑148 (MQ=255)
                                                                                                                              catcaAGCCGCGGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGt                  >  2:61757/1‑148 (MQ=255)
                                                                                                                               atcaAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTg                 >  1:3230929/1‑148 (MQ=255)
                                                                                                                               atcaAGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTg                 >  1:2720608/1‑148 (MQ=255)
                                                                                                                                  aaGCCGCTGGATTTCTATAACCTACCGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGaaa              <  2:244137/148‑1 (MQ=255)
                                                                                                                                  aaGCCGCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGaaa              >  2:1149663/1‑148 (MQ=255)
                                                                                                                                       gCTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGcc         <  1:1803220/148‑1 (MQ=255)
                                                                                                                                        cTGGATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGccc        <  1:433353/148‑1 (MQ=255)
                                                                                                                                          ggATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtc      >  2:2513073/1‑148 (MQ=255)
                                                                                                                                          ggATTTCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtc      >  2:2513109/1‑148 (MQ=255)
                                                                                                                                             tttCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtcgtc   >  1:155883/1‑148 (MQ=255)
                                                                                                                                              ttCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtcgtcg  >  1:2512616/1‑148 (MQ=255)
                                                                                                                                              ttCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtcgtcg  <  1:3052101/148‑1 (MQ=255)
                                                                                                                                              ttCTATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCtcgtcg  >  1:2585648/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGACGGCATCGCCACGCTGAAAGAGATCAAAGCGTTAAACCCGGCAATTCCGGTGCTGATTATGACTGCGTACTCCAGCGTCGAGACGGCGGTAGAGGCACTGAAAACTGGGGCGCTGGATTATCTCATCAAGCCGCTGGATTTCGATAACCTACAGGCGACGCTGGAAAAAGCGCTCGCGCATACGCACAGTATTGATGCTGAAACACCTGCGGTGACTGCCAGCCAGTTCGGTATGGTCGGTAAAAGCCCGGCGATGCAACACCTGCTCAGTGAAATCGCCCTCGTCG  >  CP009273/4193433‑4193722

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: