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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,062,4240GT99.3% 997.5 / ‑1.8 279I71I (ATC→ATAyihQalpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (129/148);  total (129/150)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 279 total aligned reads displayed.

CGCAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGAAAATG  >  CP009273/4062278‑4062570
                                                                                                                                                  |                                                                                                                                                  
cgcAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATg                                                                                                                                                 <  2:4043603/150‑1 (MQ=255)
 gcAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATgg                                                                                                                                                <  1:1153603/150‑1 (MQ=255)
  cAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGc                                                                                                                                               >  1:747957/1‑150 (MQ=255)
  cAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGc                                                                                                                                               >  2:2366123/1‑150 (MQ=255)
   aGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCg                                                                                                                                              >  2:289473/1‑150 (MQ=255)
     ccAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTc                                                                                                                                            <  1:2559981/150‑1 (MQ=255)
      cAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCg                                                                                                                                           >  1:2003476/1‑150 (MQ=255)
      cAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCg                                                                                                                                           >  2:4237/1‑150 (MQ=255)
      cAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCg                                                                                                                                           >  2:3603609/1‑150 (MQ=255)
      cAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCg                                                                                                                                           <  1:2483421/150‑1 (MQ=255)
        gATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGt                                                                                                                                         >  1:2086369/1‑150 (MQ=255)
           aCGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAg                                                                                                                                      >  1:325495/1‑150 (MQ=255)
             ggTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAgcg                                                                                                                                    <  2:1832488/150‑1 (MQ=255)
              gTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAgcgc                                                                                                                                   <  1:507153/150‑1 (MQ=255)
                 gtgGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAAt                                                                                                                                <  1:1398878/150‑1 (MQ=255)
                     ttAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTtc                                                                                                                            <  2:2560267/150‑1 (MQ=255)
                         ggTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCt                                                                                                                        >  1:271110/1‑150 (MQ=255)
                             gTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAg                                                                                                                    <  1:3454530/150‑1 (MQ=255)
                              tCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGt                                                                                                                   <  2:723494/150‑1 (MQ=255)
                                 ttttGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTa                                                                                                                >  2:1076013/1‑150 (MQ=255)
                                      tAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCttt                                                                                                           <  2:870587/150‑1 (MQ=255)
                                       aGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTa                                                                                                          >  1:3064887/1‑150 (MQ=255)
                                       aGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTa                                                                                                          <  1:1050926/150‑1 (MQ=255)
                                       aGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGAGCTAATGTCAGAACAACGGCTGAAATGAATTAAACAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTa                                                                                                          >  2:827629/1‑150 (MQ=255)
                                         ttCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAAt                                                                                                        <  2:2092900/150‑1 (MQ=255)
                                          tCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATg                                                                                                       <  1:4111221/150‑1 (MQ=255)
                                            cAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCt                                                                                                     <  2:3698269/150‑1 (MQ=255)
                                            cAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCt                                                                                                     <  2:1044896/150‑1 (MQ=255)
                                            cAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCt                                                                                                     >  2:3833035/1‑150 (MQ=255)
                                             aGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTg                                                                                                    >  2:2693113/1‑150 (MQ=255)
                                             aGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTg                                                                                                    >  2:1098145/1‑150 (MQ=255)
                                                aaTAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGAGCAATTTTATCCTGTAGTTTATCTTTAATGCTGaaa                                                                                                 >  1:1638948/1‑150 (MQ=255)
                                                  tAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAAtt                                                                                               <  2:4074430/150‑1 (MQ=255)
                                                   aaaCGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCATCCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTa                                                                                              <  2:883412/150‑1 (MQ=255)
                                                   aaaCGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTa                                                                                              >  2:574268/1‑150 (MQ=255)
                                                    aaCGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTAc                                                                                             >  1:1539645/1‑150 (MQ=255)
                                                       gCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcgc                                                                                          >  1:17242/1‑150 (MQ=255)
                                                         cccTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCgg                                                                                        >  1:2700439/1‑150 (MQ=255)
                                                            tGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAAc                                                                                     >  2:790867/1‑150 (MQ=255)
                                                             gATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACa                                                                                    <  2:1811283/150‑1 (MQ=255)
                                                               tcgtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACata                                                                                  <  1:3586249/150‑1 (MQ=255)
                                                                cgtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACatat                                                                                 <  1:3135187/150‑1 (MQ=255)
                                                                cgtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACatat                                                                                 <  2:2801061/150‑1 (MQ=255)
                                                                 gtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATc                                                                                <  1:2167043/150‑1 (MQ=255)
                                                                 gtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATc                                                                                <  1:682614/150‑1 (MQ=255)
                                                                 gtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATc                                                                                <  1:2770871/150‑1 (MQ=255)
                                                                      cAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGatat                                                                           >  2:4009827/1‑150 (MQ=255)
                                                                       agagATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGTAATTACCGCGGAACATATCGATATc                                                                          >  2:3069516/1‑150 (MQ=255)
                                                                       agagATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGACCATATCGATATc                                                                          >  2:223464/1‑150 (MQ=255)
                                                                       agagATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATc                                                                          >  1:1363136/1‑150 (MQ=255)
                                                                        gagaTATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATcc                                                                         <  2:1127051/150‑1 (MQ=255)
                                                                        gagaTATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATcc                                                                         <  2:1537914/150‑1 (MQ=255)
                                                                          gaTATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGc                                                                       >  2:1788972/1‑150 (MQ=255)
                                                                           atatTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCt                                                                      <  2:1248309/150‑1 (MQ=255)
                                                                             atTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCtat                                                                    <  1:1068100/150‑1 (MQ=255)
                                                                             atTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCtat                                                                    >  2:2181187/1‑150 (MQ=255)
                                                                              ttCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCtata                                                                   <  1:3780400/150‑1 (MQ=255)
                                                                               tCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATAc                                                                  <  1:2529733/150‑1 (MQ=255)
                                                                               tCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGAGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATAc                                                                  <  1:4002766/150‑1 (MQ=255)
                                                                                 aGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCt                                                                >  1:1414928/1‑150 (MQ=255)
                                                                                     tGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGc                                                            >  2:1078223/1‑150 (MQ=255)
                                                                                         gcTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAAt                                                        >  1:397690/1‑150 (MQ=255)
                                                                                         gcTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAAt                                                        >  1:397651/1‑150 (MQ=255)
                                                                                          cTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATc                                                       >  2:438453/1‑150 (MQ=255)
                                                                                          cTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATc                                                       >  1:2873212/1‑150 (MQ=255)
                                                                                           tAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATcc                                                      >  2:1368665/1‑150 (MQ=255)
                                                                                            aaTGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCa                                                     <  1:2866229/150‑1 (MQ=255)
                                                                                            aaTGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCa                                                     <  1:1943759/150‑1 (MQ=255)
                                                                                            aaTGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCa                                                     <  2:1202072/150‑1 (MQ=255)
                                                                                                tCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATaaa                                                 >  2:2740325/1‑150 (MQ=255)
                                                                                                  aGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACa                                               <  1:1826904/150‑1 (MQ=255)
                                                                                                   gAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACaa                                              >  2:357903/1‑150 (MQ=255)
                                                                                                   gAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACaa                                              <  1:2848160/150‑1 (MQ=255)
                                                                                                   gAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACaa                                              <  2:2208272/150‑1 (MQ=255)
                                                                                                   gAACCACGGCTGAAATAAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACaa                                              >  2:3426125/1‑150 (MQ=255)
                                                                                                      ccACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGa                                           <  2:3471953/150‑1 (MQ=255)
                                                                                                         cGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTa                                        >  2:250016/1‑150 (MQ=255)
                                                                                                           gCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATc                                      <  1:1444657/150‑1 (MQ=255)
                                                                                                             tGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCtt                                    >  2:970911/1‑150 (MQ=255)
                                                                                                             tGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCtt                                    >  1:1252106/1‑150 (MQ=255)
                                                                                                             tGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCtt                                    >  1:2182165/1‑150 (MQ=255)
                                                                                                             tGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCtt                                    <  1:912136/150‑1 (MQ=255)
                                                                                                              gAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCttt                                   <  1:3476157/150‑1 (MQ=255)
                                                                                                                  tGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATGTCGATATCCGCTATACCTGAGCCACTCCATAAACAAGGAGTAGCTTTGcgt                               >  1:3274434/1‑148 (MQ=255)
                                                                                                                   gAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTAt                              <  2:1057269/150‑1 (MQ=255)
                                                                                                                      ttAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATggg                           <  1:429143/150‑1 (MQ=255)
                                                                                                                        aaccaaccATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGtt                         >  1:1013238/1‑150 (MQ=255)
                                                                                                                         accaaccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTa                        <  1:4167576/150‑1 (MQ=255)
                                                                                                                          ccaaccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTaa                       <  2:166367/150‑1 (MQ=255)
                                                                                                                            aaccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGCACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTaaaa                     >  2:3783589/1‑150 (MQ=255)
                                                                                                                            aaccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTaaaa                     >  2:3156742/1‑150 (MQ=255)
                                                                                                                             accaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAAt                    <  1:3349221/150‑1 (MQ=255)
                                                                                                                               caTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATaa                  >  2:3677056/1‑150 (MQ=255)
                                                                                                                                aTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAg                 >  2:2814989/1‑150 (MQ=255)
                                                                                                                                aTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAg                 >  2:1422052/1‑150 (MQ=255)
                                                                                                                                   cGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACg              <  1:729203/150‑1 (MQ=255)
                                                                                                                                    ggTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgt             >  1:2739556/1‑150 (MQ=255)
                                                                                                                                          tGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGa       <  1:1508674/150‑1 (MQ=255)
                                                                                                                                            gCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGaaa     <  2:1755181/150‑1 (MQ=255)
                                                                                                                                              tGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGAAAAt   >  2:1579484/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: