Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,559,507 T→A 32.0% 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,5070TA32.0% 211.5 / 147.3 150N162Y (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 (47/55);  new base A (25/23);  total (72/78)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

TCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTC  >  CP009273/1559362‑1559654
                                                                                                                                                 |                                                                                                                                                   
tCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATtgc                                                                                                                                                   <  1:3222096/148‑1 (MQ=255)
 ccTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATtgca                                                                                                                                                  <  1:3058964/148‑1 (MQ=255)
 ccTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATagca                                                                                                                                                  >  2:3303755/1‑148 (MQ=255)
   tGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATtgcact                                                                                                                                                >  2:1031538/1‑148 (MQ=255)
         cgggcggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACaa                                                                                                                                          <  1:1634853/148‑1 (MQ=255)
          gggcggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAat                                                                                                                                         <  2:2870505/148‑1 (MQ=255)
            gcggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCTCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAatat                                                                                                                                       <  2:2379757/148‑1 (MQ=255)
               gTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGc                                                                                                                                    >  2:1931503/1‑148 (MQ=255)
                         gCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATc                                                                                                                          <  1:1845685/148‑1 (MQ=255)
                          cAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCTGATCg                                                                                                                         >  2:2366490/1‑148 (MQ=255)
                           aaCTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGa                                                                                                                        <  2:561242/148‑1 (MQ=255)
                            aCTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAg                                                                                                                       >  2:2577502/1‑148 (MQ=255)
                            aCTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAg                                                                                                                       >  1:3207630/1‑148 (MQ=255)
                             cTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAgg                                                                                                                      >  2:1252681/1‑148 (MQ=255)
                              tgttgtAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGa                                                                                                                     >  2:722125/1‑148 (MQ=255)
                                ttgtACACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACa                                                                                                                   <  2:1725015/148‑1 (MQ=255)
                                   tAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATc                                                                                                                <  1:1956335/148‑1 (MQ=255)
                                   tAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATc                                                                                                                >  2:2253620/1‑148 (MQ=255)
                                   tAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATc                                                                                                                >  2:2757190/1‑148 (MQ=255)
                                     aaCGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCAc                                                                                                              >  2:3142256/1‑148 (MQ=255)
                                      aCGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCAcc                                                                                                             >  1:2162233/1‑148 (MQ=255)
                                      aCGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCAcc                                                                                                             >  1:453218/1‑148 (MQ=255)
                                      aCGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCAcc                                                                                                             <  1:2293116/148‑1 (MQ=255)
                                       cGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcg                                                                                                            >  2:2399599/1‑148 (MQ=255)
                                         aaTGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggt                                                                                                          >  1:2730057/1‑148 (MQ=255)
                                         aaTGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggt                                                                                                          >  1:2729822/1‑148 (MQ=255)
                                         aaTGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggt                                                                                                          <  1:762979/148‑1 (MQ=255)
                                          aTGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtc                                                                                                         >  1:3008842/1‑148 (MQ=255)
                                           tGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtcg                                                                                                        >  1:1848850/1‑148 (MQ=255)
                                           tGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtcg                                                                                                        >  1:3251184/1‑148 (MQ=255)
                                            gCGTTTATCGTCAGGGAATTCAGGATTGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCTGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtcgg                                                                                                       <  2:2971190/148‑1 (MQ=255)
                                             cGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCTCTGTCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtcggt                                                                                                      <  1:1185097/148‑1 (MQ=255)
                                             cGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtcggt                                                                                                      >  1:1409171/1‑148 (MQ=255)
                                               gTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGgtcggtcg                                                                                                    <  2:2555629/148‑1 (MQ=255)
                                                  aTCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCgaga                                                                                                 >  2:1489937/1‑148 (MQ=255)
                                                     ggCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATgg                                                                                              <  2:3153799/148‑1 (MQ=255)
                                                     ggCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATgg                                                                                              <  1:3214719/148‑1 (MQ=255)
                                                      gCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGt                                                                                             <  1:2198565/148‑1 (MQ=255)
                                                      gCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGt                                                                                             <  2:3267440/148‑1 (MQ=255)
                                                      gCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGt                                                                                             <  2:1992186/148‑1 (MQ=255)
                                                      gCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGt                                                                                             <  2:1630001/148‑1 (MQ=255)
                                                        aGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCa                                                                                           <  1:338709/148‑1 (MQ=255)
                                                          ggAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAAc                                                                                         >  1:2462408/1‑148 (MQ=255)
                                                                 aGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCAcc                                                                                                             >  1:1295968/1‑121 (MQ=255)
                                                                 aGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCAcc                                                                                                             <  2:1295968/121‑1 (MQ=255)
                                                                 aGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCaata                                                                                  <  1:1046305/148‑1 (MQ=255)
                                                                 aGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCaata                                                                                  <  1:1046190/148‑1 (MQ=255)
                                                                   gAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCaataat                                                                                >  2:484283/1‑148 (MQ=255)
                                                                    aaTGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATc                                                                               <  2:861235/148‑1 (MQ=255)
                                                                            ggAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAAtt                                                                           >  1:3297453/1‑144 (MQ=255)
                                                                            ggAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAAtt                                                                           >  1:3298185/1‑144 (MQ=255)
                                                                            ggAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAAtt                                                                           <  2:3298185/144‑1 (MQ=255)
                                                                             gAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCgg                                                                      >  1:2767030/1‑148 (MQ=255)
                                                                             gAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCgg                                                                      <  1:291504/148‑1 (MQ=255)
                                                                             gAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCgg                                                                      >  2:62704/1‑148 (MQ=255)
                                                                             gAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCgg                                                                      >  1:367080/1‑148 (MQ=255)
                                                                              aGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCggg                                                                     >  2:2134182/1‑148 (MQ=255)
                                                                                  tATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTg                                                                 <  2:679595/148‑1 (MQ=255)
                                                                                   aTCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGt                                                                <  1:1015244/148‑1 (MQ=255)
                                                                                   aTCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGt                                                                <  1:1015383/148‑1 (MQ=255)
                                                                                              tACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTgcgc                                                     >  1:1620485/1‑148 (MQ=255)
                                                                                               aCCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGcc                                                    >  2:1059150/1‑148 (MQ=255)
                                                                                                ccGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCa                                                   <  2:227903/148‑1 (MQ=255)
                                                                                                  gCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCAtt                                                 <  1:1767707/148‑1 (MQ=255)
                                                                                                   cTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCAttt                                                <  1:554794/148‑1 (MQ=255)
                                                                                                      gCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTa                                             <  1:511933/148‑1 (MQ=255)
                                                                                                      gCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTa                                             <  1:889460/148‑1 (MQ=255)
                                                                                                        ttCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACg                                           <  2:150674/148‑1 (MQ=255)
                                                                                                              aaTGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCa                                     >  1:2291672/1‑148 (MQ=255)
                                                                                                               aTGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCAt                                    >  2:1197463/1‑148 (MQ=255)
                                                                                                               aTGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCAt                                    >  2:1639999/1‑148 (MQ=255)
                                                                                                               aTGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCAt                                    >  1:2594991/1‑148 (MQ=255)
                                                                                                                  cAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATccc                                 >  2:3050039/1‑148 (MQ=255)
                                                                                                                     tAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTGTTTGCGCCCTTTCTACGCTGGCATCCCGTa                              >  2:1747820/1‑148 (MQ=255)
                                                                                                                     tAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTa                              >  2:1178504/1‑148 (MQ=255)
                                                                                                                      aaCCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTAc                             <  2:2607671/148‑1 (MQ=255)
                                                                                                                      aaCCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTAc                             <  2:2607675/148‑1 (MQ=255)
                                                                                                                      aaCCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGNCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTAc                             <  1:1682627/148‑1 (MQ=255)
                                                                                                                          aaaCATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAgc                         >  2:2193459/1‑148 (MQ=255)
                                                                                                                            aCATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAgcgc                       <  2:426031/148‑1 (MQ=255)
                                                                                                                            aCATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAgcgc                       <  2:2091186/148‑1 (MQ=255)
                                                                                                                             cATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCAcc                      >  2:2504400/1‑148 (MQ=255)
                                                                                                                               tttCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAg                    >  2:1290268/1‑148 (MQ=255)
                                                                                                                                ttCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAgg                   >  2:2955638/1‑148 (MQ=255)
                                                                                                                                ttCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAgg                   >  1:1168388/1‑148 (MQ=255)
                                                                                                                                 tCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGt                  <  1:1767312/148‑1 (MQ=255)
                                                                                                                                    gTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGtaat               <  1:2948437/148‑1 (MQ=255)
                                                                                                                                      aaaTGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGtaataa             >  2:2716466/1‑148 (MQ=255)
                                                                                                                                        aTGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAAc           >  2:1490881/1‑148 (MQ=255)
                                                                                                                                        aTGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAAc           <  1:71071/148‑1 (MQ=255)
                                                                                                                                        aTGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAAc           <  1:1352800/148‑1 (MQ=255)
                                                                                                                                         tGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTCATAAACt          <  1:1112873/148‑1 (MQ=255)
                                                                                                                                         tGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTTTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACt          <  2:429533/148‑1 (MQ=255)
                                                                                                                                          gcgcGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCTGTCGAGATGGTCAATGGCAANAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACtt         >  1:1681228/1‑148 (MQ=255)
                                                                                                                                             cgATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTcc      >  1:373132/1‑148 (MQ=255)
                                                                                                                                               aTTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTcccc    >  2:3182158/1‑148 (MQ=255)
                                                                                                                                               aTTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTcccc    <  1:62704/148‑1 (MQ=255)
                                                                                                                                               aTTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTcccc    <  1:62705/148‑1 (MQ=255)
                                                                                                                                                ttGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCt   <  2:2651267/148‑1 (MQ=255)
                                                                                                                                                 tgcactgcacAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTc  >  1:1571549/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTC  >  CP009273/1559362‑1559654

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: