Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,542,876 T→G 27.5% N185T (AAC→ACC)  malQ ← 4‑alpha‑glucanotransferase (amylomaltase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,542,8760TG27.5% 395.1 / 290.6 305N185T (AAC→ACC) malQ4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base T (104/117);  new base G (38/46);  total (142/163)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 309 total aligned reads displayed.

GCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCC  >  CP009273/3542727‑3543025
                                                                                                                                                     |                                                                                                                                                      
   aCTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCa                                                                                                                                                     <  1:2267740/150‑1 (MQ=255)
        caccaGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCa                                                                                                                                                >  1:2218481/1‑150 (MQ=255)
        caccaGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCa                                                                                                                                                >  1:629832/1‑150 (MQ=255)
         accaGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCaa                                                                                                                                               >  1:880998/1‑150 (MQ=255)
          ccaGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAAt                                                                                                                                              <  2:1353028/150‑1 (MQ=255)
            aGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGa                                                                                                                                            >  2:116000/1‑150 (MQ=255)
            aGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGa                                                                                                                                            <  2:3139708/150‑1 (MQ=255)
             ggCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGaa                                                                                                                                           >  2:3809476/1‑150 (MQ=255)
              gCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAAc                                                                                                                                          <  2:2534656/150‑1 (MQ=255)
               ccTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAACTCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACg                                                                                                                                         <  2:2585600/150‑1 (MQ=255)
               ccTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACg                                                                                                                                         <  1:1066476/150‑1 (MQ=255)
               ccTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACg                                                                                                                                         <  2:1915585/150‑1 (MQ=255)
                  gAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGAcc                                                                                                                                      >  1:200096/1‑150 (MQ=255)
                    gCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCg                                                                                                                                    >  1:287736/1‑150 (MQ=255)
                         ttCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACg                                                                                                                               <  2:886675/150‑1 (MQ=255)
                           cGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGtt                                                                                                                             >  2:2825935/1‑150 (MQ=255)
                             cTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGtttt                                                                                                                           <  2:3198635/150‑1 (MQ=255)
                             cTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGtttt                                                                                                                           <  1:3673409/150‑1 (MQ=255)
                             cTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGtttt                                                                                                                           <  2:3932063/150‑1 (MQ=255)
                               aaGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGc                                                                                                                         >  1:2245139/1‑150 (MQ=255)
                                aGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGcc                                                                                                                        <  2:3678816/150‑1 (MQ=255)
                                 gATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCca                                                                                                                       <  1:4332448/150‑1 (MQ=255)
                                  aTGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCcac                                                                                                                      >  1:3033060/1‑150 (MQ=255)
                                   tGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCcaca                                                                                                                     >  1:1161975/1‑150 (MQ=255)
                                     gAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATc                                                                                                                   >  1:1140192/1‑150 (MQ=255)
                                        aTCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATccac                                                                                                                <  2:4107616/150‑1 (MQ=255)
                                        aTCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATccac                                                                                                                >  1:2806462/1‑150 (MQ=255)
                                        aTCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATccac                                                                                                                >  1:1853125/1‑150 (MQ=255)
                                        aTCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATccac                                                                                                                >  1:1798782/1‑150 (MQ=255)
                                          cTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATccacca                                                                                                              <  1:1586752/150‑1 (MQ=255)
                                              aaCGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGAAAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCAt                                                                                                          >  1:3246539/1‑150 (MQ=255)
                                                cGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCg                                                                                                        <  1:3902166/150‑1 (MQ=255)
                                                  gCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCt                                                                                                      <  1:4255212/150‑1 (MQ=255)
                                                   cGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCtt                                                                                                     >  2:811212/1‑150 (MQ=255)
                                                   cGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCtt                                                                                                     >  2:1256524/1‑150 (MQ=255)
                                                    gTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGTGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCttt                                                                                                    <  2:4308370/150‑1 (MQ=255)
                                                    gTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCttt                                                                                                    >  1:2869573/1‑150 (MQ=255)
                                                    gTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCttt                                                                                                    <  2:583811/150‑1 (MQ=255)
                                                      tAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTga                                                                                                  >  1:2989461/1‑150 (MQ=255)
                                                       aaCGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTgag                                                                                                 <  1:2950742/150‑1 (MQ=255)
                                                         cGTCGATATAAATCACATTCAGCCCACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTAGCCACATACACAAGCATCGCTTAGAGAt                                                                                               >  2:2331036/1‑150 (MQ=255)
                                                          gTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATc                                                                                              <  2:158030/150‑1 (MQ=255)
                                                            cGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGc                                                                                            >  1:2635410/1‑150 (MQ=255)
                                                               tataAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCaa                                                                                         >  1:784505/1‑150 (MQ=255)
                                                                ataAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCaaa                                                                                        >  2:2363564/1‑150 (MQ=255)
                                                                 taAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCaaaa                                                                                       >  2:3513963/1‑150 (MQ=255)
                                                                  aaaTCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAt                                                                                      >  1:793500/1‑150 (MQ=255)
                                                                    aTCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATcc                                                                                    <  1:2015245/150‑1 (MQ=255)
                                                                     tCACATTCAGCCAACGGCGAGAAGACGGGCTGTAGGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATccc                                                                                   <  2:1199466/150‑1 (MQ=255)
                                                                        caTTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCaa                                                                                >  1:2502621/1‑150 (MQ=255)
                                                                            cAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATAcc                                                                            >  2:145257/1‑150 (MQ=255)
                                                                            cAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATAcc                                                                            >  2:412214/1‑150 (MQ=255)
                                                                              gCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATAcccc                                                                          >  1:3375334/1‑150 (MQ=255)
                                                                                cAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGACAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAg                                                                        <  2:120493/150‑1 (MQ=255)
                                                                                   cGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGttt                                                                    <  1:3272177/150‑1 (MQ=255)
                                                                                   cGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGttt                                                                    <  2:2726030/150‑1 (MQ=255)
                                                                                    ggCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGtttt                                                                   <  1:3861563/150‑1 (MQ=255)
                                                                                    ggCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGtttt                                                                   >  2:2152275/1‑150 (MQ=255)
                                                                                       gagaAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGttttttt                                                                <  1:1193846/150‑1 (MQ=255)
                                                                                       gagaAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGttttttt                                                                >  1:1004883/1‑150 (MQ=255)
                                                                                        agaagaCGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTc                                                               <  1:4385436/150‑1 (MQ=255)
                                                                                         gaagaCGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTcc                                                              <  1:1602806/150‑1 (MQ=255)
                                                                                         gaagaCGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTcc                                                              >  2:2044418/1‑150 (MQ=255)
                                                                                             aCGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATc                                                          <  2:1545332/150‑1 (MQ=255)
                                                                                              cGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCg                                                         >  2:766467/1‑150 (MQ=255)
                                                                                               gggCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGc                                                        <  1:829944/150‑1 (MQ=255)
                                                                                                     tATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGATTCGGGTTCATGACAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCACAATACCCCAGTTTTTTTCCGATCGCAGCGta                                                  <  2:3252821/150‑1 (MQ=255)
                                                                                                       tGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGTACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATACCCAATACCCCAGTTTTTTTCCGATCGCAGCGtata                                                >  2:3034303/1‑150 (MQ=255)
                                                                                                       tGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGtata                                                <  2:3590255/150‑1 (MQ=255)
                                                                                                        gggCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATaa                                               <  2:2213690/150‑1 (MQ=255)
                                                                                                          gCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAg                                             <  1:760271/150‑1 (MQ=255)
                                                                                                          gCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAg                                             <  1:2568361/150‑1 (MQ=255)
                                                                                                          gCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAg                                             <  1:3765639/150‑1 (MQ=255)
                                                                                                             ggCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTg                                          >  1:1495512/1‑150 (MQ=255)
                                                                                                              gcgcTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATCCCCCAGTTTTTTTTCCGATCGCAGCGTATAAAGCTg                                          <  2:344044/150‑1 (MQ=255)
                                                                                                              gcgcTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTTCCGATCGCAGCGTATAAAGCTg                                          <  2:3469185/150‑1 (MQ=255)
                                                                                                               cgcTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGaa                                        >  1:4377934/1‑150 (MQ=255)
                                                                                                                gcTCTCCGGATGTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCTCGTTTTGACACATCCAACAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCTGCGTATAAAGCTGAAc                                       <  1:330320/150‑1 (MQ=255)
                                                                                                                 ctctCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACg                                      >  1:3577015/1‑150 (MQ=255)
                                                                                                                      cGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGc                                 >  2:4442494/1‑150 (MQ=255)
                                                                                                                         aTTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCAcc                              <  2:1057666/150‑1 (MQ=255)
                                                                                                                         aTTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCAcc                              <  1:4101359/150‑1 (MQ=255)
                                                                                                                              ccGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCcaca                         <  2:681914/150‑1 (MQ=255)
                                                                                                                              ccGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCcaca                         <  1:4523703/150‑1 (MQ=255)
                                                                                                                                ggATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGc                       >  1:956109/1‑150 (MQ=255)
                                                                                                                                 gATAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCt                      >  2:2908823/1‑150 (MQ=255)
                                                                                                                                  aTAGAGCGCATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCtt                     <  1:2420944/150‑1 (MQ=255)
                                                                                                                                   tAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCttt                    >  1:3107733/1‑150 (MQ=255)
                                                                                                                                   tAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCttt                    <  1:1654630/150‑1 (MQ=255)
                                                                                                                                    agagCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCtttt                   >  1:600221/1‑150 (MQ=255)
                                                                                                                                        cgcATGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGttt               <  2:2235967/150‑1 (MQ=255)
                                                                                                                                          cATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTAt             >  2:1402394/1‑150 (MQ=255)
                                                                                                                                           aTGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTc           <  2:1720766/150‑1 (MQ=255)
                                                                                                                                            tGAATCGGGGTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTc           <  1:1295563/150‑1 (MQ=255)
                                                                                                                                             gAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCa          <  2:1809478/150‑1 (MQ=255)
                                                                                                                                               aTCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGc        >  2:4002788/1‑150 (MQ=255)
                                                                                                                                               aTCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGc        <  2:83945/150‑1 (MQ=255)
                                                                                                                                                  ggggTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAAc     <  1:730904/150‑1 (MQ=255)
                                                                                                                                                  gggTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAAc     >  1:2132376/1‑150 (MQ=255)
                                                                                                                                                   ggTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACg    <  1:233516/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCC  >  CP009273/3542727‑3543025

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: