Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,572,520 A→T 92.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,5200AT92.8% 566.7 / 27.1 181Y116N (TAC→AAC) hsdMDNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base A (5/8);  new base T (73/95);  total (78/103)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 183 total aligned reads displayed.

TTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGAACTGCAA  >  CP009273/4572375‑4572664
                                                                                                                                                 |                                                                                                                                                
ttCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTcc                                                                                                                                                >  1:2808557/1‑148 (MQ=255)
   agcagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGt                                                                                                                                             <  1:466447/148‑1 (MQ=255)
   agcagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGt                                                                                                                                             <  2:1179419/148‑1 (MQ=255)
    gcagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTc                                                                                                                                            <  1:3279740/148‑1 (MQ=255)
     cagATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTcc                                                                                                                                           <  2:1025187/148‑1 (MQ=255)
          gAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGa                                                                                                                                      >  1:3231108/1‑148 (MQ=255)
           aataatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGaa                                                                                                                                     >  1:2700650/1‑148 (MQ=255)
             taatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATc                                                                                                                                   <  1:102419/148‑1 (MQ=255)
             taatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATc                                                                                                                                   <  1:102316/148‑1 (MQ=255)
              aatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCTAGCGAATcc                                                                                                                                  >  2:2300068/1‑148 (MQ=255)
              aatGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATcc                                                                                                                                  <  2:1881339/148‑1 (MQ=255)
               atGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCa                                                                                                                                 <  2:1388486/148‑1 (MQ=255)
                tGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCag                                                                                                                                >  2:2358692/1‑147 (MQ=255)
                  gTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCatat                                                                                                                              <  2:3259543/148‑1 (MQ=255)
                   ttttAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATAtt                                                                                                                             <  2:2473494/148‑1 (MQ=255)
                    tttAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTg                                                                                                                            >  1:3224772/1‑148 (MQ=255)
                    tttAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTg                                                                                                                            >  2:457961/1‑148 (MQ=255)
                     ttAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGc                                                                                                                           <  1:3050905/148‑1 (MQ=255)
                     ttAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGc                                                                                                                           >  2:386128/1‑148 (MQ=255)
                        aTCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGa                                                                                                                        >  2:762665/1‑148 (MQ=255)
                          cAGCGGAAGAGGGGAGAAGTACTGGCCTGCACAAGACTTGGTTTCATTCTCGTTCTTCTGAAACAGCCCATCGTACATATCGCAGAAGTCATCGCTCGACTTACCGGGCGCGCCGTTGTTCAAGTCCAGCGAATCCATATTGCTGAca                                                                                                                      >  2:1157415/1‑147 (MQ=255)
                                aCGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCg                                                                                                                <  1:1112566/148‑1 (MQ=255)
                                aCGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCg                                                                                                                >  1:389047/1‑148 (MQ=255)
                                aCGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCg                                                                                                                >  2:1537209/1‑148 (MQ=255)
                                 cgcgGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCgg                                                                                                               <  2:2385938/148‑1 (MQ=255)
                                      ggTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTAtt                                                                                                          <  2:1794175/148‑1 (MQ=255)
                                         gAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGt                                                                                                       <  1:654858/148‑1 (MQ=255)
                                         gAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGt                                                                                                       >  1:2061903/1‑148 (MQ=255)
                                             tACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCg                                                                                                   >  2:2522546/1‑148 (MQ=255)
                                               cTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGc                                                                                                 <  2:1934763/148‑1 (MQ=255)
                                               cTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGAGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGc                                                                                                 <  1:3361172/148‑1 (MQ=255)
                                                tGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCt                                                                                                >  2:1693232/1‑148 (MQ=255)
                                                  gCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCg                                                                                              >  1:2937035/1‑148 (MQ=255)
                                                          cAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtggt                                                                                      >  2:379827/1‑148 (MQ=255)
                                                          cAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtggt                                                                                      >  2:105688/1‑148 (MQ=255)
                                                           aGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtggtg                                                                                     >  2:346422/1‑148 (MQ=255)
                                                            gACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGAtggtgg                                                                                    >  2:302116/1‑148 (MQ=255)
                                                             aCTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATggtggt                                                                                   >  2:3046577/1‑148 (MQ=255)
                                                                tGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACt                                                                                >  2:2711583/1‑148 (MQ=255)
                                                                 ggTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTa                                                                               <  1:2255247/148‑1 (MQ=255)
                                                                   tttCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAAc                                                                             <  2:2237329/148‑1 (MQ=255)
                                                                   tttCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAAc                                                                             <  1:2190147/148‑1 (MQ=255)
                                                                    ttCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACa                                                                            <  1:1432109/148‑1 (MQ=255)
                                                                    ttCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACa                                                                            <  1:2022366/148‑1 (MQ=255)
                                                                      cATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACAtt                                                                          <  1:1569087/148‑1 (MQ=255)
                                                                           gcgTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaa                                                                     >  1:3534430/1‑148 (MQ=255)
                                                                            cgTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaa                                                                    <  2:2808557/148‑1 (MQ=255)
                                                                            cgTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGGTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaa                                                                    <  2:2618046/148‑1 (MQ=255)
                                                                              ttcttcTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGaaaaa                                                                  <  1:2522546/148‑1 (MQ=255)
                                                                               tcttctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAAc                                                                 >  2:1840298/1‑148 (MQ=255)
                                                                                 ttctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTg                                                               <  1:3513145/148‑1 (MQ=255)
                                                                                  tctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGc                                                              <  1:2472412/148‑1 (MQ=255)
                                                                                  tctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGc                                                              >  2:2027443/1‑148 (MQ=255)
                                                                                  tctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGc                                                              <  1:2949319/148‑1 (MQ=255)
                                                                                  tctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGc                                                              >  2:3366231/1‑148 (MQ=255)
                                                                                   ctGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGcc                                                             >  2:1120236/1‑148 (MQ=255)
                                                                                     gCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTg                                                           <  1:3037120/148‑1 (MQ=255)
                                                                                      cAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGt                                                          <  2:2863148/148‑1 (MQ=255)
                                                                                       aaCAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTa                                                         <  2:590139/148‑1 (MQ=255)
                                                                                       aaCAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTa                                                         <  2:588799/148‑1 (MQ=255)
                                                                                        aCAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTAc                                                        >  1:1462303/1‑148 (MQ=255)
                                                                                         catcccTTCGTACATATCGCAGAAGTCATCTCGCGACTTACAGTGCGCGCAGTTGTACAATTCCAGCGAATCCATATTGCTGACAAGTGAGGTTATTTGTTTCGGCTCGGTGATTGTGGTACTAACATTAAGAAAAACTGACTTTAcc                                                       >  2:652046/4‑148 (MQ=255)
                                                                                           gCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAg                                                     <  1:998789/148‑1 (MQ=255)
                                                                                           gCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAg                                                     >  2:1999086/1‑148 (MQ=255)
                                                                                           gCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAg                                                     >  2:1999070/1‑148 (MQ=255)
                                                                                            cccTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGc                                                    <  2:2990828/148‑1 (MQ=255)
                                                                                            cccTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGc                                                    <  2:589646/148‑1 (MQ=255)
                                                                                            cccTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAAATGCCTGTACCAGc                                                    >  1:1885703/1‑148 (MQ=255)
                                                                                              cTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtt                                                  <  2:3032769/148‑1 (MQ=255)
                                                                                               ttCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCttt                                                 >  1:68727/1‑148 (MQ=255)
                                                                                                tCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCtttt                                                >  1:932152/1‑148 (MQ=255)
                                                                                                 cGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCttttt                                               <  2:1149771/148‑1 (MQ=255)
                                                                                                    aCATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTc                                            <  2:1642101/148‑1 (MQ=255)
                                                                                                      atatCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCAt                                          <  1:3031745/148‑1 (MQ=255)
                                                                                                      atatCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCAt                                          <  1:1184751/148‑1 (MQ=255)
                                                                                                          cGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTc                                      >  2:213659/1‑148 (MQ=255)
                                                                                                            ccGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGc                                    <  2:1247287/148‑1 (MQ=255)
                                                                                                            ccGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGc                                    <  2:3023519/148‑1 (MQ=255)
                                                                                                             cGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGcc                                   <  2:68727/148‑1 (MQ=255)
                                                                                                                 gTCATCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCTAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTaaa                               <  1:2300068/148‑1 (MQ=255)
                                                                                                                    aTCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGc                            <  1:2598218/148‑1 (MQ=255)
                                                                                                                    aTCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGc                            <  1:1528096/148‑1 (MQ=255)
                                                                                                                     tCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCa                           <  2:1544940/148‑1 (MQ=255)
                                                                                                                     tCGCGCGACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCa                           >  2:3531262/1‑148 (MQ=255)
                                                                                                                           gACTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCa                     <  2:1964602/148‑1 (MQ=255)
                                                                                                                            aCTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAt                    >  1:1739289/1‑148 (MQ=255)
                                                                                                                            aCTTACCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAt                    >  1:847047/1‑148 (MQ=255)
                                                                                                                             cTTACCGTGCGCGCCGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAtt                   <  1:386128/148‑1 (MQ=255)
                                                                                                                                aCCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCGTCTTAGCCTAAATGCACGAGCAttttt                >  1:953704/1‑148 (MQ=255)
                                                                                                                                aCCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAttttt                <  1:174527/148‑1 (MQ=255)
                                                                                                                                aCCGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCAttttt                <  1:1330278/148‑1 (MQ=255)
                                                                                                                                 ccGTGCGCGCCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTc               <  1:1537209/148‑1 (MQ=255)
                                                                                                                                       gcgcCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGa         <  1:3172382/148‑1 (MQ=255)
                                                                                                                                        cgcCGTTGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaa        <  1:481918/148‑1 (MQ=255)
                                                                                                                                           cGTTGTACCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactg     >  1:2716350/1‑148 (MQ=255)
                                                                                                                                             ttGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgca   >  2:839309/1‑148 (MQ=255)
                                                                                                                                             ttGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgca   >  1:168552/1‑148 (MQ=255)
                                                                                                                                             ttGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgca   >  1:2142023/1‑148 (MQ=255)
                                                                                                                                              tGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgcaa  >  2:1446230/1‑148 (MQ=255)
                                                                                                                                              tGTTCCAGTCCAGCGAATCCATATTGCTGACCAGTGCGGTTATTTGTTTCGGCTCGGTGATGGTGGTACTAACATTATGAAAAACTGCCTGTACCAGCTTTTTGTCATCTTCGCCTAAATGCACGAGCATTTTTCGGTAGaactgcaa  <  1:3114039/148‑1 (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: ‑