Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,605,105 C→T 44.8% intergenic (‑168/+59) uxaB ← / ← yneF altronate oxidoreductase, NAD‑dependent/putative membrane‑bound diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,605,1050CT44.8% 116.6 / inf 239intergenic (‑168/+59)uxaB/yneFaltronate oxidoreductase, NAD‑dependent/putative membrane‑bound diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (73/59);  new base T (48/59);  total (121/118)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 242 total aligned reads displayed.

TAGTCATACAACCTGTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAGACACGTATCGGCTTCGACCAT  >  CP009273/1604956‑1605253
                                                                                                                                                     |                                                                                                                                                    
tAGTCATACAACCTGTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTc                                                                                                                                                      <  1:3700069/150‑1 (MQ=255)
 aGTCATACAACCTGTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTcc                                                                                                                                                     <  2:4592662/150‑1 (MQ=255)
  gTCATACAACCTGTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCt                                                                                                                                                    <  1:3249702/150‑1 (MQ=255)
       aCAACCTGTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCATTTTCTGGATc                                                                                                                                               <  2:199197/150‑1 (MQ=255)
        cAACCTGTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCa                                                                                                                                              <  1:2769684/150‑1 (MQ=255)
        cAACCTGTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCa                                                                                                                                              <  1:1057675/150‑1 (MQ=255)
         aaCCTGTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCAt                                                                                                                                             >  1:1426395/1‑150 (MQ=255)
         aaCCTGTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCAt                                                                                                                                             <  1:2557287/150‑1 (MQ=255)
            cTGTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATaaa                                                                                                                                          <  1:96985/150‑1 (MQ=255)
              gTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGt                                                                                                                                        <  2:2786707/150‑1 (MQ=255)
              gTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGt                                                                                                                                        <  2:3391519/150‑1 (MQ=255)
               tttGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTg                                                                                                                                       <  1:899172/150‑1 (MQ=255)
                  gAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTa                                                                                                                                    <  1:1375025/150‑1 (MQ=255)
                    aTTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGa                                                                                                                                  <  2:4527221/150‑1 (MQ=255)
                       ggTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAaca                                                                                                                               <  1:4395202/150‑1 (MQ=255)
                          aCGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCAAGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACAt                                                                                                                            >  2:375953/1‑150 (MQ=255)
                           cGACAGGTTAGCAAGCTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACAtt                                                                                                                           >  2:206228/1‑150 (MQ=255)
                           cGACAGGTTAGCAAGCTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACAtt                                                                                                                           >  2:1611185/1‑150 (MQ=255)
                           cGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACAtt                                                                                                                           <  1:303853/150‑1 (MQ=255)
                            gACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTg                                                                                                                          <  1:1393238/150‑1 (MQ=255)
                            gACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTg                                                                                                                          >  2:747651/1‑150 (MQ=255)
                               aGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCAt                                                                                                                       >  1:843716/1‑150 (MQ=255)
                               aGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCAt                                                                                                                       >  2:2190103/1‑150 (MQ=255)
                                 gTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTc                                                                                                                     >  2:3137045/1‑150 (MQ=255)
                                 gTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTc                                                                                                                     >  1:1771104/1‑150 (MQ=255)
                                  ttAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCa                                                                                                                    <  2:2797058/150‑1 (MQ=255)
                                  ttAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATAATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCa                                                                                                                    <  1:345712/150‑1 (MQ=255)
                                     gCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAAt                                                                                                                 <  2:2679637/150‑1 (MQ=255)
                                          cTTTAATACGCCGAACCCCTGTTTTGATCAGCTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATcgcgc                                                                                                            >  1:3889189/1‑150 (MQ=255)
                                          cTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATcgcgc                                                                                                            <  1:3963175/150‑1 (MQ=255)
                                           tttAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATcgcgcg                                                                                                           <  1:1634619/150‑1 (MQ=255)
                                            ttAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCAGTCAAATCGCGCGt                                                                                                          >  2:4333973/1‑150 (MQ=255)
                                             tAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTa                                                                                                         >  1:1390486/1‑150 (MQ=255)
                                              aaTACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGAGCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTaa                                                                                                        >  2:1523521/1‑150 (MQ=255)
                                                tACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATg                                                                                                      <  2:361735/150‑1 (MQ=255)
                                                tACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATg                                                                                                      >  1:1086353/1‑150 (MQ=255)
                                                 aCGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGa                                                                                                     <  2:4351848/150‑1 (MQ=255)
                                                    ccGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATa                                                                                                  <  1:4089636/150‑1 (MQ=255)
                                                    ccGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATa                                                                                                  <  2:3090031/150‑1 (MQ=255)
                                                    ccGAACACCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATa                                                                                                  <  2:1595530/150‑1 (MQ=255)
                                                     cGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATaa                                                                                                 >  1:2218623/1‑150 (MQ=255)
                                                      gAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATaaa                                                                                                >  2:1501055/1‑150 (MQ=255)
                                                      gAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATaaa                                                                                                >  1:1681468/1‑150 (MQ=255)
                                                        aCCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGa                                                                                              <  2:802280/150‑1 (MQ=255)
                                                         ccccTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGAt                                                                                             >  1:2015946/1‑150 (MQ=255)
                                                          cccTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATg                                                                                            >  2:4257244/1‑150 (MQ=255)
                                                            cTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTc                                                                                          <  2:1850198/150‑1 (MQ=255)
                                                             tGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCa                                                                                         >  1:1802050/1‑150 (MQ=255)
                                                               ttttGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGa                                                                                       >  2:3497603/1‑150 (MQ=255)
                                                                tttGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGAc                                                                                      <  2:571458/150‑1 (MQ=255)
                                                                  tGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACaa                                                                                    <  1:716345/150‑1 (MQ=255)
                                                                   gATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAAc                                                                                   <  2:1096071/150‑1 (MQ=255)
                                                                     tCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACtt                                                                                 >  2:2353140/1‑150 (MQ=255)
                                                                       aaCTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTcc                                                                               >  2:1036965/1‑150 (MQ=255)
                                                                        aCTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCt                                                                              >  1:2864333/1‑150 (MQ=255)
                                                                        aCTCCTGATGAATAATTAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCt                                                                              <  2:2962855/150‑1 (MQ=255)
                                                                         cTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTTCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTc                                                                             >  2:3541619/1‑150 (MQ=255)
                                                                           ccTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCAc                                                                           >  1:2842180/1‑150 (MQ=255)
                                                                            cTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCAcc                                                                          <  2:1230155/150‑1 (MQ=255)
                                                                                tgatTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTaa                                                                      <  1:4532266/150‑1 (MQ=255)
                                                                                  atTAATGAGCAGTTTTTTGATTAAAGTGTGGCGCGGATCATTGTTTAATCGAGTAAAAACGCCTTTTCCTGGATCTTAAAGTGGTAGAACACATTGCATTCAAATCTCGCTTTATGAATAAAGATGTCATACAACTTCCTCACCGTAACg                                                                    <  2:2105637/150‑1 (MQ=255)
                                                                                   ttAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGc                                                                   >  1:4494490/1‑150 (MQ=255)
                                                                                    tAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTACCTCACCGTAACGCa                                                                  >  1:524770/1‑150 (MQ=255)
                                                                                     aaTGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCAt                                                                 >  1:3012425/1‑150 (MQ=255)
                                                                                      aTGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATa                                                                >  1:1525089/1‑150 (MQ=255)
                                                                                      aTGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATa                                                                <  1:4370899/150‑1 (MQ=255)
                                                                                      aTGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATa                                                                <  2:4165624/150‑1 (MQ=255)
                                                                                       tGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAg                                                               <  2:483072/150‑1 (MQ=255)
                                                                                       tGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAg                                                               <  1:204833/150‑1 (MQ=255)
                                                                                        gAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGt                                                              <  2:3795878/150‑1 (MQ=255)
                                                                                         aGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTg                                                             <  2:4330773/150‑1 (MQ=255)
                                                                                          acAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCATTTTTCTGGATCATAAAGTGGTATAACACATTGAATTCAAATAGCGCGTAATGAATAAAGATGTCAGACAACTTACACACCGTAACGCATAGAGc                                                            >  2:3927383/2‑150 (MQ=255)
                                                                                           cAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCt                                                           <  1:2044791/150‑1 (MQ=255)
                                                                                                 tATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACg                                                     >  2:3552182/1‑150 (MQ=255)
                                                                                                    gagaAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGtt                                                  <  1:1121201/150‑1 (MQ=255)
                                                                                                     agaAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTg                                                 >  2:2646994/1‑150 (MQ=255)
                                                                                                      gaAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTgc                                                <  1:1902252/150‑1 (MQ=255)
                                                                                                       aaaaGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTgcg                                               <  1:568561/150‑1 (MQ=255)
                                                                                                         aaGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCGGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGcc                                             <  2:3846268/150‑1 (MQ=255)
                                                                                                          aGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGccc                                            >  1:4523512/1‑150 (MQ=255)
                                                                                                            tgtgGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCAt                                          >  2:243590/1‑150 (MQ=255)
                                                                                                                    ggATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTttcttc                                  >  2:362043/1‑150 (MQ=255)
                                                                                                                    ggATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTttcttc                                  >  2:950734/1‑150 (MQ=255)
                                                                                                                       tCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAg                               <  1:1441392/150‑1 (MQ=255)
                                                                                                                       tCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAg                               <  1:4584106/150‑1 (MQ=255)
                                                                                                                            gTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATa                          >  2:1454569/1‑150 (MQ=255)
                                                                                                                            gTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATa                          <  2:1356755/150‑1 (MQ=255)
                                                                                                                             tttAATCGAGGAAAAACGCCTTGTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATAc                         <  2:544723/150‑1 (MQ=255)
                                                                                                                               tAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAg                       <  1:2134647/150‑1 (MQ=255)
                                                                                                                                  tCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAGaca                    >  1:2423802/1‑150 (MQ=255)
                                                                                                                                  tCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAGaca                    >  1:4429442/1‑150 (MQ=255)
                                                                                                                                   cGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAGacac                   >  1:2813740/1‑150 (MQ=255)
                                                                                                                                     aGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTTCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATATTTCTTCGAGCGATACAGACACGt                 >  2:3288373/1‑150 (MQ=255)
                                                                                                                                      ggAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAGACACGTa                >  2:4423920/1‑150 (MQ=255)
                                                                                                                                      ggAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAGACACGTa                >  2:4423758/1‑150 (MQ=255)
                                                                                                                                          aaaCGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAGACACGTATCgg            >  2:2138740/1‑150 (MQ=255)
                                                                                                                                            aCGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAGACACGTATCGGCt          >  1:637876/1‑150 (MQ=255)
                                                                                                                                             cGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAGACACGTATCGGCtt         >  2:964332/1‑150 (MQ=255)
                                                                                                                                              gCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAGACACGTATCGGCTTc        <  2:1099612/150‑1 (MQ=255)
                                                                                                                                                 ttttCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAGACACGTATCGGCTTCGAc     >  2:3976638/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAGTCATACAACCTGTTTGAATTGGTACGACAGGTTAGCAAACTTTAATACGCCGAACCCCTGTTTTGATCAACTCCTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCACCGTAACGCATAGTGCTGGTACGGTTGCGCCCATCTTTCTTCGAGCGATACAGACACGTATCGGCTTCGACCAT  >  CP009273/1604956‑1605253

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: