Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,265,064 C→A 90.3% intergenic (‑95/‑335) yjcB ← / → yjcC putative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,265,0640CA90.3% 650.1 / 61.2 226intergenic (‑95/‑335)yjcB/yjcCputative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (9/13);  new base A (85/119);  total (94/132)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 234 total aligned reads displayed.

GCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTC  >  CP009273/4264917‑4265212
                                                                                                                                                   |                                                                                                                                                    
gCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAaaac                                                                                                                                                    >  2:3196847/1‑150 (MQ=255)
 cACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAaaca                                                                                                                                                   >  2:544002/1‑150 (MQ=255)
 cACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAAGTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCGGGCTATCGCGGGGCTTGCAGAAaaca                                                                                                                                                   <  2:1034620/150‑1 (MQ=255)
    tAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACaaaa                                                                                                                                                >  1:3327066/1‑150 (MQ=255)
    tAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACaaaa                                                                                                                                                <  2:4031995/150‑1 (MQ=255)
      agaagTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAAt                                                                                                                                              >  1:1923892/1‑150 (MQ=255)
      agaagTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAAt                                                                                                                                              <  2:3144207/150‑1 (MQ=255)
       gaagTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATg                                                                                                                                             <  2:821146/150‑1 (MQ=255)
       gaagTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGGAGCCATCAAACAAGTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATg                                                                                                                                             <  2:1648541/150‑1 (MQ=255)
        aagTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGa                                                                                                                                            <  2:4540302/150‑1 (MQ=255)
          gTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGaaa                                                                                                                                          <  2:4125344/150‑1 (MQ=255)
          gTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGaaa                                                                                                                                          <  2:2969615/150‑1 (MQ=255)
          gTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGaaa                                                                                                                                          >  2:4167086/1‑150 (MQ=255)
           ttGCCAGCTAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAac                                                                                                                                         <  1:812552/150‑1 (MQ=255)
           ttGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAac                                                                                                                                         <  1:1822876/150‑1 (MQ=255)
             gCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAacac                                                                                                                                       >  1:3970272/1‑150 (MQ=255)
              ccAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACAct                                                                                                                                      >  1:716584/1‑150 (MQ=255)
               cAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACActc                                                                                                                                     >  2:4019745/1‑150 (MQ=255)
                aGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACActct                                                                                                                                    >  1:2786473/1‑150 (MQ=255)
                    aagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCtgttt                                                                                                                                <  1:1285872/150‑1 (MQ=255)
                    aagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCtgttt                                                                                                                                <  2:2514321/150‑1 (MQ=255)
                    aagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCtgttt                                                                                                                                >  2:190756/1‑150 (MQ=255)
                    aagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCtgttg                                                                                                                                >  2:3391447/1‑149 (MQ=255)
                     agaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCtgtttg                                                                                                                               <  1:205984/150‑1 (MQ=255)
                      gaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCtgtttgt                                                                                                                              <  2:2574550/150‑1 (MQ=255)
                      gaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCtgtttgt                                                                                                                              >  2:3253899/1‑150 (MQ=255)
                      gaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCtgtttgt                                                                                                                              <  1:3523017/150‑1 (MQ=255)
                      gaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCtgtttgt                                                                                                                              <  1:955362/150‑1 (MQ=255)
                       aagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCtgtttgtt                                                                                                                             <  2:3043394/150‑1 (MQ=255)
                        agaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTgtttgttt                                                                                                                            <  1:48481/150‑1 (MQ=255)
                            ccACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTTCACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCAtt                                                                                                                        <  1:2385060/150‑1 (MQ=255)
                             cACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTa                                                                                                                       <  1:1042076/150‑1 (MQ=255)
                             cACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTa                                                                                                                       <  2:1352252/150‑1 (MQ=255)
                             cACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTa                                                                                                                       <  2:4470918/150‑1 (MQ=255)
                              aCGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTaa                                                                                                                      <  2:1464038/150‑1 (MQ=255)
                               cGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAAt                                                                                                                     <  1:339529/150‑1 (MQ=255)
                                 ccAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAAttt                                                                                                                   <  2:1220952/150‑1 (MQ=255)
                                   aGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTtg                                                                                                                 <  2:43107/150‑1 (MQ=255)
                                     tggtgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTtgtg                                                                                                               <  1:115517/150‑1 (MQ=255)
                                        tgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAAc                                                                                                            >  2:2655629/1‑150 (MQ=255)
                                           gggTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACtat                                                                                                         <  1:1142100/150‑1 (MQ=255)
                                             gTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTatat                                                                                                       <  1:4329724/150‑1 (MQ=255)
                                                 ccATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATcaca                                                                                                   >  2:2537461/1‑150 (MQ=255)
                                                   aTCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAat                                                                                                 <  1:1171540/150‑1 (MQ=255)
                                                    tCAAACAAGTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAatt                                                                                                <  2:1763334/150‑1 (MQ=255)
                                                     cAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAattg                                                                                               >  2:4650535/1‑150 (MQ=255)
                                                         cAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAattgattg                                                                                           >  2:1306523/1‑150 (MQ=255)
                                                           aTTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGAttgtt                                                                                         >  1:1969563/1‑150 (MQ=255)
                                                            ttCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGAttgttt                                                                                        <  1:3002888/150‑1 (MQ=255)
                                                               ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATtgtttgtt                                                                                     >  1:907070/1‑150 (MQ=255)
                                                               ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATtgtttgtt                                                                                     >  1:908082/1‑150 (MQ=255)
                                                                 ctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAg                                                                                   <  2:500120/150‑1 (MQ=255)
                                                                  tGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGc                                                                                  >  1:1809363/1‑150 (MQ=255)
                                                                      aaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCata                                                                              >  2:1476928/1‑150 (MQ=255)
                                                                        gagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCGGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATAtt                                                                            >  2:1769302/1‑150 (MQ=255)
                                                                        gagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATAtt                                                                            <  1:2091339/150‑1 (MQ=255)
                                                                        gagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATAtt                                                                            <  1:2905944/150‑1 (MQ=255)
                                                                        gagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATAtt                                                                            >  1:2475440/1‑150 (MQ=255)
                                                                         agCCGATGCCCTCGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTa                                                                           <  2:789068/150‑1 (MQ=255)
                                                                            cGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGc                                                                        >  1:3862659/1‑150 (MQ=255)
                                                                             gATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGcc                                                                       <  1:612797/150‑1 (MQ=255)
                                                                               tGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGt                                                                     >  1:4364120/1‑150 (MQ=255)
                                                                                   cTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACt                                                                 >  1:4628917/1‑150 (MQ=255)
                                                                                     ggACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACttt                                                               >  2:2832269/1‑150 (MQ=255)
                                                                                      gACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACtttt                                                              >  2:2679495/1‑150 (MQ=255)
                                                                                       aCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTa                                                             <  1:1306523/150‑1 (MQ=255)
                                                                                       aCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTa                                                             >  1:3700734/1‑150 (MQ=255)
                                                                                            gcggcgTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCt                                                        >  2:1132504/1‑150 (MQ=255)
                                                                                               gcgTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTa                                                     <  1:1467815/150‑1 (MQ=255)
                                                                                                 gttgttCACGTGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCTTTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACa                                                   <  2:4103787/150‑1 (MQ=255)
                                                                                                 gttgttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACa                                                   >  1:2167040/1‑150 (MQ=255)
                                                                                                  ttgttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAg                                                  <  1:2542633/150‑1 (MQ=255)
                                                                                                  ttgttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAg                                                  <  2:2446153/150‑1 (MQ=255)
                                                                                                  ttgttCACGGGAGGAGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAg                                                  <  2:592763/150‑1 (MQ=255)
                                                                                                     ttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGAtt                                               >  1:1352363/1‑150 (MQ=255)
                                                                                                     ttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGAtt                                               >  2:4745211/1‑150 (MQ=255)
                                                                                                       cACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAAGTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAt                                             >  1:3787834/1‑150 (MQ=255)
                                                                                                          gggAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgtg                                          <  2:2765506/150‑1 (MQ=255)
                                                                                                             aGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGtt                                       >  1:4164230/1‑150 (MQ=255)
                                                                                                             aGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGtt                                       <  1:1628076/150‑1 (MQ=255)
                                                                                                             aGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATGATGTGGtt                                       >  1:2918429/1‑150 (MQ=255)
                                                                                                                cGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGttttt                                    >  2:3595395/1‑150 (MQ=255)
                                                                                                                 gTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTc                                   <  2:3919433/150‑1 (MQ=255)
                                                                                                                 gTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTc                                   <  1:4060290/150‑1 (MQ=255)
                                                                                                                    gTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGt                                <  1:1126332/150‑1 (MQ=255)
                                                                                                                     tACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTg                               <  1:2102195/150‑1 (MQ=255)
                                                                                                                          tctGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCAtt                          <  2:2935568/150‑1 (MQ=255)
                                                                                                                            tGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTaa                        >  2:3792439/1‑150 (MQ=255)
                                                                                                                             ggCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCTGTGGCATTAAg                       >  1:4475177/1‑150 (MQ=255)
                                                                                                                               cTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAggg                     <  2:4626521/150‑1 (MQ=255)
                                                                                                                                tATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                    >  1:3479119/1‑150 (MQ=255)
                                                                                                                                 aTCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCa                   >  2:3120441/1‑150 (MQ=255)
                                                                                                                                 aTCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCa                   >  2:3120676/1‑150 (MQ=255)
                                                                                                                                 aTCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCa                   <  1:4737636/150‑1 (MQ=255)
                                                                                                                                 aTCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCa                   <  1:166720/150‑1 (MQ=255)
                                                                                                                                    gcgGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACTGATTATGTGGTTTTTCAGTGGCATTAAGGGCatat                <  1:190756/150‑1 (MQ=255)
                                                                                                                                    gcgGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCatat                <  2:2397068/150‑1 (MQ=255)
                                                                                                                                            tGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCccgc        >  2:2201123/1‑150 (MQ=255)
                                                                                                                                              cAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCccgccg      >  1:2605130/1‑150 (MQ=255)
                                                                                                                                                  acacAAAAATTAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCctc  <  1:162901/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTC  >  CP009273/4264917‑4265212

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: