Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 215,833 G→T 10.1% N83K (AAC→AAA tsaA ← tRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273215,8330GT10.1% 782.0 / 41.6 228N83K (AAC→AAAtsaAtRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (96/108);  new base T (13/10);  total (110/118)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 233 total aligned reads displayed.

TACCGGCGTACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAATGAGATGCAGTTC  >  CP009273/215686‑215979
                                                                                                                                                   |                                                                                                                                                  
tACCGGCGTACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCg                                                                                                                                                    <  2:1116249/148‑1 (MQ=255)
tACCGGCGTACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCg                                                                                                                                                    >  1:1387065/1‑148 (MQ=255)
 aCCGGCGTACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGt                                                                                                                                                   >  1:1561612/1‑148 (MQ=255)
   cGGCGTACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTa                                                                                                                                                 <  1:350577/148‑1 (MQ=255)
   cGGCGTACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTa                                                                                                                                                 <  2:2267491/148‑1 (MQ=255)
     gCGTACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTAcc                                                                                                                                               >  1:2370132/1‑148 (MQ=255)
     gCGTACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCAGTCTGGCGTTAcc                                                                                                                                               >  1:1768009/1‑148 (MQ=255)
      cGTACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCTTTAccg                                                                                                                                              >  1:2095006/1‑148 (MQ=255)
      cGTACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTAccg                                                                                                                                              >  2:2651888/1‑148 (MQ=255)
      cGTACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTAccg                                                                                                                                              >  2:1822666/1‑148 (MQ=255)
        tACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCCATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCTTTAccgcc                                                                                                                                            >  1:2020449/1‑148 (MQ=255)
         aCCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTAccgccg                                                                                                                                           <  1:253566/148‑1 (MQ=255)
         aCCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTAccgccg                                                                                                                                           >  1:1352668/1‑148 (MQ=255)
                 ccAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCgg                                                                                                                                   <  2:879769/148‑1 (MQ=255)
                 ccAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCgg                                                                                                                                   <  2:879786/148‑1 (MQ=255)
                    gATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTgg                                                                                                                                >  1:42988/1‑148 (MQ=255)
                        cAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGc                                                                                                                            <  1:2748889/148‑1 (MQ=255)
                              gCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTc                                                                                                                      <  1:773387/148‑1 (MQ=255)
                                    cTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCTTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGc                                                                                                                <  1:1822410/148‑1 (MQ=255)
                                       cAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACCTGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACAGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAg                                                                                                             >  2:1250773/1‑148 (MQ=255)
                                       cAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAg                                                                                                             >  1:3046220/1‑148 (MQ=255)
                                        aGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAgc                                                                                                            <  2:624787/148‑1 (MQ=255)
                                               cGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTc                                                                                                     <  1:1687671/148‑1 (MQ=255)
                                               cGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTc                                                                                                     <  2:1745767/148‑1 (MQ=255)
                                                  tGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCAt                                                                                                  >  1:1832166/1‑148 (MQ=255)
                                                  tGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCAt                                                                                                  >  1:1520689/1‑148 (MQ=255)
                                                    tCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCTTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCg                                                                                                <  1:2719295/148‑1 (MQ=255)
                                                    tCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCg                                                                                                <  1:1034136/148‑1 (MQ=255)
                                                    tCTTTATGGCAAACAACCTCTTTCAGCTCTACAAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGATCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCg                                                                                                >  2:475846/1‑148 (MQ=255)
                                                     cTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGt                                                                                               >  1:2227556/1‑148 (MQ=255)
                                                      tttATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGtt                                                                                              >  2:2859184/1‑148 (MQ=255)
                                                      tttATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGtt                                                                                              <  1:1555615/148‑1 (MQ=255)
                                                        tATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTtg                                                                                            <  2:1802016/148‑1 (MQ=255)
                                                          tGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTtgat                                                                                          <  1:1893390/148‑1 (MQ=255)
                                                          tGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTtgat                                                                                          <  1:2819232/148‑1 (MQ=255)
                                                           ggCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTtgatg                                                                                         >  1:2691224/1‑148 (MQ=255)
                                                             cAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGaa                                                                                       >  2:360601/1‑148 (MQ=255)
                                                                acaacCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGa                                                                                    <  2:3094103/148‑1 (MQ=255)
                                                                acaacCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGa                                                                                    >  2:2361876/1‑148 (MQ=255)
                                                                acaacCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGa                                                                                    <  2:1287957/148‑1 (MQ=255)
                                                                 caacCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCTTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGAc                                                                                   >  2:494266/1‑148 (MQ=255)
                                                                     ctctTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGaaa                                                                               <  1:2864878/148‑1 (MQ=255)
                                                                       ctTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAg                                                                             >  2:2603235/1‑148 (MQ=255)
                                                                       ctTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAg                                                                             >  2:252425/1‑148 (MQ=255)
                                                                         ttCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGa                                                                           >  2:892782/1‑148 (MQ=255)
                                                                          tCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGAt                                                                          <  1:1360302/148‑1 (MQ=255)
                                                                           cAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATc                                                                         <  1:765597/148‑1 (MQ=255)
                                                                            aGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATcc                                                                        <  2:1520961/148‑1 (MQ=255)
                                                                             gCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCa                                                                       >  2:706808/1‑148 (MQ=255)
                                                                             gCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCa                                                                       <  1:1268850/148‑1 (MQ=255)
                                                                                ctACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATaa                                                                    <  2:1767297/148‑1 (MQ=255)
                                                                                ctACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATaa                                                                    <  1:2618498/148‑1 (MQ=255)
                                                                                ctACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATaa                                                                    <  1:2215646/148‑1 (MQ=255)
                                                                                  aCCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAAt                                                                  <  2:2895246/148‑1 (MQ=255)
                                                                                  aCCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAAt                                                                  <  2:3220965/148‑1 (MQ=255)
                                                                                  aCCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAAt                                                                  <  2:1352658/148‑1 (MQ=255)
                                                                                    cAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATgg                                                                >  1:476756/1‑148 (MQ=255)
                                                                                    cAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATgg                                                                <  1:1216300/148‑1 (MQ=255)
                                                                                    cAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATgg                                                                >  1:3199390/1‑148 (MQ=255)
                                                                                        gACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGa                                                            <  2:1660305/148‑1 (MQ=255)
                                                                                         aCAGGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGaa                                                           <  2:1712101/148‑1 (MQ=255)
                                                                                           aTGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGACGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACg                                                         <  1:474395/148‑1 (MQ=255)
                                                                                            tGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGc                                                        <  1:2700917/148‑1 (MQ=255)
                                                                                               cAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTc                                                     >  2:2662992/1‑148 (MQ=255)
                                                                                                 aTTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCa                                                   >  2:1370203/1‑148 (MQ=255)
                                                                                                 aTTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCa                                                   >  2:1887031/1‑148 (MQ=255)
                                                                                                   tGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAgg                                                 >  1:899236/1‑148 (MQ=255)
                                                                                                       ttAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCcga                                             <  1:1822666/148‑2 (MQ=255)
                                                                                                       ttAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCcga                                             <  1:2651888/148‑2 (MQ=255)
                                                                                                           ggCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAAc                                         >  1:166442/1‑148 (MQ=255)
                                                                                                              ggAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGc                                      <  2:1387009/148‑1 (MQ=255)
                                                                                                               gAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCg                                     >  1:2946282/1‑148 (MQ=255)
                                                                                                               gAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCg                                     >  1:2946099/1‑148 (MQ=255)
                                                                                                                aaaGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGt                                    <  1:596098/148‑1 (MQ=255)
                                                                                                                    tAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGc                                >  1:2414908/1‑148 (MQ=255)
                                                                                                                       agCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTg                             >  2:3030511/1‑148 (MQ=255)
                                                                                                                         cGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGt                           >  1:251768/1‑148 (MQ=255)
                                                                                                                           tgtTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTg                         <  1:252428/148‑1 (MQ=255)
                                                                                                                           tgtTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTg                         >  1:3130412/1‑148 (MQ=255)
                                                                                                                           tgtTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAAGGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTg                         <  1:252425/148‑1 (MQ=255)
                                                                                                                              tGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAg                      <  1:2904544/148‑1 (MQ=255)
                                                                                                                                cGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAggg                    <  2:1031673/148‑1 (MQ=255)
                                                                                                                                cGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAggg                    <  2:24576/148‑1 (MQ=255)
                                                                                                                                cGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAggg                    <  2:1275179/148‑1 (MQ=255)
                                                                                                                                   aaaCCCCCATTCTGGCTTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGc                 >  2:2973657/1‑148 (MQ=255)
                                                                                                                                    aaCCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCa                <  2:2993538/148‑1 (MQ=255)
                                                                                                                                    aaCACCCATTCTGGCATTACCGCCAAGACGAGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGTATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGcc                >  1:1226771/1‑147 (MQ=255)
                                                                                                                                     aCCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCaa               >  1:1887482/1‑148 (MQ=255)
                                                                                                                                      cccccATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAAt              >  2:602509/1‑148 (MQ=255)
                                                                                                                                       ccccATTCTGGCTTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAATg             <  1:26709/148‑1 (MQ=255)
                                                                                                                                       ccccATTCTGGCTTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAATg             >  1:3172850/1‑148 (MQ=255)
                                                                                                                                       ccccATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAATg             <  2:2750574/148‑1 (MQ=255)
                                                                                                                                         ccATTCTGGCTTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAATgag           <  2:2095006/148‑1 (MQ=255)
                                                                                                                                         ccATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAATgag           >  1:1744018/1‑148 (MQ=255)
                                                                                                                                          cATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAATgaga          <  1:2241620/148‑1 (MQ=255)
                                                                                                                                            ttCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAATGAGATg        <  2:42988/148‑1 (MQ=255)
                                                                                                                                             tCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAATGAGATGc       <  2:2152433/148‑1 (MQ=255)
                                                                                                                                             tCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAAGGAGCAATGAGATGc       >  1:2425433/1‑148 (MQ=255)
                                                                                                                                               tGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAATGAGATGCAg     <  1:693180/148‑1 (MQ=255)
                                                                                                                                               tGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAATGAGATGCAg     <  2:1446648/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TACCGGCGTACCATCGACCAGATCCAGACTGCCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGGAGCAATGAGATGCAGTTC  >  CP009273/215686‑215979

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: