Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,572,520 A→T 91.8% Y116N (TAC→AAC)  hsdM ← DNA methyltransferase M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,572,5200AT91.8% 488.4 / 28.5 158Y116N (TAC→AAC) hsdMDNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base A (3/10);  new base T (80/65);  total (83/75)
Fisher's exact test for biased strand distribution p-value = 3.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 161 total aligned reads displayed.

TTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGAACTGCAA  >  CP009273/4572375‑4572664
                                                                                                                                                 |                                                                                                                                                
ttCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTcc                                                                                                                                                >  1:404730/1‑148 (MQ=255)
ttCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTAcc                                                                                                                                                <  1:1269847/148‑1 (MQ=255)
 tCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCa                                                                                                                                               <  1:627278/148‑1 (MQ=255)
   agcagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGt                                                                                                                                             <  1:1128427/148‑1 (MQ=255)
    gcagATGAAGAATGGTTTTCATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTc                                                                                                                                            <  2:2014552/148‑1 (MQ=255)
       gATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTtccagtccag                                                                                                                                         <  2:2493444/148‑1 (MQ=255)
       gATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTtccagtccag                                                                                                                                         <  2:2143243/148‑1 (MQ=255)
         tGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCg                                                                                                                                       >  1:2210819/1‑148 (MQ=255)
         tGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCg                                                                                                                                       >  1:2210809/1‑148 (MQ=255)
         tGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCg                                                                                                                                       >  1:1644458/1‑148 (MQ=255)
            ataatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAAt                                                                                                                                    <  2:1780725/148‑1 (MQ=255)
              aatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATcc                                                                                                                                  >  2:2322085/1‑148 (MQ=255)
              aatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATcc                                                                                                                                  <  1:1642138/148‑1 (MQ=255)
                tGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCat                                                                                                                                >  2:1873701/1‑148 (MQ=255)
                 ggTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCata                                                                                                                               >  2:222279/1‑148 (MQ=255)
                 ggTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCata                                                                                                                               >  2:222357/1‑148 (MQ=255)
                  gTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCatat                                                                                                                              >  1:2533130/1‑148 (MQ=255)
                    tttAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTg                                                                                                                            <  1:2076706/148‑1 (MQ=255)
                        aTCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGa                                                                                                                        <  2:1853797/148‑1 (MQ=255)
                          cAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGAcc                                                                                                                      >  2:395329/1‑148 (MQ=255)
                          cAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGAcc                                                                                                                      >  2:1402655/1‑148 (MQ=255)
                               gACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGc                                                                                                                 >  2:790779/1‑148 (MQ=255)
                                 cgcgGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCgg                                                                                                               <  1:997603/148‑1 (MQ=255)
                                     gggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAt                                                                                                           >  1:1243688/1‑148 (MQ=255)
                                     gggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAt                                                                                                           >  1:2892476/1‑148 (MQ=255)
                                      ggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAtt                                                                                                          <  2:3193732/148‑1 (MQ=255)
                                        tGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTg                                                                                                        >  1:2446545/1‑148 (MQ=255)
                                        tGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTg                                                                                                        >  2:2059256/1‑148 (MQ=255)
                                        tGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTg                                                                                                        <  2:101103/148‑1 (MQ=255)
                                          aaGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGtt                                                                                                      <  1:1372301/148‑1 (MQ=255)
                                           aGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGttt                                                                                                     >  2:85011/1‑148 (MQ=255)
                                           aGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGttt                                                                                                     <  2:735250/148‑1 (MQ=255)
                                            gTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTc                                                                                                    <  1:2451340/148‑1 (MQ=255)
                                             tACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCg                                                                                                   >  1:1971200/1‑148 (MQ=255)
                                             tACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCg                                                                                                   >  2:160208/1‑148 (MQ=255)
                                              aCTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCgg                                                                                                  <  1:2985313/148‑1 (MQ=255)
                                              aCTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCgg                                                                                                  >  2:2680387/1‑148 (MQ=255)
                                              aCTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCgg                                                                                                  >  2:1028705/1‑148 (MQ=255)
                                                  gCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCg                                                                                              <  1:2374766/148‑1 (MQ=255)
                                                  gCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCg                                                                                              >  2:2038498/1‑148 (MQ=255)
                                                  gCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCg                                                                                              <  1:2938468/148‑1 (MQ=255)
                                                      gCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGa                                                                                          <  1:2652345/148‑1 (MQ=255)
                                                      gCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGa                                                                                          <  1:186424/148‑1 (MQ=255)
                                                        aCCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtg                                                                                        <  1:2362886/148‑1 (MQ=255)
                                                          cAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtggt                                                                                      <  2:2410753/148‑1 (MQ=255)
                                                            gACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtggtgg                                                                                    >  1:2569057/1‑148 (MQ=255)
                                                              cTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTa                                                                                  <  2:865329/148‑1 (MQ=255)
                                                              cTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTa                                                                                  >  2:437243/1‑148 (MQ=255)
                                                                  gTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTaa                                                                              >  1:2216564/1‑148 (MQ=255)
                                                                  gTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTaa                                                                              >  1:1849913/1‑148 (MQ=255)
                                                                      cATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACAtt                                                                          <  1:1537612/148‑1 (MQ=255)
                                                                      cATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACAtt                                                                          <  1:1537625/148‑1 (MQ=255)
                                                                       aTTCGCGGTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTa                                                                         <  2:2600589/148‑1 (MQ=255)
                                                                        ttCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGACTCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCGCGGTGATGGTGGTACTAACATTAt                                                                        >  2:1979164/1‑148 (MQ=255)
                                                                        ttCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTAt                                                                        >  2:311156/1‑148 (MQ=255)
                                                                        ttCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTAt                                                                        >  2:311234/1‑148 (MQ=255)
                                                                         tCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATg                                                                       >  2:2922321/1‑148 (MQ=255)
                                                                             gTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaaa                                                                   <  2:537102/148‑1 (MQ=255)
                                                                              ttcttcTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaaaa                                                                  >  1:3206615/1‑148 (MQ=255)
                                                                              ttcttcTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaaaa                                                                  >  2:2054006/1‑148 (MQ=255)
                                                                               tcttctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAAc                                                                 >  2:433653/1‑148 (MQ=255)
                                                                                cttctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACt                                                                >  1:1986701/1‑148 (MQ=255)
                                                                                  tctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGc                                                              >  1:3123010/1‑148 (MQ=255)
                                                                                  tctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGc                                                              >  2:2788132/1‑148 (MQ=255)
                                                                                  tctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGc                                                              <  1:837290/148‑1 (MQ=255)
                                                                                    tGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCt                                                            <  2:189551/148‑1 (MQ=255)
                                                                                    tGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCt                                                            >  2:2758620/1‑148 (MQ=255)
                                                                                       aaCAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTa                                                         <  2:404730/148‑1 (MQ=255)
                                                                                          aGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCa                                                      >  2:521252/1‑148 (MQ=255)
                                                                                            cccTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGc                                                    >  2:2973121/1‑148 (MQ=255)
                                                                                            cccTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGc                                                    >  2:1325433/1‑148 (MQ=255)
                                                                                              cTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtt                                                  <  1:2620361/148‑1 (MQ=255)
                                                                                              cTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtt                                                  <  1:820351/148‑1 (MQ=255)
                                                                                                tCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTAGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtttt                                                >  1:1167805/1‑148 (MQ=255)
                                                                                                  gTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTg                                              <  2:1283744/148‑1 (MQ=255)
                                                                                                  gTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTg                                              <  2:601042/148‑1 (MQ=255)
                                                                                                     cATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCa                                           <  1:1584449/148‑1 (MQ=255)
                                                                                                     cATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCa                                           <  1:2143500/148‑1 (MQ=255)
                                                                                                       tatCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATc                                         <  2:1663692/148‑1 (MQ=255)
                                                                                                         tCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTNTTCCAGTCCAGCGAATCCATATTGNTGNCCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTNTGAAAAACTGCCTGTACCAGCTTTTTGTCATCtt                                       >  1:3278948/1‑148 (MQ=255)
                                                                                                         tCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCtt                                       <  2:2210819/148‑1 (MQ=255)
                                                                                                         tCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCtt                                       <  2:2210809/148‑1 (MQ=255)
                                                                                                          cGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTc                                      >  2:2907697/1‑148 (MQ=255)
                                                                                                            ccGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGc                                    >  2:1589081/1‑148 (MQ=255)
                                                                                                            ccGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTGATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTGTTTGTCATCTTCGc                                    >  2:1705553/1‑148 (MQ=255)
                                                                                                               aaGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTa                                 >  1:2675132/1‑148 (MQ=255)
                                                                                                                 gTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTaaa                               >  2:986308/1‑148 (MQ=255)
                                                                                                                 gTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTaaa                               >  2:1705000/1‑148 (MQ=255)
                                                                                                                 gTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTaaa                               <  2:1816245/148‑1 (MQ=255)
                                                                                                                  tCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAAt                              <  1:1464929/148‑1 (MQ=255)
                                                                                                                     tCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCa                           >  1:1983436/1‑148 (MQ=255)
                                                                                                                       gcgcgACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACg                         <  2:2313298/148‑1 (MQ=255)
                                                                                                                                aCCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAttttt                >  1:3244183/1‑148 (MQ=255)
                                                                                                                                aCCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAttttt                <  1:2059256/148‑1 (MQ=255)
                                                                                                                                    tGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACAGCCTGTACCAGCTTTTTGTCATCTTAGCCTAAATGCACGAGCATTTTTCGGt            >  1:1143550/1‑148 (MQ=255)
                                                                                                                                            gttgttCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgc    >  1:963257/1‑148 (MQ=255)
                                                                                                                                             ttGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgca   <  1:81946/148‑1 (MQ=255)
                                                                                                                                              tGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgcaa  <  1:84780/148‑1 (MQ=255)
                                                                                                                                              tGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgcaa  <  1:3125683/148‑1 (MQ=255)
                                                                                                                                              tGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgcaa  >  2:878263/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGAACTGCAA  >  CP009273/4572375‑4572664

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