Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,572,520 A→T 100% 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,5200AT97.4% 626.1 / 4.6 192Y116N (TAC→AAC) hsdMDNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base T (91/96);  total (95/97)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
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 193 total aligned reads displayed.

TCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGAACTGCAA  >  CP009273/4572376‑4572664
                                                                                                                                                |                                                                                                                                                
tCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCa                                                                                                                                               >  2:3087840/1‑148 (MQ=255)
   gcagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTc                                                                                                                                            <  1:3093406/148‑1 (MQ=255)
   gcagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTc                                                                                                                                            <  1:386788/148‑1 (MQ=255)
   gcagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTc                                                                                                                                            <  1:819696/148‑1 (MQ=255)
   gcagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGAGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTc                                                                                                                                            <  1:2697719/148‑1 (MQ=255)
         gAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGa                                                                                                                                      >  2:1544787/1‑148 (MQ=255)
          aataatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGaa                                                                                                                                     >  1:2055462/1‑148 (MQ=255)
           ataatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAAt                                                                                                                                    >  2:2784821/1‑148 (MQ=255)
               tGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCat                                                                                                                                >  2:3531752/1‑148 (MQ=255)
               tGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCat                                                                                                                                >  2:1599465/1‑148 (MQ=255)
                ggTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCata                                                                                                                               >  1:1713833/1‑148 (MQ=255)
                 gTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCatat                                                                                                                              >  1:2493477/1‑148 (MQ=255)
                  ttttAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATAtt                                                                                                                             >  2:11679/1‑148 (MQ=255)
                   tttAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTg                                                                                                                            >  1:2901499/1‑148 (MQ=255)
                        tCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGAc                                                                                                                       <  1:3242020/148‑1 (MQ=255)
                          aGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCa                                                                                                                     <  1:3798880/148‑1 (MQ=255)
                           gCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAg                                                                                                                    <  2:303976/148‑1 (MQ=255)
                           gCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAg                                                                                                                    >  1:486324/1‑148 (MQ=255)
                           gCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAg                                                                                                                    >  1:484358/1‑148 (MQ=255)
                             ggACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTg                                                                                                                  <  1:3810342/148‑1 (MQ=255)
                              gACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGc                                                                                                                 >  1:3552104/1‑148 (MQ=255)
                                cgcgGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCgg                                                                                                               <  2:2402189/148‑1 (MQ=255)
                                  cgGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGtt                                                                                                             <  2:3404383/148‑1 (MQ=255)
                                    gggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAt                                                                                                           >  2:1485542/1‑148 (MQ=255)
                                    gggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAt                                                                                                           >  2:1637350/1‑148 (MQ=255)
                                    gggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATACATATTGCTGACCAGTGCGGTTAt                                                                                                           >  2:2301701/1‑148 (MQ=255)
                                     tgttaaGTACTGGCCTGCACAAGTCTTGGTTTCAGTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCTCGACGTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAtt                                                                                                          <  2:3263788/144‑1 (MQ=255)
                                     ggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAtt                                                                                                          <  1:1270611/148‑1 (MQ=255)
                                       tGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTg                                                                                                        <  2:1360859/148‑1 (MQ=255)
                                         aaGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACATCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGtt                                                                                                      <  2:3670023/148‑1 (MQ=255)
                                         aaGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGtt                                                                                                      >  1:3522147/1‑148 (MQ=255)
                                         aaGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGtt                                                                                                      >  1:672617/1‑148 (MQ=255)
                                             aCTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCgg                                                                                                  <  1:1313925/148‑1 (MQ=255)
                                             aCTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCgg                                                                                                  <  2:3766462/148‑1 (MQ=255)
                                              cTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACCTATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGc                                                                                                 <  2:2125606/148‑1 (MQ=255)
                                              cTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGc                                                                                                 <  2:2125931/148‑1 (MQ=255)
                                               tGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCt                                                                                                <  2:3519257/148‑1 (MQ=255)
                                                    tGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTg                                                                                           >  1:934709/1‑148 (MQ=255)
                                                    tGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTg                                                                                           <  1:1336014/148‑1 (MQ=255)
                                                     gCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGa                                                                                          <  2:1136015/148‑1 (MQ=255)
                                                     gCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGa                                                                                          <  2:1135868/148‑1 (MQ=255)
                                                        ccAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtgg                                                                                       >  2:1279066/1‑148 (MQ=255)
                                                        ccAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtgg                                                                                       >  2:1278739/1‑148 (MQ=255)
                                                        ccAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtgg                                                                                       <  1:2159852/148‑1 (MQ=255)
                                                         cAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtggt                                                                                      >  1:3021378/1‑148 (MQ=255)
                                                           gACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtggtgg                                                                                    >  1:2734915/1‑148 (MQ=255)
                                                             cTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTa                                                                                  <  1:241114/148‑1 (MQ=255)
                                                             cTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTa                                                                                  <  1:1009400/148‑1 (MQ=255)
                                                                ggTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTa                                                                               <  2:412616/148‑1 (MQ=255)
                                                                 gTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTaa                                                                              >  1:2183471/1‑148 (MQ=255)
                                                                 gTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTaa                                                                              >  1:2183272/1‑148 (MQ=255)
                                                                      aTTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTa                                                                         >  2:495491/1‑148 (MQ=255)
                                                                       ttCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTAt                                                                        >  2:3588967/1‑148 (MQ=255)
                                                                        tCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATg                                                                       >  1:376923/1‑148 (MQ=255)
                                                                        tCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATg                                                                       >  1:3489900/1‑148 (MQ=255)
                                                                        tCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATg                                                                       >  1:2589431/1‑148 (MQ=255)
                                                                           cgTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaa                                                                    >  2:61179/1‑148 (MQ=255)
                                                                           cgTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaa                                                                    >  2:57622/1‑148 (MQ=255)
                                                                           cgTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaa                                                                    >  1:2756475/1‑148 (MQ=255)
                                                                            gTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaaa                                                                   <  1:2315104/148‑1 (MQ=255)
                                                                             ttcttcTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaaaa                                                                  <  2:2331756/148‑1 (MQ=255)
                                                                              tcttctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAAc                                                                 >  1:1563394/1‑148 (MQ=255)
                                                                              tcttctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAAc                                                                 >  1:3689430/1‑148 (MQ=255)
                                                                                ttctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTg                                                               <  2:3676297/148‑1 (MQ=255)
                                                                                      aaCAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTa                                                         <  2:1050252/148‑1 (MQ=255)
                                                                                        cAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTAcc                                                       >  2:3571500/1‑148 (MQ=255)
                                                                                        cAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTAGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTAcc                                                       >  1:656490/1‑148 (MQ=255)
                                                                                            ccTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCt                                                   <  2:2557873/148‑1 (MQ=255)
                                                                                             cTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtt                                                  >  2:2254748/1‑148 (MQ=255)
                                                                                             cTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtt                                                  <  1:647042/148‑1 (MQ=255)
                                                                                             cTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtt                                                  >  2:683103/1‑148 (MQ=255)
                                                                                              ttCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCttt                                                 <  2:3795487/148‑1 (MQ=255)
                                                                                              ttCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCttt                                                 <  2:2925697/148‑1 (MQ=255)
                                                                                               tCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtttt                                                <  2:3552104/148‑1 (MQ=255)
                                                                                               tCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtttt                                                >  1:343785/1‑148 (MQ=255)
                                                                                               tCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtttt                                                >  1:1013058/1‑148 (MQ=255)
                                                                                                   aCATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCAGATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTc                                            >  1:1637571/1‑148 (MQ=255)
                                                                                                     atatCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCAt                                          >  1:1242865/1‑148 (MQ=255)
                                                                                                         cGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTc                                      >  2:69594/1‑148 (MQ=255)
                                                                                                         cGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTc                                      >  2:66294/1‑148 (MQ=255)
                                                                                                         cGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTc                                      >  1:3500132/1‑148 (MQ=255)
                                                                                                          gCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCg                                     <  1:714715/148‑1 (MQ=255)
                                                                                                           ccGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGc                                    <  2:2176841/148‑1 (MQ=255)
                                                                                                            cGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGcc                                   <  1:3236460/148‑1 (MQ=255)
                                                                                                               aGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTaa                                >  1:2247316/1‑148 (MQ=255)
                                                                                                               aGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTaa                                >  1:464164/1‑148 (MQ=255)
                                                                                                                          gACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCa                     <  2:3115820/148‑1 (MQ=255)
                                                                                                                          gACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCa                     <  2:1996166/148‑1 (MQ=255)
                                                                                                                           aCTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAt                    >  1:2579787/1‑148 (MQ=255)
                                                                                                                            cTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAtt                   <  1:1599465/148‑1 (MQ=255)
                                                                                                                             ttACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAttt                  <  1:1072148/148‑1 (MQ=255)
                                                                                                                              tACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAtttt                 >  1:3377434/1‑148 (MQ=255)
                                                                                                                               aCCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAttttt                <  1:2903825/148‑1 (MQ=255)
                                                                                                                                  gTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCgg             <  2:2151919/148‑1 (MQ=255)
                                                                                                                                   tGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGt            <  1:1279066/148‑1 (MQ=255)
                                                                                                                                       cgcCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaa        <  2:2329947/148‑1 (MQ=255)
                                                                                                                                       cgcCGTTGTTCCAGTCCAGCGAATCCATAGTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaa        <  1:3103910/148‑1 (MQ=255)
                                                                                                                                          cGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactg     <  1:3531743/148‑1 (MQ=255)
                                                                                                                                            ttGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgca   <  2:3001534/148‑1 (MQ=255)
                                                                                                                                             tGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgcaa  >  2:3196899/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGAACTGCAA  >  CP009273/4572376‑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: ‑