Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,808,907 C→A 100% E298* (GAA→TAA)  gshA ← glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,808,9070CA99.0% 635.5 / ‑1.6 196E298* (GAA→TAA) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (96/98);  total (98/98)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-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 201 total aligned reads displayed.

TGCTGTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATAC  >  CP009273/2808762‑2809050
                                                                                                                                                 |                                                                                                                                               
tgctgTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGt                                                                                                                                               >  1:2397804/1‑148 (MQ=255)
 gctgTGCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACagttagtt                                                                                                                                              >  1:2531433/1‑148 (MQ=255)
  ctgTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGttt                                                                                                                                             <  1:2383596/148‑1 (MQ=255)
    gTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTc                                                                                                                                           >  2:2342686/1‑148 (MQ=255)
    gTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTc                                                                                                                                           >  1:2022162/1‑148 (MQ=255)
     ttCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCa                                                                                                                                          <  2:2389232/148‑1 (MQ=255)
      tCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCaa                                                                                                                                         >  2:2733140/1‑148 (MQ=255)
       cATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGttt                                                                                                                                             >  2:450988/1‑143 (MQ=255)
         tCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCt                                                                                                                                      >  1:1736383/1‑148 (MQ=255)
          cTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTg                                                                                                                                     >  2:1356485/1‑148 (MQ=255)
          cTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTg                                                                                                                                     >  2:2281349/1‑148 (MQ=255)
               ccAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAaca                                                                                                                                >  2:1127966/1‑148 (MQ=255)
                   tCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGtt                                                                                                                            >  1:3375020/1‑148 (MQ=255)
                    cGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTg                                                                                                                           >  2:225870/1‑148 (MQ=255)
                    cGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTg                                                                                                                           <  2:1951600/148‑1 (MQ=255)
                    cGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTg                                                                                                                           <  2:120615/148‑1 (MQ=255)
                    cGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTg                                                                                                                           <  2:119665/148‑1 (MQ=255)
                    cGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTATTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTg                                                                                                                           >  2:227788/1‑148 (MQ=255)
                        gagaACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGt                                                                                                                       <  2:1272773/148‑1 (MQ=255)
                        gagaACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGt                                                                                                                       <  2:1474799/148‑1 (MQ=255)
                          gaACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTg                                                                                                                     <  2:775068/148‑1 (MQ=255)
                            aCGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGAt                                                                                                                   <  2:3127915/148‑1 (MQ=255)
                             cGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGAtt                                                                                                                  >  2:372329/1‑148 (MQ=255)
                              gggTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGAttt                                                                                                                 <  2:2063367/148‑1 (MQ=255)
                                         cAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTa                                                                                                      <  1:1904697/148‑1 (MQ=255)
                                         cAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTa                                                                                                      <  2:2547105/148‑1 (MQ=255)
                                           gCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTAcc                                                                                                    >  2:536547/1‑148 (MQ=255)
                                             gAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGt                                                                                                  <  1:1349897/148‑1 (MQ=255)
                                             gAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGt                                                                                                  <  1:3179819/148‑1 (MQ=255)
                                               aCGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCGGCCTCTTACCGTCt                                                                                                <  1:1190134/148‑1 (MQ=255)
                                               aCGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCt                                                                                                >  2:2371006/1‑148 (MQ=255)
                                               aCGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCt                                                                                                >  1:3042181/1‑148 (MQ=255)
                                               aCGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCt                                                                                                >  2:171626/1‑148 (MQ=255)
                                                   aCTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCTGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTtct                                                                                            <  2:558533/148‑1 (MQ=255)
                                                   aCTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTtct                                                                                            >  2:3494921/1‑148 (MQ=255)
                                                    cTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTtctc                                                                                           >  1:1536132/1‑148 (MQ=255)
                                                     ttCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCa                                                                                          >  1:127642/1‑148 (MQ=255)
                                                     ttCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCa                                                                                          >  2:1248400/1‑148 (MQ=255)
                                                         atatatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAc                                                                                      >  2:1914730/1‑148 (MQ=255)
                                                         atatatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAc                                                                                      >  2:18098/1‑148 (MQ=255)
                                                            tatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCaa                                                                                   >  1:3538533/1‑148 (MQ=255)
                                                             atTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAAt                                                                                  <  2:259808/148‑1 (MQ=255)
                                                                cAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCtt                                                                               <  1:1621336/148‑1 (MQ=255)
                                                                cAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCtt                                                                               >  1:3355130/1‑148 (MQ=255)
                                                                  aTGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTa                                                                                                      >  1:1970323/1‑123 (MQ=255)
                                                                  aTGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTa                                                                                                      >  1:170185/1‑123 (MQ=255)
                                                                  aTGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTcg                                                                             >  2:2158819/1‑148 (MQ=255)
                                                                        ccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAAACGGCGCGGACACTTAGTTTTCAAACAGCAACACGTTGCTGTTGATTTGCATCCTCTTACCCTCATTCTCAAAACCAATCTTCGCGTAct                                                                       >  1:719026/1‑148 (MQ=255)
                                                                         cACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTActc                                                                      <  1:659387/148‑1 (MQ=255)
                                                                           cGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCtt                                                                    <  2:2564558/148‑1 (MQ=255)
                                                                             tAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTcc                                                                  >  1:3460759/1‑148 (MQ=255)
                                                                             tAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTcc                                                                  >  2:407392/1‑148 (MQ=255)
                                                                              aaCAGCGCATCAGAAGGCGACTCGCCGCTGCGTGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCg                                                                 <  2:355123/148‑1 (MQ=255)
                                                                              aaCAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCg                                                                 <  1:2281349/148‑1 (MQ=255)
                                                                               aCAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGa                                                                <  2:3172471/148‑1 (MQ=255)
                                                                                cAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGAt                                                               >  2:2167821/1‑148 (MQ=255)
                                                                                  gcgcATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATgg                                                             <  2:1163483/148‑1 (MQ=255)
                                                                                     cATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATATGCAGCCTCTTACCGTCTTTCTCAATAACAATCTTCGCGTACTCTTCCGATGGCGt                                                          >  2:3381577/1‑148 (MQ=255)
                                                                                      aTCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGtt                                                         >  1:3122743/1‑148 (MQ=255)
                                                                                      aTCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGtt                                                         <  1:262628/148‑1 (MQ=255)
                                                                                       tCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAACAATCTTCGCGTACTCTTCCGATGGCGttt                                                        >  1:1583198/1‑148 (MQ=255)
                                                                                        cAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGtttt                                                       >  2:240077/1‑148 (MQ=255)
                                                                                          gAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGa                                                     >  1:1037008/1‑148 (MQ=255)
                                                                                             ggCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTg                                                  <  2:2786763/148‑1 (MQ=255)
                                                                                             ggCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTg                                                  >  1:3325518/1‑148 (MQ=255)
                                                                                               cGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGcc                                                >  1:1972622/1‑148 (MQ=255)
                                                                                                 aCTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTg                                              >  1:459177/1‑148 (MQ=255)
                                                                                                   tCGCCGCTGCGGGTAACGCGTTTTGGGCGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGtt                                            <  1:1327276/148‑1 (MQ=255)
                                                                                                      ccgcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTaa                                         <  2:1657115/148‑1 (MQ=255)
                                                                                                      ccgcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTaa                                         <  1:1226257/148‑1 (MQ=255)
                                                                                                       cgcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAg                                        <  2:2531433/148‑1 (MQ=255)
                                                                                                         cTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGc                                      <  1:1291902/148‑1 (MQ=255)
                                                                                                               gTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAATGCCCGCTa                                >  1:3323020/1‑148 (MQ=255)
                                                                                                                tAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTAc                               >  2:901110/1‑148 (MQ=255)
                                                                                                                  aCGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGt                             <  1:2597023/148‑1 (MQ=255)
                                                                                                                  aCGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGt                             >  2:376018/1‑148 (MQ=255)
                                                                                                                  aCGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGt                             <  2:576259/148‑1 (MQ=255)
                                                                                                                     cgTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCATCCGATGGCGTTTTGATTGCCTGTTTAATGCCCGCTACGTACt                          >  1:886234/1‑148 (MQ=255)
                                                                                                                       ttttGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCg                        <  1:1120343/148‑1 (MQ=255)
                                                                                                                       ttttGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCg                        <  1:3445425/148‑1 (MQ=255)
                                                                                                                       ttttGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCg                        >  2:3055744/1‑148 (MQ=255)
                                                                                                                        tttGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGt                       >  1:2307090/1‑148 (MQ=255)
                                                                                                                           ggACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaaa                    >  2:101495/1‑148 (MQ=255)
                                                                                                                            gACGAATCGGCGCGTACAGTTAGTTTTCAATCTTCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAg                   <  2:719026/148‑1 (MQ=255)
                                                                                                                            gACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAg                   <  1:2733140/148‑1 (MQ=255)
                                                                                                                             aCGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                  <  1:227041/148‑1 (MQ=255)
                                                                                                                             aCGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                  <  1:1179322/148‑1 (MQ=255)
                                                                                                                                aaTCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCg               >  2:565465/1‑148 (MQ=255)
                                                                                                                                 aTCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGt              <  1:3559594/148‑1 (MQ=255)
                                                                                                                                  tCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGtt             >  1:2760855/1‑148 (MQ=255)
                                                                                                                                   cGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTg            <  1:2158819/148‑1 (MQ=255)
                                                                                                                                     gcgcGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGaa          >  1:1735894/1‑148 (MQ=255)
                                                                                                                                      cgcgTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAg         <  1:3152591/148‑1 (MQ=255)
                                                                                                                                         gTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTa      <  2:200472/148‑1 (MQ=255)
                                                                                                                                          tACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa     >  2:162857/1‑148 (MQ=255)
                                                                                                                                           aCAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt    <  1:2078959/148‑1 (MQ=255)
                                                                                                                                           aCAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt    <  1:2078911/148‑1 (MQ=255)
                                                                                                                                           aCAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt    >  1:3330951/1‑148 (MQ=255)
                                                                                                                                             agttagttTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATAc  >  1:2229599/1‑148 (MQ=255)
                                                                                                                                             agttagttTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATAc  >  1:1145785/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TGCTGTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATAC  >  CP009273/2808762‑2809050

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