Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,523,345 G→T 7.7% A86E (GCG→GAG)  yjhQ ← putative acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,523,3450GT 7.7% 969.6 / 37.3 260A86E (GCG→GAG) yjhQputative acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (112/128);  new base T (4/16);  total (116/144)
Fisher's exact test for biased strand distribution p-value = 3.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 272 total aligned reads displayed.

CTCCGGCACAGGGTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAAGAGGCC  >  CP009273/4523196‑4523494
                                                                                                                                                     |                                                                                                                                                     
ctcCGGCACAGGGTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCg                                                                                                                                                       <  2:4224346/150‑1 (MQ=255)
ctcCGGCACAGGGTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCg                                                                                                                                                       <  2:4223527/150‑1 (MQ=255)
  ccGGCACAGGGTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGcc                                                                                                                                                     <  2:3035013/150‑1 (MQ=255)
   cGGCACAGGGTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCa                                                                                                                                                    <  2:115756/150‑1 (MQ=255)
   cGGCACAGGGTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCa                                                                                                                                                    <  2:2526304/150‑1 (MQ=255)
   cGGCACAGGGTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCa                                                                                                                                                    <  2:325049/150‑1 (MQ=255)
        caGGGTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGa                                                                                                                                               <  2:1287315/150‑1 (MQ=255)
        caGGGTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGa                                                                                                                                               <  2:1289032/150‑1 (MQ=255)
         aGGGTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTCTTCCGGTATGACCGCCAGGGGAg                                                                                                                                              >  2:1334904/1‑150 (MQ=255)
          gggTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGc                                                                                                                                             <  2:2585888/150‑1 (MQ=255)
          gggTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGc                                                                                                                                             <  2:3870613/150‑1 (MQ=255)
            gTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGa                                                                                                                                           <  2:1185233/150‑1 (MQ=255)
             ttCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAg                                                                                                                                          <  1:15842/150‑1 (MQ=255)
                gAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGAt                                                                                                                                       <  2:3302484/150‑1 (MQ=255)
                gAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGAt                                                                                                                                       <  1:2593218/150‑1 (MQ=255)
                gAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGAt                                                                                                                                       >  2:4030451/1‑150 (MQ=255)
                   cccATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGAtgtg                                                                                                                                    <  1:3614376/150‑1 (MQ=255)
                     cATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCa                                                                                                                                  <  2:2570627/150‑1 (MQ=255)
                      aTGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCAt                                                                                                                                 <  1:2837393/150‑1 (MQ=255)
                      aTGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGAGGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCAt                                                                                                                                 <  2:3344474/150‑1 (MQ=255)
                       tGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATc                                                                                                                                >  1:4459330/1‑150 (MQ=255)
                        gTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCTCCAGGGGAGCGAGGATGTGCATCa                                                                                                                               <  1:3995824/150‑1 (MQ=255)
                        gTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCa                                                                                                                               >  1:23337/1‑150 (MQ=255)
                        gTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCa                                                                                                                               <  2:624108/150‑1 (MQ=255)
                            ggggTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACgg                                                                                                                           >  1:3248682/1‑150 (MQ=255)
                              ggTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTg                                                                                                                         <  1:770814/150‑1 (MQ=255)
                                tAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGaa                                                                                                                       >  1:3994331/1‑150 (MQ=255)
                                 aaTAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAAt                                                                                                                      >  1:2025598/1‑150 (MQ=255)
                                    cggTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCa                                                                                                                   <  2:1916503/149‑1 (MQ=255)
                                     ggTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCAt                                                                                                                  <  1:3746135/150‑1 (MQ=255)
                                      gtggcgtggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCAtt                                                                                                                 <  2:2926715/150‑1 (MQ=255)
                                       tggcgtggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCAttt                                                                                                                <  1:370783/150‑1 (MQ=255)
                                       tggcgtggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCAttt                                                                                                                <  1:2241658/150‑1 (MQ=255)
                                        ggcgtggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTtc                                                                                                               >  2:2283804/1‑150 (MQ=255)
                                          cgtggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTtctc                                                                                                             <  1:4240565/150‑1 (MQ=255)
                                          cgtggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCCCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTtctc                                                                                                             >  2:1762224/1‑150 (MQ=255)
                                            tggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCt                                                                                                           <  2:2141682/150‑1 (MQ=255)
                                             ggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCtt                                                                                                          <  2:4236475/150‑1 (MQ=255)
                                              gcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCttt                                                                                                         <  1:2042560/150‑1 (MQ=255)
                                               ccGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTg                                                                                                        <  2:2619204/150‑1 (MQ=255)
                                                cGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCATTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGa                                                                                                       <  2:4146178/150‑1 (MQ=255)
                                                  agaACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCTCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAg                                                                                                     <  1:1069652/150‑1 (MQ=255)
                                                    aaCAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGt                                                                                                   >  1:4284270/1‑150 (MQ=255)
                                                        aaaaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGcc                                                                                               <  2:480172/150‑1 (MQ=255)
                                                         aaaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGccc                                                                                              <  2:2381373/150‑1 (MQ=255)
                                                         aaaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGccc                                                                                              >  2:1599842/1‑150 (MQ=255)
                                                               tCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaa                                                                                        <  2:3764483/150‑1 (MQ=255)
                                                                cTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCTCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaa                                                                                       <  1:2637564/150‑1 (MQ=255)
                                                                 tGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaat                                                                                      <  1:639341/150‑1 (MQ=255)
                                                                  ggCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCTCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaata                                                                                     <  2:198866/150‑1 (MQ=255)
                                                                     aGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaa                                                                                  >  1:3323916/1‑150 (MQ=255)
                                                                      gCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAAGCGTCCTCACACGCACATCCCTTGGTATTCCGGTATGACCGCCAGGGCAGAGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaat                                                                                 >  1:2927077/1‑150 (MQ=255)
                                                                       cccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaata                                                                                <  1:3008387/150‑1 (MQ=255)
                                                                        ccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAatat                                                                               <  1:1445356/150‑1 (MQ=255)
                                                                        ccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAatat                                                                               <  1:4231761/150‑1 (MQ=255)
                                                                         cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATg                                                                              <  1:4482906/150‑1 (MQ=255)
                                                                         cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATg                                                                              >  2:3714115/1‑150 (MQ=255)
                                                                         cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATg                                                                              <  1:4416624/150‑1 (MQ=255)
                                                                            cAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGccc                                                                           <  1:953222/150‑1 (MQ=255)
                                                                                  cAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCTCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGc                                                                     >  2:3462663/1‑150 (MQ=255)
                                                                                     gTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCttt                                                                  <  1:1182053/150‑1 (MQ=255)
                                                                                      tGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTg                                                                 <  2:62490/150‑1 (MQ=255)
                                                                                      tGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCAGTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTg                                                                 <  2:489545/150‑1 (MQ=255)
                                                                                        cTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGcc                                                               <  1:1038960/150‑1 (MQ=255)
                                                                                         tCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCt                                                              <  1:2728932/150‑1 (MQ=255)
                                                                                          cGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCtt                                                             >  2:4369028/1‑150 (MQ=255)
                                                                                           gATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTc                                                            >  2:4209330/1‑150 (MQ=255)
                                                                                           gATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTc                                                            >  2:4209574/1‑150 (MQ=255)
                                                                                           gATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTc                                                            >  1:3499657/1‑150 (MQ=255)
                                                                                            aTGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCa                                                           >  1:2270646/1‑150 (MQ=255)
                                                                                            aTGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCa                                                           <  1:1687171/150‑1 (MQ=255)
                                                                                                 aGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCg                                                      >  1:3661882/1‑150 (MQ=255)
                                                                                                   tGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCggg                                                    >  2:2554378/1‑150 (MQ=255)
                                                                                                    gCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGc                                                   >  1:2245028/1‑150 (MQ=255)
                                                                                                      gTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCTCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTa                                                 <  1:3144463/150‑1 (MQ=255)
                                                                                                      gTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTa                                                 <  1:1587910/150‑1 (MQ=255)
                                                                                                      gTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTa                                                 <  2:1884758/150‑1 (MQ=255)
                                                                                                       tATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCtcg                                                >  2:1735029/1‑148 (MQ=255)
                                                                                                          cAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTagcag                                             >  2:141730/1‑150 (MQ=255)
                                                                                                           aaCCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCTCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTagcagc                                            <  1:1798671/150‑1 (MQ=255)
                                                                                                              cGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCTCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGAt                                         >  1:3608971/1‑150 (MQ=255)
                                                                                                                 cctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAAc                                      >  2:1506053/1‑150 (MQ=255)
                                                                                                                  ctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAAcg                                     <  2:100466/150‑1 (MQ=255)
                                                                                                                    cccACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAAcgcg                                   >  2:1867305/1‑150 (MQ=255)
                                                                                                                      cACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCggg                                 >  1:3574231/1‑150 (MQ=255)
                                                                                                                        cGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGgc                               <  2:682412/150‑1 (MQ=255)
                                                                                                                           ccATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGgcgcg                            >  1:1412638/1‑150 (MQ=255)
                                                                                                                           ccATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGgcgcg                            >  1:3733463/1‑150 (MQ=255)
                                                                                                                           ccATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGgcgcg                            >  1:3715976/1‑150 (MQ=255)
                                                                                                                               cccTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTCCATAAAATATGCCCGACGGCTTTGCCTGCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCt                        >  2:2487624/1‑150 (MQ=255)
                                                                                                                                   tGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCa                    >  2:2756512/1‑150 (MQ=255)
                                                                                                                                      tATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATcc                 <  2:2035102/150‑1 (MQ=255)
                                                                                                                                        ttCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTc               >  1:1393391/1‑150 (MQ=255)
                                                                                                                                           cGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAg            <  2:4457595/150‑1 (MQ=255)
                                                                                                                                             gTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTa          >  2:2448537/1‑150 (MQ=255)
                                                                                                                                                 gACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAaga      >  1:4170574/1‑150 (MQ=255)
                                                                                                                                                  aCCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAagag     <  2:3573503/150‑1 (MQ=255)
                                                                                                                                                    cGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAAGAGGc   <  1:1879/150‑1 (MQ=255)
                                                                                                                                                    cGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAAGAGGc   <  2:2244609/150‑1 (MQ=255)
                                                                                                                                                     gCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAAGAGGcc  >  2:3912234/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCCGGCACAGGGTTCGAACCCATGTCGGGGGTAATAGGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCTAGCAGCGATAACGCGGGGCGCGCGCTTTCATCCTCAAGTAAAGAGGCC  >  CP009273/4523196‑4523494

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: