Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,062,424 G→T 90.9% I71I (ATC→ATA)  yihQ ← alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,062,4240GT90.9% 890.8 / 75.6 298I71I (ATC→ATA) yihQalpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base G (13/14);  new base T (91/180);  total (104/194)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 298 total aligned reads displayed.

CGCAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGAAAATG  >  CP009273/4062278‑4062570
                                                                                                                                                  |                                                                                                                                                  
 gctgCCATATACGGTTGTGGTTATGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATgg                                                                                                                                                <  1:2674437/147‑1 (MQ=255)
 gcAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATgg                                                                                                                                                <  2:3304534/150‑1 (MQ=255)
 gcAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATgg                                                                                                                                                <  1:654870/150‑1 (MQ=255)
  cAGCCAGATAGGGTTGTGTTTATGGTTGTCTTTTTGTAGTTCCAGCAATACACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGc                                                                                                                                               <  1:2917264/150‑1 (MQ=255)
   aGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCg                                                                                                                                              >  2:926720/1‑150 (MQ=255)
     ccAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTc                                                                                                                                            <  2:2478391/150‑1 (MQ=255)
     ccAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTc                                                                                                                                            <  2:4168650/150‑1 (MQ=255)
       aGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCgg                                                                                                                                          >  1:1015190/1‑150 (MQ=255)
       aGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCgg                                                                                                                                          <  1:1771175/150‑1 (MQ=255)
        gATACGGTTGTGTTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGt                                                                                                                                         <  2:3613536/150‑1 (MQ=255)
         aTACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTa                                                                                                                                        >  2:4331515/1‑150 (MQ=255)
         aTACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGCCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTa                                                                                                                                        >  2:1930335/1‑150 (MQ=255)
          tACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTaa                                                                                                                                       >  2:3014748/1‑150 (MQ=255)
              gTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAgcgc                                                                                                                                   >  2:3143417/1‑150 (MQ=255)
                 gtgGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAAt                                                                                                                                <  2:706392/150‑1 (MQ=255)
                  tgGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCCGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAAtt                                                                                                                               <  1:535237/150‑1 (MQ=255)
                  tgGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCCGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAAtt                                                                                                                               <  1:536275/150‑1 (MQ=255)
                   ggTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAAttt                                                                                                                              >  1:25377/1‑150 (MQ=255)
                      tAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTtct                                                                                                                           >  1:1468376/1‑150 (MQ=255)
                         ggTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCt                                                                                                                        >  1:2899643/1‑150 (MQ=255)
                              tCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGt                                                                                                                   <  1:3659517/150‑1 (MQ=255)
                                gTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGttt                                                                                                                 <  2:693408/150‑1 (MQ=255)
                                  tttGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTAt                                                                                                               >  1:3605941/1‑150 (MQ=255)
                                     gTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCtt                                                                                                            >  2:4217225/1‑150 (MQ=255)
                                     gTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCtt                                                                                                            <  1:860031/150‑1 (MQ=255)
                                     gTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCtt                                                                                                            >  1:1516914/1‑150 (MQ=255)
                                       aGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTa                                                                                                          >  2:955519/1‑150 (MQ=255)
                                       aGTACCAGCAATAAACGCCCCTGATCGTCGGCATAGATATTCAGCGTGGCGCTAATGACAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGTCGTCTGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTa                                                                                                          >  2:1060723/1‑150 (MQ=255)
                                         ttCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAAt                                                                                                        <  2:4108205/150‑1 (MQ=255)
                                            cAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCt                                                                                                     >  2:277293/1‑150 (MQ=255)
                                             aGCAATACACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTg                                                                                                    <  2:2953971/150‑1 (MQ=255)
                                             aGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTg                                                                                                    >  1:827936/1‑150 (MQ=255)
                                                aaTAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGaaa                                                                                                 >  1:1396275/1‑150 (MQ=255)
                                                  tAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAAtt                                                                                               <  2:597804/150‑1 (MQ=255)
                                                  tAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAAtt                                                                                               >  1:2011281/1‑150 (MQ=255)
                                                      cGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATTCTTTAATTAAcg                                                                                           >  2:3642342/1‑150 (MQ=255)
                                                      cGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcg                                                                                           <  2:1165583/150‑1 (MQ=255)
                                                      cGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcg                                                                                           <  1:4237935/150‑1 (MQ=255)
                                                         cccTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCgg                                                                                        >  2:903276/1‑150 (MQ=255)
                                                         cccTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCgg                                                                                        >  2:2600870/1‑150 (MQ=255)
                                                          ccTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGa                                                                                       <  2:1253403/150‑1 (MQ=255)
                                                          ccTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGa                                                                                       <  1:611232/150‑1 (MQ=255)
                                                           cTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGaa                                                                                      >  1:4203425/1‑150 (MQ=255)
                                                            tGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTACGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAAc                                                                                     >  1:1187880/1‑150 (MQ=255)
                                                            tGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAAc                                                                                     >  1:3392458/1‑150 (MQ=255)
                                                               tcgtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACata                                                                                  <  2:1976175/150‑1 (MQ=255)
                                                               tcgtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACata                                                                                  <  2:1881222/150‑1 (MQ=255)
                                                                cgtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACatat                                                                                 <  1:1815766/150‑1 (MQ=255)
                                                                 gtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATc                                                                                <  1:1614810/150‑1 (MQ=255)
                                                                 gtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATc                                                                                <  2:1480080/150‑1 (MQ=255)
                                                                  tcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCg                                                                               >  1:1209423/1‑150 (MQ=255)
                                                                   cgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGa                                                                              <  1:2807851/150‑1 (MQ=255)
                                                                        gagaTATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATcc                                                                         >  2:88950/1‑150 (MQ=255)
                                                                         agaTATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCg                                                                        >  1:1547732/1‑150 (MQ=255)
                                                                           atatTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCt                                                                      <  2:4034136/150‑1 (MQ=255)
                                                                           atatTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCt                                                                      <  2:2087803/150‑1 (MQ=255)
                                                                            tatTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCta                                                                     <  1:2429152/150‑1 (MQ=255)
                                                                               tCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATAc                                                                  <  2:227297/150‑1 (MQ=255)
                                                                                    gTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAg                                                             >  1:805691/1‑150 (MQ=255)
                                                                                     tGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGc                                                            <  2:308399/150‑1 (MQ=255)
                                                                                            aaTGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCa                                                     >  2:3553207/1‑150 (MQ=255)
                                                                                            aaTGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCa                                                     <  1:3462895/150‑1 (MQ=255)
                                                                                            aaTGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCa                                                     >  1:1425923/1‑150 (MQ=255)
                                                                                              tGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATa                                                   <  1:1369691/150‑1 (MQ=255)
                                                                                              tGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATa                                                   <  2:3903511/150‑1 (MQ=255)
                                                                                              tGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATa                                                   <  2:2747205/150‑1 (MQ=255)
                                                                                                tCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATaaa                                                 >  2:3957752/1‑150 (MQ=255)
                                                                                                tCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATaaa                                                 <  1:1880945/150‑1 (MQ=255)
                                                                                                 cAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAAc                                                <  1:2301467/150‑1 (MQ=255)
                                                                                                  aGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACa                                               <  1:968002/150‑1 (MQ=255)
                                                                                                     aCCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAgg                                            <  2:3726675/150‑1 (MQ=255)
                                                                                                      ccACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGa                                           >  2:1362945/1‑150 (MQ=255)
                                                                                                        aCGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGAtt                                         <  1:3275145/150‑1 (MQ=255)
                                                                                                            cTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCt                                     <  1:29593/150‑1 (MQ=255)
                                                                                                             tGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCtt                                    <  1:926720/150‑1 (MQ=255)
                                                                                                             tGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCtt                                    <  2:320706/150‑1 (MQ=255)
                                                                                                              gAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCttt                                   <  1:746334/150‑1 (MQ=255)
                                                                                                               aaaTGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTg                                  <  1:203083/150‑1 (MQ=255)
                                                                                                                   gAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTAt                              <  2:1757718/150‑1 (MQ=255)
                                                                                                                    aaTTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATg                             <  2:1015190/150‑1 (MQ=255)
                                                                                                                         accaaccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTa                        <  1:4331515/150‑1 (MQ=255)
                                                                                                                            aaccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTaaaa                     >  2:3794169/1‑150 (MQ=255)
                                                                                                                            aaccaTCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTaaaa                     <  1:3211835/150‑1 (MQ=255)
                                                                                                                             accaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAAt                    <  1:3259494/150‑1 (MQ=255)
                                                                                                                              ccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATa                   <  2:3386129/150‑1 (MQ=255)
                                                                                                                                aTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAg                 <  1:3867142/150‑1 (MQ=255)
                                                                                                                                  ccGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGAc               <  1:3981345/150‑1 (MQ=255)
                                                                                                                                  ccGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGAc               <  2:629643/150‑1 (MQ=255)
                                                                                                                                    ggTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgt             <  2:724623/150‑1 (MQ=255)
                                                                                                                                    ggTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgt             >  1:1655403/1‑150 (MQ=255)
                                                                                                                                     gTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgtt            <  2:3343051/150‑1 (MQ=255)
                                                                                                                                     gTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgtt            <  1:1645820/150‑1 (MQ=255)
                                                                                                                                      tGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCGATGGGTTAAAATAAGACgttg           >  1:4026820/1‑150 (MQ=255)
                                                                                                                                        aCTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgttgtt         <  1:194147/150‑1 (MQ=255)
                                                                                                                                         cTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGttgttg        <  1:416308/150‑1 (MQ=255)
                                                                                                                                         cTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGttgttg        <  1:3014748/150‑1 (MQ=255)
                                                                                                                                         cTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGttgttg        <  2:1516914/150‑1 (MQ=255)
                                                                                                                                            gCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGaaa     <  2:235388/150‑1 (MQ=255)
                                                                                                                                            gCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGaaa     >  2:1048988/1‑150 (MQ=255)
                                                                                                                                               gACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGAAAATg  >  2:4236045/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGCAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGAAAATG  >  CP009273/4062278‑4062570

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