Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,841,044 C→A 61.8% A519A (GCG→GCT hycC ← hydrogenase 3, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,841,0440CA61.8% 278.4 / inf 408A519A (GCG→GCThycChydrogenase 3, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (79/77);  new base A (124/128);  total (203/205)
Fisher's exact test for biased strand distribution p-value = 8.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 409 total aligned reads displayed.

CATTTGCGTAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTA  >  CP009273/2840895‑2841192
                                                                                                                                                     |                                                                                                                                                    
 aTTTGCGTAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCCTTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCg                                                                                                                                                     >  2:4512662/1‑150 (MQ=255)
 aTTTGCGTAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCg                                                                                                                                                     <  2:4035944/150‑1 (MQ=255)
    tGCGTAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCa                                                                                                                                                  <  1:229408/150‑1 (MQ=255)
     gCGTAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCAt                                                                                                                                                 <  2:4864482/150‑1 (MQ=255)
      cGTAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGTGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATa                                                                                                                                                <  2:4450304/150‑1 (MQ=255)
      cGTAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGTGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATa                                                                                                                                                <  2:1999070/150‑1 (MQ=255)
      cGTAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATa                                                                                                                                                >  2:3944494/1‑150 (MQ=255)
            ttCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAatgatg                                                                                                                                          <  1:382554/150‑1 (MQ=255)
            ttCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAatgatg                                                                                                                                          <  1:564787/150‑1 (MQ=255)
             tCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAAtgatga                                                                                                                                         <  1:3274403/150‑1 (MQ=255)
                gCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATg                                                                                                                                      <  1:2677567/150‑1 (MQ=255)
                  aCCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGgc                                                                                                                                    >  2:2267442/1‑150 (MQ=255)
                  aCCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGgc                                                                                                                                    >  1:4289469/1‑150 (MQ=255)
                  aCCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGgc                                                                                                                                    >  2:885559/1‑150 (MQ=255)
                    cGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGgcag                                                                                                                                  <  2:3580704/150‑1 (MQ=255)
                       cgcAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGcagcag                                                                                                                               >  2:303007/1‑150 (MQ=255)
                         cAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCg                                                                                                                             <  2:174274/150‑1 (MQ=255)
                           aaCGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCggg                                                                                                                           >  1:4549359/1‑150 (MQ=255)
                           aaCGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCggg                                                                                                                           <  1:4826307/150‑1 (MQ=255)
                           aaCGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCggg                                                                                                                           >  1:1495189/1‑150 (MQ=255)
                            aCGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGggc                                                                                                                          <  2:2827659/150‑1 (MQ=255)
                             cGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGggca                                                                                                                         <  2:1067635/150‑1 (MQ=255)
                             cGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGggca                                                                                                                         <  2:1225462/150‑1 (MQ=255)
                                    ttCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATc                                                                                                                  <  2:3817257/150‑1 (MQ=255)
                                       aCCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATcagc                                                                                                               <  2:3992324/150‑1 (MQ=255)
                                       aCCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATcagc                                                                                                               <  2:2672369/150‑1 (MQ=255)
                                        ccGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATcagca                                                                                                              <  2:1617321/150‑1 (MQ=255)
                                         cGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATcagcag                                                                                                             >  1:1981459/1‑150 (MQ=255)
                                         cGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATcagcag                                                                                                             <  1:418914/150‑1 (MQ=255)
                                             aTGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAAc                                                                                                         <  1:1000478/150‑1 (MQ=255)
                                             aTGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAAc                                                                                                         >  1:4105478/1‑150 (MQ=255)
                                                 cAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGa                                                                                                     >  2:1898471/1‑150 (MQ=255)
                                                  aaaaCCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGat                                                                                                    <  2:727070/150‑1 (MQ=255)
                                                   aaaCCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGatc                                                                                                   <  1:152648/150‑1 (MQ=255)
                                                     aCCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGatcat                                                                                                 >  1:189033/1‑150 (MQ=255)
                                                     aCCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCACAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGatcat                                                                                                 <  2:2423318/150‑1 (MQ=255)
                                                        gTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCgg                                                                                              <  2:4194175/150‑1 (MQ=255)
                                                         tGAGCGGTAATCACCATTGATTTTACGTGGTCGTAACCGCACACCCAGGCCGCACAGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCgga                                                                                             >  1:2814855/1‑149 (MQ=255)
                                                           aGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTg                                                                                           <  2:4629895/150‑1 (MQ=255)
                                                            gCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTga                                                                                          <  1:505712/150‑1 (MQ=255)
                                                             cGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTgag                                                                                         <  1:1745488/150‑1 (MQ=255)
                                                               gTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGaa                                                                                       <  1:337845/150‑1 (MQ=255)
                                                               gTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGaa                                                                                       >  1:4599064/1‑150 (MQ=255)
                                                                tAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGaaa                                                                                      >  2:3723535/1‑150 (MQ=255)
                                                                 aaTCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAAc                                                                                     <  2:2335079/150‑1 (MQ=255)
                                                                  aTCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACg                                                                                    >  2:1174100/1‑150 (MQ=255)
                                                                       cATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCAACGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGATTGGCAGCAGAGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACggtggt                                                                               >  2:1395810/1‑150 (MQ=255)
                                                                          tGATTTTTCGTGGTCGTAACCGCACACCCCGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGtt                                                                            <  2:4424693/150‑1 (MQ=255)
                                                                          tGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGtt                                                                            <  2:950951/150‑1 (MQ=255)
                                                                           gATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTa                                                                           >  2:4101590/1‑150 (MQ=255)
                                                                           gATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTa                                                                           >  1:89477/1‑150 (MQ=255)
                                                                           gATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTa                                                                           >  1:4538476/1‑150 (MQ=255)
                                                                            aTTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAg                                                                          >  2:4851556/1‑150 (MQ=255)
                                                                             tttttCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGc                                                                         <  1:2377498/150‑1 (MQ=255)
                                                                             tttttCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGc                                                                         <  1:93268/150‑1 (MQ=255)
                                                                                ttCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAgg                                                                      >  1:1220168/1‑150 (MQ=255)
                                                                                ttCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAgg                                                                      >  2:4723303/1‑150 (MQ=255)
                                                                                ttCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAgg                                                                      <  2:2920395/150‑1 (MQ=255)
                                                                                   gTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTc                                                                   <  1:4697214/150‑1 (MQ=255)
                                                                                   gTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTc                                                                   <  1:3519740/150‑1 (MQ=255)
                                                                                    tGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTcc                                                                  <  1:4272465/150‑1 (MQ=255)
                                                                                     ggTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCa                                                                 <  2:1471595/150‑1 (MQ=255)
                                                                                        cGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCg                                                              >  1:1923065/1‑150 (MQ=255)
                                                                                              cGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCagag                                                        <  2:3054393/150‑1 (MQ=255)
                                                                                               gCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAgg                                                       >  2:1599736/1‑150 (MQ=255)
                                                                                               gCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAgg                                                       >  2:4614630/1‑150 (MQ=255)
                                                                                                 acacCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTa                                                     >  1:583601/1‑150 (MQ=255)
                                                                                                   acCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAca                                                   <  2:3445735/150‑1 (MQ=255)
                                                                                                            gCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCagag                                          <  2:2614352/150‑1 (MQ=255)
                                                                                                              aCCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCagagag                                        <  1:4044161/150‑1 (MQ=255)
                                                                                                                 gcgGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCAt                                     <  1:4021943/150‑1 (MQ=255)
                                                                                                                 gcgGGAACGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCAt                                     <  2:2389763/150‑1 (MQ=255)
                                                                                                                     gAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGc                                 <  2:2691837/150‑1 (MQ=255)
                                                                                                                      aaCGCGATGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCa                                >  1:3638625/1‑150 (MQ=255)
                                                                                                                            aTGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGcc                          >  1:105151/1‑150 (MQ=255)
                                                                                                                            aTGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGcc                          >  1:2473235/1‑150 (MQ=255)
                                                                                                                             tGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCa                         >  2:1960288/1‑150 (MQ=255)
                                                                                                                             tGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCa                         <  1:303164/150‑1 (MQ=255)
                                                                                                                             tGGCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCa                         >  2:676832/1‑150 (MQ=255)
                                                                                                                              ggCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCAc                        >  1:4146682/1‑150 (MQ=255)
                                                                                                                              ggCAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCAc                        <  1:209015/150‑1 (MQ=255)
                                                                                                                                cAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACgg                      >  1:1571215/1‑150 (MQ=255)
                                                                                                                                cAAACGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGAGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACgg                      >  2:1811008/1‑150 (MQ=255)
                                                                                                                                 aaaCGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGgc                     >  2:1738510/1‑150 (MQ=255)
                                                                                                                                 aaaCGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGgc                     >  2:2890729/1‑150 (MQ=255)
                                                                                                                                 aaaCGATCGCCTTTGCAAAGAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGgc                     <  2:4341598/150‑1 (MQ=255)
                                                                                                                                    cGATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCgca                  >  1:3735485/1‑150 (MQ=255)
                                                                                                                                     gATCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCgcag                 <  1:4457862/150‑1 (MQ=255)
                                                                                                                                      aTCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCgcagc                >  2:4392898/1‑150 (MQ=255)
                                                                                                                                       tCGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCgcagca               <  1:1146161/150‑1 (MQ=255)
                                                                                                                                        cGCCTTTGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCaa              >  2:4064184/1‑150 (MQ=255)
                                                                                                                                           cTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAacac           <  2:3245056/150‑1 (MQ=255)
                                                                                                                                            tttGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACAcc          <  1:4715680/150‑1 (MQ=255)
                                                                                                                                              tGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGc        >  1:2982877/1‑150 (MQ=255)
                                                                                                                                              tGCAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGc        >  2:1345583/1‑150 (MQ=255)
                                                                                                                                                cAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCa      <  2:2173723/150‑1 (MQ=255)
                                                                                                                                                cAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCa      <  1:1410778/150‑1 (MQ=255)
                                                                                                                                                cAAATAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCa      >  1:3786271/1‑150 (MQ=255)
                                                                                                                                                 aaaTAGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCaa     >  1:4575286/1‑150 (MQ=255)
                                                                                                                                                  aaTCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAAt    >  1:680596/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CATTTGCGTAGTTTCAGCACCGGCGCAAACGCCTGTTTCACCGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTA  >  CP009273/2840895‑2841192

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: