Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,542,744 T→A 5.4% Q229L (CAG→CTG)  malQ ← 4‑alpha‑glucanotransferase (amylomaltase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,542,7440TA 5.4% 1152.6 / 41.0 334Q229L (CAG→CTG) malQ4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base T (145/171);  new base A (11/7);  total (156/178)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 336 total aligned reads displayed.

CCATCTGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGAC  >  CP009273/3542595‑3542893
                                                                                                                                                     |                                                                                                                                                     
 cATCTGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTg                                                                                                                                                      >  2:2661560/1‑150 (MQ=255)
 cATCTGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTg                                                                                                                                                      <  2:1687824/150‑1 (MQ=255)
 cATCTGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTg                                                                                                                                                      <  2:3933314/150‑1 (MQ=255)
 cATCTGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTg                                                                                                                                                      <  2:3593029/150‑1 (MQ=255)
    cTGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGc                                                                                                                                                   >  1:1042952/1‑150 (MQ=255)
     tGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGcc                                                                                                                                                  <  1:4394539/150‑1 (MQ=255)
       cTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCctc                                                                                                                                                <  1:966929/150‑1 (MQ=255)
         catcatcaCGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCtt                                                                                                                                              >  1:2293411/1‑150 (MQ=255)
          atcatcaCGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTc                                                                                                                                             <  1:2444213/150‑1 (MQ=255)
           tcatcaCGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCg                                                                                                                                            <  1:1582291/150‑1 (MQ=255)
                aCGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAg                                                                                                                                       >  2:3667325/1‑150 (MQ=255)
                 cGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGa                                                                                                                                      >  1:3946313/1‑150 (MQ=255)
                     gcgcGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGa                                                                                                                                  <  1:2788609/150‑1 (MQ=255)
                       gcgAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGaaa                                                                                                                                >  2:13325/1‑150 (MQ=255)
                       gcgAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGaaa                                                                                                                                <  2:4415489/150‑1 (MQ=255)
                        cgAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAAt                                                                                                                               <  1:222569/150‑1 (MQ=255)
                         gAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATc                                                                                                                              <  2:3649399/150‑1 (MQ=255)
                         gAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATc                                                                                                                              >  2:1256520/1‑150 (MQ=255)
                          aaaCCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCt                                                                                                                             >  1:116881/1‑150 (MQ=255)
                          aaaCCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCt                                                                                                                             >  2:1394597/1‑150 (MQ=255)
                           aaCCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCtt                                                                                                                            >  1:788841/1‑150 (MQ=255)
                            aCCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTc                                                                                                                           <  1:2916244/150‑1 (MQ=255)
                            aCCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTc                                                                                                                           >  2:2229624/1‑150 (MQ=255)
                            aCCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTc                                                                                                                           >  1:2963646/1‑150 (MQ=255)
                              cTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCaa                                                                                                                         >  1:2005015/1‑150 (MQ=255)
                              cTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCaa                                                                                                                         >  2:4778342/1‑150 (MQ=255)
                                 tCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACgg                                                                                                                      >  2:3777612/1‑150 (MQ=255)
                                 tCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACgg                                                                                                                      <  1:505074/150‑1 (MQ=255)
                                      gCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTa                                                                                                                 >  2:1627604/1‑150 (MQ=255)
                                       ccATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGTAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTaa                                                                                                                >  1:2528171/1‑150 (MQ=255)
                                       ccATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGTAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTaa                                                                                                                >  1:3451526/1‑150 (MQ=255)
                                       ccATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTaa                                                                                                                <  2:723190/150‑1 (MQ=255)
                                           tCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTc                                                                                                            >  1:1330039/1‑150 (MQ=255)
                                           tCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTc                                                                                                            >  1:1485117/1‑150 (MQ=255)
                                           tCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTc                                                                                                            <  1:680324/150‑1 (MQ=255)
                                             gTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGa                                                                                                          <  1:4046269/150‑1 (MQ=255)
                                              tAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGat                                                                                                         >  1:896472/1‑150 (MQ=255)
                                               aaTGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGata                                                                                                        >  1:1384800/1‑150 (MQ=255)
                                                    tGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAAt                                                                                                   <  2:1752473/150‑1 (MQ=255)
                                                    tGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAAt                                                                                                   >  1:4144653/1‑150 (MQ=255)
                                                     gTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATc                                                                                                  <  2:2978922/150‑1 (MQ=255)
                                                      tcttttttAGGGCGGTATCCGTGGAGTTATCGACCGAGTCGGCTTCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCTACTGCCTCCAGGCCTGAGCCTCTTCGCTAAGTTGTTTTTATTCAAGGGCGTTATCGTCGATATATATca                                                                                                 <  2:151949/147‑1 (MQ=255)
                                                       cATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATcac                                                                                                <  1:3695371/150‑1 (MQ=255)
                                                       cATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATcac                                                                                                >  2:3255061/1‑150 (MQ=255)
                                                        aTTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATcaca                                                                                               >  1:4302358/1‑150 (MQ=255)
                                                         tttttAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACAt                                                                                              <  1:1368867/150‑1 (MQ=255)
                                                           tttAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTc                                                                                            >  1:2928836/1‑150 (MQ=255)
                                                            ttAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCa                                                                                           <  2:445200/150‑1 (MQ=255)
                                                               gggCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGcc                                                                                        >  1:184975/1‑150 (MQ=255)
                                                                ggCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCa                                                                                       <  2:3061137/150‑1 (MQ=255)
                                                                ggCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCa                                                                                       <  1:268841/150‑1 (MQ=255)
                                                                  cGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAAc                                                                                     >  2:2089444/1‑150 (MQ=255)
                                                                  cGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAAc                                                                                     <  1:4092741/150‑1 (MQ=255)
                                                                  cGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAAc                                                                                     >  1:3662571/1‑150 (MQ=255)
                                                                    gTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACgg                                                                                   <  1:1918823/150‑1 (MQ=255)
                                                                     tAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGc                                                                                  <  2:4584847/150‑1 (MQ=255)
                                                                         cGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCgaga                                                                              >  2:1873820/1‑150 (MQ=255)
                                                                          gTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGagaa                                                                             >  2:2588737/1‑150 (MQ=255)
                                                                            ggAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGagaaga                                                                           <  2:1238732/150‑1 (MQ=255)
                                                                            ggAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGagaaga                                                                           >  1:4597766/1‑150 (MQ=255)
                                                                            ggAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGagaaga                                                                           <  2:31889/150‑1 (MQ=255)
                                                                             gAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGAc                                                                          >  2:2840028/1‑150 (MQ=255)
                                                                               gTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACgg                                                                        >  2:2888048/1‑150 (MQ=255)
                                                                                tAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACggg                                                                       >  1:1257362/1‑150 (MQ=255)
                                                                                  aTCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCt                                                                     >  2:2331719/1‑150 (MQ=255)
                                                                                    cGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGt                                                                   <  2:4761081/150‑1 (MQ=255)
                                                                                    cGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGt                                                                   >  2:3514052/1‑150 (MQ=255)
                                                                                    cGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGt                                                                   >  2:1108105/1‑150 (MQ=255)
                                                                                     gACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTa                                                                  >  2:414811/1‑150 (MQ=255)
                                                                                     gACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTa                                                                  >  1:4432697/1‑150 (MQ=255)
                                                                                      aCCCTGTCGGCAGCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTAt                                                                 <  1:433068/150‑1 (MQ=255)
                                                                                      aCCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTAt                                                                 >  2:1839501/1‑150 (MQ=255)
                                                                                          aGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGc                                                             >  2:2914011/1‑150 (MQ=255)
                                                                                            tCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTg                                                           >  1:3265343/1‑150 (MQ=255)
                                                                                                   cgcgcgCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGctct                                                    <  2:3674087/150‑1 (MQ=255)
                                                                                                      gcgcCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCg                                                 >  1:3085297/1‑150 (MQ=255)
                                                                                                      gcgcCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGAGGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCg                                                 <  2:2721502/150‑1 (MQ=255)
                                                                                                         ccTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGAt                                              >  2:1955894/1‑150 (MQ=255)
                                                                                                          cTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGAtt                                             >  1:1432397/1‑150 (MQ=255)
                                                                                                           tGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGAttt                                            >  2:471129/1‑150 (MQ=255)
                                                                                                            gTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTg                                           <  2:3807074/150‑1 (MQ=255)
                                                                                                              tGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGcc                                         >  2:4184358/1‑150 (MQ=255)
                                                                                                                cAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCgg                                       >  2:2501672/1‑150 (MQ=255)
                                                                                                                cAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCgg                                       >  2:4406333/1‑150 (MQ=255)
                                                                                                                 aGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGa                                      >  1:3045377/1‑150 (MQ=255)
                                                                                                                   cGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGTGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATa                                    <  1:1599818/150‑1 (MQ=255)
                                                                                                                      cTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATagag                                 >  1:3300747/1‑150 (MQ=255)
                                                                                                                       tGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAgc                                >  2:1327840/1‑150 (MQ=255)
                                                                                                                          tGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCa                             <  1:1747064/150‑1 (MQ=255)
                                                                                                                               ggTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAAt                        <  1:714599/150‑1 (MQ=255)
                                                                                                                               ggTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAAt                        <  2:432500/150‑1 (MQ=255)
                                                                                                                               ggTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAAt                        <  2:3032704/150‑1 (MQ=255)
                                                                                                                                     cAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGtt                  <  1:3558129/150‑1 (MQ=255)
                                                                                                                                        cTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAg               >  1:3706921/1‑150 (MQ=255)
                                                                                                                                            caccaGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCa           <  2:659424/150‑1 (MQ=255)
                                                                                                                                            caccaGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCa           <  1:4772626/150‑1 (MQ=255)
                                                                                                                                            caccaGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGAGTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCa           >  2:800849/1‑150 (MQ=255)
                                                                                                                                              ccaGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAAt         <  1:3546233/150‑1 (MQ=255)
                                                                                                                                              ccaGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAAt         >  2:3501219/1‑150 (MQ=255)
                                                                                                                                                     tGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGAc  <  1:341585/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCATCTGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGAC  >  CP009273/3542595‑3542893

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: ‑