Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,559,507 T→A 40.3% N162Y (AAT→TAT)  dosP ← oxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,559,5070TA40.3% 176.8 / 298.8 233N162Y (AAT→TAT) dosPoxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator
Reads supporting (aligned to +/- strand):  ref base T (70/69);  new base A (44/50);  total (114/119)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

TCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCC  >  CP009273/1559359‑1559652
                                                                                                                                                    |                                                                                                                                                 
   tCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATtgca                                                                                                                                                >  2:1834731/1‑149 (MQ=255)
        aTCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATagcactgca                                                                                                                                           <  1:585520/149‑1 (MQ=255)
        aTCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGAGtgcactgca                                                                                                                                           >  1:2681643/1‑149 (MQ=255)
             gggcggTTTTCCAGAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAata                                                                                                                                      <  1:2862581/149‑1 (MQ=255)
              ggcggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTCAATGCGCGATAGCACTGCACAatat                                                                                                                                     >  1:1675404/1‑149 (MQ=255)
              ggcggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAatat                                                                                                                                     <  1:2857393/149‑1 (MQ=255)
              ggcggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAatat                                                                                                                                     <  1:3750154/149‑1 (MQ=255)
                cggTTTTCCATAGCAACTGTTGTAAACGAAGGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATgg                                                                                                                                   <  1:1468849/149‑1 (MQ=255)
                 ggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGc                                                                                                                                  >  2:3530803/1‑149 (MQ=255)
                 ggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGc                                                                                                                                  >  2:3794851/1‑149 (MQ=255)
                 ggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGc                                                                                                                                  >  2:468165/1‑149 (MQ=255)
                 ggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGc                                                                                                                                  >  2:1752304/1‑149 (MQ=255)
                  gTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCg                                                                                                                                 >  2:413552/1‑149 (MQ=255)
                  gTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCg                                                                                                                                 >  1:1422374/1‑149 (MQ=255)
                  gTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCg                                                                                                                                 >  2:2434222/1‑149 (MQ=255)
                   ttttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGt                                                                                                                                >  1:607571/1‑149 (MQ=255)
                   ttttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGt                                                                                                                                >  2:2803855/1‑149 (MQ=255)
                   ttttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGt                                                                                                                                >  2:1434028/1‑149 (MQ=255)
                   ttttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGt                                                                                                                                >  2:4084515/1‑149 (MQ=255)
                      tCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTcc                                                                                                                             <  1:1182845/149‑1 (MQ=255)
                         aTAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGa                                                                                                                          <  2:1058924/149‑1 (MQ=255)
                             cAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGa                                                                                                                      <  1:3451456/149‑1 (MQ=255)
                              aaCTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAg                                                                                                                     <  2:2582229/149‑1 (MQ=255)
                              aaCTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAg                                                                                                                     >  1:2467957/1‑149 (MQ=255)
                               aCTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAgg                                                                                                                    >  1:511939/1‑149 (MQ=255)
                                cTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGa                                                                                                                   >  2:3528299/1‑149 (MQ=255)
                                 tgttgtAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGAc                                                                                                                  <  1:4025547/149‑1 (MQ=255)
                                 tgttgtAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGAc                                                                                                                  >  1:2795548/1‑149 (MQ=255)
                                     gtAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATc                                                                                                              <  2:657376/149‑1 (MQ=255)
                                     gtAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATc                                                                                                              <  2:2320326/149‑1 (MQ=255)
                                     gtAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATc                                                                                                              >  1:783614/1‑149 (MQ=255)
                                       aaaCGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCAc                                                                                                            >  1:2609768/1‑149 (MQ=255)
                                          cGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcgg                                                                                                         <  2:427715/149‑1 (MQ=255)
                                          cGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcgg                                                                                                         >  2:1370565/1‑149 (MQ=255)
                                           gAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggt                                                                                                        <  1:3158362/149‑1 (MQ=255)
                                           gAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggt                                                                                                        >  2:589707/1‑149 (MQ=255)
                                            aaTGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtc                                                                                                       <  2:3941053/149‑1 (MQ=255)
                                            aaTGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtc                                                                                                       >  1:3284061/1‑149 (MQ=255)
                                               gCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtcggt                                                                                                    >  1:1711728/1‑149 (MQ=255)
                                               gCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtcggt                                                                                                    <  1:919268/149‑1 (MQ=255)
                                                 ggTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGgtcggtcg                                                                                                  >  2:3598120/1‑149 (MQ=255)
                                                  gTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCga                                                                                                 >  2:1495237/1‑149 (MQ=255)
                                                    tATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCgaga                                                                                               <  1:2428180/149‑1 (MQ=255)
                                                          cAGGGAATTCAGGAAGGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCa                                                                                         <  2:93615/149‑1 (MQ=255)
                                                           aGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCaa                                                                                        <  2:83072/149‑1 (MQ=255)
                                                           aGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCaa                                                                                        <  2:192100/149‑1 (MQ=255)
                                                           aGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCaa                                                                                        <  1:2739609/149‑1 (MQ=255)
                                                           aGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCaa                                                                                        >  1:2989495/1‑149 (MQ=255)
                                                            gggAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGTGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAAc                                                                                       <  1:838400/149‑1 (MQ=255)
                                                            gggAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAAc                                                                                       <  2:524739/149‑1 (MQ=255)
                                                             ggAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACg                                                                                      >  2:1293832/1‑149 (MQ=255)
                                                               aaTTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGc                                                                                    <  2:2398714/149‑1 (MQ=255)
                                                                    aGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCaataa                                                                               >  2:295656/1‑149 (MQ=255)
                                                                     ggAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCaataat                                                                              <  2:3284061/149‑1 (MQ=255)
                                                                       aaTGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCa                                                                            >  2:1664169/1‑149 (MQ=255)
                                                                       aaTGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCa                                                                            >  2:3831700/1‑149 (MQ=255)
                                                                        aTGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCaa                                                                           <  1:762079/149‑1 (MQ=255)
                                                                        aTGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCaa                                                                           >  2:594742/1‑149 (MQ=255)
                                                                         tGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAAt                                                                          <  2:769389/149‑1 (MQ=255)
                                                                         tGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAAt                                                                          <  1:1134073/149‑1 (MQ=255)
                                                                            tCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGt                                                                       <  2:1408506/149‑1 (MQ=255)
                                                                             cAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTc                                                                      >  2:2784464/1‑149 (MQ=255)
                                                                               ggAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCgg                                                                    <  2:3828874/149‑1 (MQ=255)
                                                                                gAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCggg                                                                   <  1:1752304/149‑1 (MQ=255)
                                                                                gAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCggg                                                                   <  1:3794851/149‑1 (MQ=255)
                                                                                   tgtATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCt                                                                >  2:3424152/1‑149 (MQ=255)
                                                                                         gggCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCtt                                                          >  1:3314794/1‑149 (MQ=255)
                                                                                           gCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCtttt                                                        <  1:3357417/149‑1 (MQ=255)
                                                                                                aTACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTgcgc                                                   >  1:2034092/1‑149 (MQ=255)
                                                                                                aTACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAAGCAATTGTCGGGTCTGTTCTTTTTgcgc                                                   >  1:3271602/1‑149 (MQ=255)
                                                                                                   ccGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCAt                                                <  1:3666276/149‑1 (MQ=255)
                                                                                                     gCTGGCTTCTCTAATGCAGTAACCTAACATTTCGGTTAATGCGGGATAGGAGTGCACTTTTTGTCGTTGGGGATCGAGGACTATCACCGGTCGGTCTAGATGTTCAACGGCTATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCtttt                                              <  1:2907433/149‑1 (MQ=255)
                                                                                                     gCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCAttt                                              <  1:1959805/149‑1 (MQ=255)
                                                                                                     gCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCAttt                                              <  2:2391686/149‑1 (MQ=255)
                                                                                                      cTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTc                                             >  2:1215232/1‑149 (MQ=255)
                                                                                                         gCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTAc                                          >  1:1326805/1‑149 (MQ=255)
                                                                                                          cTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATACTCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACg                                         >  2:1593480/1‑149 (MQ=255)
                                                                                                                   tGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATcc                                <  2:2467957/149‑1 (MQ=255)
                                                                                                                    gCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATccc                               >  2:2656765/1‑149 (MQ=255)
                                                                                                                      aGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGt                             >  2:231364/1‑149 (MQ=255)
                                                                                                                            cAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAgc                       >  2:2832621/1‑149 (MQ=255)
                                                                                                                               aCATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGcc                    <  2:2604897/149‑1 (MQ=255)
                                                                                                                                 aTTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAg                  >  2:2706000/1‑149 (MQ=255)
                                                                                                                                 aTTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAg                  <  1:1834731/149‑1 (MQ=255)
                                                                                                                                 aTTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAg                  <  1:2970722/149‑1 (MQ=255)
                                                                                                                                 aTTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAg                  >  1:3684479/1‑149 (MQ=255)
                                                                                                                                 aTTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAg                  >  1:3122556/1‑149 (MQ=255)
                                                                                                                                  tttCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAgg                 <  1:1665229/149‑1 (MQ=255)
                                                                                                                                  tttCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGACGACAATCACGGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAgg                 >  1:4035190/1‑149 (MQ=255)
                                                                                                                                     cGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGtaa              <  2:2031671/149‑1 (MQ=255)
                                                                                                                                      ggTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGtaat             >  2:3035426/1‑149 (MQ=255)
                                                                                                                                      ggTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCTGTCGAGATGGTCAACGTCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACTCTGGCATCCCGTACCAGCGCCAGGtaat             <  1:2187513/149‑1 (MQ=255)
                                                                                                                                       ggaattGCTCGAGTGCACGGCAAATTATGGCGTTCCGGATCGAGGACAATCACAGGTCGGTAGAGATGGTCAAGGGCAATAATAAATTGTCGGGTCTGTTCTTTTTGCGACATTTGTACGCTGGCATCCCGTACCAGCGCCAGGtaata            <  2:2259836/144‑1 (MQ=255)
                                                                                                                                       gTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGtaata            >  1:1364633/1‑149 (MQ=255)
                                                                                                                                         aaaTGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATaaa          >  1:2508639/1‑149 (MQ=255)
                                                                                                                                           aTGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACt        >  1:1513392/1‑149 (MQ=255)
                                                                                                                                           aTGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACt        >  1:1900458/1‑149 (MQ=255)
                                                                                                                                           aTGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACt        >  1:1900891/1‑149 (MQ=255)
                                                                                                                                             gcgcGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACttt      <  1:413552/149‑1 (MQ=255)
                                                                                                                                                cgATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTccc   >  2:3559922/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCC  >  CP009273/1559359‑1559652

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: