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,9070CA98.4% 603.8 / ‑0.4 186E298* (GAA→TAA) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base A (95/88);  total (98/88)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 187 total aligned reads displayed.

TGCTGTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACC  >  CP009273/2808762‑2809051
                                                                                                                                                 |                                                                                                                                                
tgctgTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGt                                                                                                                                                <  2:959808/148‑1 (MQ=255)
tgctgTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGt                                                                                                                                                >  2:1302035/1‑148 (MQ=255)
     ttCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCa                                                                                                                                           >  1:1031643/1‑148 (MQ=255)
       cATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAAt                                                                                                                                         <  1:791057/148‑1 (MQ=255)
         tCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCt                                                                                                                                       >  2:797833/1‑148 (MQ=255)
          cTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTg                                                                                                                                      >  2:2419562/1‑148 (MQ=255)
          cTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTg                                                                                                                                      >  1:2192002/1‑148 (MQ=255)
             cacCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCaa                                                                                                                                   >  2:2220007/1‑148 (MQ=255)
               ccAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAaca                                                                                                                                 <  2:1911192/148‑1 (MQ=255)
                   tCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGtt                                                                                                                             >  1:1080917/1‑148 (MQ=255)
                   tCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGtt                                                                                                                             >  1:162263/1‑148 (MQ=255)
                    cGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTg                                                                                                                            <  2:763922/148‑1 (MQ=255)
                    cGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTg                                                                                                                            >  1:3458560/1‑148 (MQ=255)
                       cGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTg                                                                                                                         >  1:1393474/1‑148 (MQ=255)
                       cGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTg                                                                                                                         <  2:2570211/148‑1 (MQ=255)
                        gagaACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGt                                                                                                                        <  1:1894075/148‑1 (MQ=255)
                          gaACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTg                                                                                                                      >  2:2329509/1‑148 (MQ=255)
                           aaCGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGa                                                                                                                     <  1:726045/148‑1 (MQ=255)
                            aCGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGAt                                                                                                                    <  2:1165633/148‑1 (MQ=255)
                            aCGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGAt                                                                                                                    <  2:567423/148‑1 (MQ=255)
                             cGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGAtt                                                                                                                   <  2:575484/148‑1 (MQ=255)
                               ggTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTg                                                                                                                 >  1:2441360/1‑148 (MQ=255)
                               ggTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTg                                                                                                                 >  1:1446572/1‑148 (MQ=255)
                               ggTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTg                                                                                                                 <  1:550083/148‑1 (MQ=255)
                               ggTTCATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTg                                                                                                                 >  2:219867/1‑148 (MQ=255)
                                     tGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCct                                                                                                           >  2:1463890/1‑148 (MQ=255)
                                       tCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCctct                                                                                                         >  1:2730658/1‑148 (MQ=255)
                                           gCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTAcc                                                                                                     >  2:392636/1‑148 (MQ=255)
                                           gCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTAcc                                                                                                     >  2:74448/1‑148 (MQ=255)
                                             gAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGt                                                                                                   <  2:3431146/148‑1 (MQ=255)
                                              aaCGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTc                                                                                                  <  2:113754/148‑1 (MQ=255)
                                               aCGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCt                                                                                                 >  1:1561753/1‑148 (MQ=255)
                                               aCGCACTTCAAGATAGTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCt                                                                                                 <  2:2368314/148‑1 (MQ=255)
                                                cGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCtt                                                                                                <  2:2183340/148‑1 (MQ=255)
                                                cGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCtt                                                                                                >  1:2832079/1‑148 (MQ=255)
                                                  cACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTtc                                                                                              >  2:684035/1‑148 (MQ=255)
                                                   tcTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTtct                                                                                             <  1:3345834/147‑1 (MQ=255)
                                                   aCTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTtct                                                                                             <  1:3345661/148‑1 (MQ=255)
                                                       cAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAAt                                                                                         >  2:474926/1‑148 (MQ=255)
                                                        aaTATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATa                                                                                        <  1:2373528/148‑1 (MQ=255)
                                                        aaTATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATa                                                                                        <  2:3544308/148‑1 (MQ=255)
                                                         atatatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAc                                                                                       >  2:3212210/1‑148 (MQ=255)
                                                          tatatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAcc                                                                                      >  2:2525261/1‑148 (MQ=255)
                                                          tatatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAcc                                                                                      >  2:2977296/1‑148 (MQ=255)
                                                          tatatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAcc                                                                                      >  1:426783/1‑148 (MQ=255)
                                                             atTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAAt                                                                                   <  2:1174555/148‑1 (MQ=255)
                                                                cAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATATTCAACACGTTGCTGTTGATTTGCAGCCTCATACCGTCATTCTCAATACCAATCtt                                                                                >  1:2794527/1‑148 (MQ=255)
                                                                  aTGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTcg                                                                              >  1:2233775/1‑148 (MQ=255)
                                                                  aTGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAAAACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTcg                                                                              >  1:2398899/1‑148 (MQ=255)
                                                                   tGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTcgc                                                                             >  1:2159103/1‑148 (MQ=255)
                                                                   tGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTcgc                                                                             >  1:980433/1‑148 (MQ=255)
                                                                        ccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTAct                                                                        <  1:3163/148‑1 (MQ=255)
                                                                        ccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTAct                                                                        <  1:940407/148‑1 (MQ=255)
                                                                            gTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTc                                                                    >  2:3300484/1‑148 (MQ=255)
                                                                            gTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTc                                                                    >  1:1853966/1‑148 (MQ=255)
                                                                            gTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACAGTCTTTCTCAATACCAATCTTCGCGTACTCTTc                                                                    >  2:3300501/1‑148 (MQ=255)
                                                                               aCAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGa                                                                 <  2:2965940/148‑1 (MQ=255)
                                                                               aCAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGa                                                                 >  1:3039662/1‑148 (MQ=255)
                                                                                cAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGAt                                                                >  2:2084283/1‑148 (MQ=255)
                                                                                  gcgcATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATgg                                                              >  1:693840/1‑148 (MQ=255)
                                                                                  gcgcATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGGACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATgg                                                              >  1:712110/1‑148 (MQ=255)
                                                                                   cgcATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGTAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGc                                                             >  1:2295419/1‑148 (MQ=255)
                                                                                    gcATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCg                                                            <  1:308588/148‑1 (MQ=255)
                                                                                     cATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGt                                                           >  1:169871/1‑148 (MQ=255)
                                                                                     cATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGt                                                           >  1:169307/1‑148 (MQ=255)
                                                                                       tCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGttt                                                         <  2:1976161/148‑1 (MQ=255)
                                                                                          gAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGa                                                      <  2:162263/148‑1 (MQ=255)
                                                                                          gAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGa                                                      <  2:3051106/148‑1 (MQ=255)
                                                                                           aaGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGAt                                                     >  2:943194/1‑148 (MQ=255)
                                                                                             ggCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTg                                                   <  2:3292448/148‑1 (MQ=255)
                                                                                               cGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGcc                                                 <  2:1273258/148‑1 (MQ=255)
                                                                                                      ccgcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTaa                                          <  1:1436978/148‑1 (MQ=255)
                                                                                                         cTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGc                                       <  1:2419562/148‑1 (MQ=255)
                                                                                                         cTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGc                                       >  2:2558219/1‑148 (MQ=255)
                                                                                                         cTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGc                                       <  1:2711911/148‑1 (MQ=255)
                                                                                                          tGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGcc                                      >  1:2035936/1‑148 (MQ=255)
                                                                                                          tGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGcc                                      <  1:2589832/148‑1 (MQ=255)
                                                                                                          tGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGcc                                      <  1:3199029/148‑1 (MQ=255)
                                                                                                             gggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGc                                   <  2:61462/148‑1 (MQ=255)
                                                                                                             gggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGc                                   <  1:277266/148‑1 (MQ=255)
                                                                                                              ggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCt                                  >  1:2082125/1‑148 (MQ=255)
                                                                                                              ggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCt                                  >  1:559017/1‑148 (MQ=255)
                                                                                                                   cgcgTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTa                             >  1:1126032/1‑148 (MQ=255)
                                                                                                                        tttGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGt                        >  2:2747416/1‑148 (MQ=255)
                                                                                                                        tttGGACGAATCGGCGCGTACAGTTAGTTTTCAAACTGCAACACCTTGCTGTTGATTTGCAGCCTCTTACCGTCTATATCAATACCAATCTTCGCCTACTCTTCACAAGGCCTTTTGATTGCCTGTGTAATGCACGCTACGTACTCGt                        >  1:2150165/1‑148 (MQ=255)
                                                                                                                         ttGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTa                       <  1:2679891/148‑1 (MQ=255)
                                                                                                                            gACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAg                    <  2:144844/148‑1 (MQ=255)
                                                                                                                            gACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAg                    >  1:871369/1‑148 (MQ=255)
                                                                                                                             aCGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                   >  1:1999073/1‑148 (MQ=255)
                                                                                                                                aaTCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCg                <  1:3469779/148‑1 (MQ=255)
                                                                                                                                aaTCGGCGCGTACAGTTAGTTTTCAATCTGCAACAAGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCg                <  2:643562/148‑1 (MQ=255)
                                                                                                                                   cGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTg             <  1:3212210/148‑1 (MQ=255)
                                                                                                                                   cGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTg             <  2:2049640/148‑1 (MQ=255)
                                                                                                                                     gcgcGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGaa           <  1:757404/148‑1 (MQ=255)
                                                                                                                                      cgcgTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAg          <  2:1653096/148‑1 (MQ=255)
                                                                                                                                      cgagTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAg          <  2:2423644/145‑1 (MQ=255)
                                                                                                                                       gcgTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAgg         <  1:219867/148‑1 (MQ=255)
                                                                                                                                         gTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAATGTa       <  2:1393474/148‑1 (MQ=255)
                                                                                                                                          tACAGTTAGTTGTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa      <  2:1851433/148‑1 (MQ=255)
                                                                                                                                           aCAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt     >  1:1233732/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGCTGTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACC  >  CP009273/2808762‑2809051

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