Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,542,744 T→A 11.3% 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,7440TA11.3% 1039.6 / 104.0 327Q229L (CAG→CTG) malQ4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base T (136/154);  new base A (7/30);  total (143/184)
Fisher's exact test for biased strand distribution p-value = 1.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 329 total aligned reads displayed.

CATCTGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGAC  >  CP009273/3542596‑3542893
                                                                                                                                                    |                                                                                                                                                     
cATCTGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTg                                                                                                                                                      <  1:2967783/150‑1 (MQ=255)
cATCTGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTg                                                                                                                                                      <  1:128464/150‑1 (MQ=255)
  tCTGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAg                                                                                                                                                    <  1:2073854/150‑1 (MQ=255)
      cTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCctc                                                                                                                                                <  2:3206392/150‑1 (MQ=255)
          tcatcaCGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCg                                                                                                                                            <  2:1540014/150‑1 (MQ=255)
                 gTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGAt                                                                                                                                     <  1:1313572/150‑1 (MQ=255)
                   tGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATgg                                                                                                                                   <  1:185752/150‑1 (MQ=255)
                      gcgAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGaaa                                                                                                                                <  1:237379/150‑1 (MQ=255)
                       cgAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAAt                                                                                                                               <  2:4307368/150‑1 (MQ=255)
                        gAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATc                                                                                                                              >  1:3614459/1‑150 (MQ=255)
                         aaaCCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCt                                                                                                                             <  2:2802867/150‑1 (MQ=255)
                           aCCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTc                                                                                                                           <  1:815070/150‑1 (MQ=255)
                            ccTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCa                                                                                                                          >  2:3959539/1‑150 (MQ=255)
                              tttCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAAc                                                                                                                        <  1:1009715/150‑1 (MQ=255)
                              tttCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAAc                                                                                                                        <  1:1111774/150‑1 (MQ=255)
                                tCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACgg                                                                                                                      >  2:3064465/1‑150 (MQ=255)
                                 ccACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGc                                                                                                                     <  2:464125/150‑1 (MQ=255)
                                   aCGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGt                                                                                                                   <  1:765544/150‑1 (MQ=255)
                                    cGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGtt                                                                                                                  <  2:2296170/150‑1 (MQ=255)
                                    cGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGtt                                                                                                                  >  1:4189172/1‑150 (MQ=255)
                                     gCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTa                                                                                                                 >  2:1886776/1‑150 (MQ=255)
                                     gCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTa                                                                                                                 >  1:2649960/1‑150 (MQ=255)
                                     gCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTa                                                                                                                 <  2:4110908/150‑1 (MQ=255)
                                      ccATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTaa                                                                                                                >  2:540795/1‑150 (MQ=255)
                                      ccATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTaa                                                                                                                <  1:3196098/150‑1 (MQ=255)
                                        tttCGTAATGCTGTCATTTTTCGGGGGGTAACCGTGGATTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACg                                                                                                              <  1:734099/149‑1 (MQ=255)
                                          tCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTc                                                                                                            <  2:2981247/150‑1 (MQ=255)
                                           cGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCg                                                                                                           >  2:69088/1‑150 (MQ=255)
                                           cGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCg                                                                                                           >  2:810589/1‑150 (MQ=255)
                                             tAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGat                                                                                                         <  1:4107379/150‑1 (MQ=255)
                                             tAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGat                                                                                                         >  1:83268/1‑150 (MQ=255)
                                             tAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGat                                                                                                         <  2:4274749/150‑1 (MQ=255)
                                              aaTGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGata                                                                                                        <  1:3897717/150‑1 (MQ=255)
                                                 gCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATaa                                                                                                     >  1:874540/1‑150 (MQ=255)
                                                  cTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATaaa                                                                                                    >  1:2719294/1‑150 (MQ=255)
                                                   tGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAAt                                                                                                   >  2:1183613/1‑150 (MQ=255)
                                                   tGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAAt                                                                                                   <  1:117427/150‑1 (MQ=255)
                                                    gTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATc                                                                                                  <  1:2200209/150‑1 (MQ=255)
                                                     tCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATca                                                                                                 <  2:1341902/150‑1 (MQ=255)
                                                     tCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATca                                                                                                 >  1:1090897/1‑150 (MQ=255)
                                                     tCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCAACCAGGCCTGAGCCTCTTCGCTAAGATGTAAATCTTCAACGGCGTTAACGTCGATATAAATca                                                                                                 <  1:90955/150‑1 (MQ=255)
                                                      cATTTAAAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGAGTGGTCGGCAACTGCCACCAGTCCTGAGCCTCTACGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATcac                                                                                                >  1:2267853/1‑150 (MQ=255)
                                                       aTTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATcaca                                                                                               >  2:869976/1‑150 (MQ=255)
                                                         ttttAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACAtt                                                                                             <  2:3866042/150‑1 (MQ=255)
                                                         ttttAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACAtt                                                                                             <  2:1678020/150‑1 (MQ=255)
                                                           ttAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCa                                                                                           >  1:1762906/1‑150 (MQ=255)
                                                            tAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAg                                                                                          <  1:491214/150‑1 (MQ=255)
                                                            tAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAg                                                                                          <  1:3054342/150‑1 (MQ=255)
                                                             aGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGc                                                                                         >  2:4181949/1‑150 (MQ=255)
                                                              gggCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGcc                                                                                        <  2:1281436/150‑1 (MQ=255)
                                                               ggCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCa                                                                                       <  1:1247515/150‑1 (MQ=255)
                                                                 cGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAAc                                                                                     <  1:4341153/150‑1 (MQ=255)
                                                                 cGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAAc                                                                                     >  1:2810496/1‑150 (MQ=255)
                                                                  ggTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACg                                                                                    <  2:2488756/150‑1 (MQ=255)
                                                                  ggTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACg                                                                                    >  2:1439388/1‑150 (MQ=255)
                                                                     aaCCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCg                                                                                 >  2:1346248/1‑150 (MQ=255)
                                                                     aaCCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCg                                                                                 <  2:2563883/150‑1 (MQ=255)
                                                                        cGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCgaga                                                                              >  2:3887925/1‑150 (MQ=255)
                                                                         gTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGagaa                                                                             >  1:3709469/1‑150 (MQ=255)
                                                                          tGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAGATCACATTCAGCCAACGGCGagaag                                                                            <  2:827892/150‑1 (MQ=255)
                                                                          tGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGagaag                                                                            <  2:3289385/150‑1 (MQ=255)
                                                                            gtgTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCGGTTGCGTGGTCGGCAACTGCCACCTGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGAc                                                                          <  1:2957603/148‑1 (MQ=255)
                                                                            gAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGAc                                                                          <  2:833057/150‑1 (MQ=255)
                                                                                 aTCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCt                                                                     >  2:1073209/1‑150 (MQ=255)
                                                                                   cGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGt                                                                   <  1:137291/150‑1 (MQ=255)
                                                                                    gACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTa                                                                  >  2:1086150/1‑150 (MQ=255)
                                                                                       ccAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTGAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATgg                                                               >  2:1608418/1‑150 (MQ=255)
                                                                                        cAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATggg                                                              <  2:3527201/150‑1 (MQ=255)
                                                                                         aGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGc                                                             >  2:4236061/1‑150 (MQ=255)
                                                                                           tCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTg                                                           <  1:2380445/150‑1 (MQ=255)
                                                                                             ggCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGgc                                                         >  2:1646016/1‑150 (MQ=255)
                                                                                             ggCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGgc                                                         >  1:145888/1‑150 (MQ=255)
                                                                                                aTCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGct                                                      <  2:590185/150‑1 (MQ=255)
                                                                                                   gcgcgcCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCtctc                                                   <  1:628936/150‑1 (MQ=255)
                                                                                                      cgcCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCgg                                                >  2:3329637/1‑150 (MQ=255)
                                                                                                      cgcCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCgg                                                >  1:176494/1‑150 (MQ=255)
                                                                                                         cTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGAtt                                             >  1:2026707/1‑150 (MQ=255)
                                                                                                          tGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGAttt                                            <  1:3440691/150‑1 (MQ=255)
                                                                                                          tGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGAttt                                            <  2:1975642/150‑1 (MQ=255)
                                                                                                           gTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTg                                           >  1:2859091/1‑150 (MQ=255)
                                                                                                           gTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTg                                           <  2:2058544/150‑1 (MQ=255)
                                                                                                              gCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGGGCCTCTTCGCTAAGATGGAAATCTTCAACGGAGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCg                                        >  2:337941/1‑150 (MQ=255)
                                                                                                                aGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGa                                      >  2:4348029/1‑150 (MQ=255)
                                                                                                                          gCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCAt                            <  1:1099409/150‑1 (MQ=255)
                                                                                                                           cGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATg                           <  2:2193891/150‑1 (MQ=255)
                                                                                                                              ggTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAAt                        <  2:287840/150‑1 (MQ=255)
                                                                                                                                    cAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGtt                  <  1:4294593/150‑1 (MQ=255)
                                                                                                                                     aaCTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTc                 <  1:3235725/150‑1 (MQ=255)
                                                                                                                                      aCTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCTGGATAGAGCGCATGAATCGGGTTCa                <  1:1279398/150‑1 (MQ=255)
                                                                                                                                      aCTGCCACCAGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCa                <  2:2580347/150‑1 (MQ=255)
                                                                                                                                       cTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAg               >  2:1860932/1‑150 (MQ=255)
                                                                                                                                       cTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAg               <  1:1641414/150‑1 (MQ=255)
                                                                                                                                        tGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAgg              <  1:3398356/150‑1 (MQ=255)
                                                                                                                                        tGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAgg              <  1:1205356/150‑1 (MQ=255)
                                                                                                                                          ccaccaGGCCAGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGcc            >  1:4013224/1‑150 (MQ=255)
                                                                                                                                               aGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGa       >  1:251549/1‑150 (MQ=255)
                                                                                                                                               aGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGa       >  1:3549764/1‑150 (MQ=255)
                                                                                                                                                cgCCTGAGCCTCTGCGCTAAGATGGAAATCATCTACGGCGTTAACGTCGATATAAATCACATTCAGCCAACTGCGAGAAGACGGGCTGTATTGGCAGGCGCTCTCCGGATTTTCCGGATAGACCGCATGACTCGGGTTGAGCCCAAAGac      >  2:2177950/2‑149 (MQ=255)
                                                                                                                                                  ccTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACg    >  2:3806637/1‑150 (MQ=255)
                                                                                                                                                    tGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGAc  >  1:1538718/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATCTGCTCATCATCACGTTGCGCGAAACCTTTCCACGCCATTCGTAATGCTGTCATTTTTAGGGCGGTAACCGTGGAGTAATCGACCCAGTCGGCATCGCGCGCCTGTTGCAGCGTCTGTTGCGTGGTCGGCAACTGCCACCAGGCCTGAGCCTCTTCGCTAAGATGGAAATCTTCAACGGCGTTAACGTCGATATAAATCACATTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGAC  >  CP009273/3542596‑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: ‑