Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,842,438 A→C 5.4% intergenic (+24/+135) ascB → / ← hycI cryptic 6‑phospho‑beta‑glucosidase/protease involved in processing C‑terminal end of HycE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,842,4380AC 5.4% 337.6 / 6.1 122intergenic (+24/+135)ascB/hycIcryptic 6‑phospho‑beta‑glucosidase/protease involved in processing C‑terminal end of HycE
Reads supporting (aligned to +/- strand):  ref base A (53/62);  new base C (4/3);  total (57/65)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

Only 100 of 122 total aligned reads displayed.

CGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTC  >  NC_000913/2842292‑2842584
                                                                                                                                                  |                                                                                                                                                  
cccTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACTGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAAc                                                                                                                                                  <  1:1025031/147‑1 (MQ=255)
  cTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGc                                                                                                                                                <  1:183073/149‑1 (MQ=255)
  cTACGGCTTTGGCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAACAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTCACAGTGCCGGATGCGGCGTGAACGc                                                                                                                                                <  1:295061/149‑1 (MQ=255)
   ttcGGCTTTGTCTTTGTTGACCGTGACGACTCAGGCAACGGTACGAGGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAACTGATTGCGAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGcc                                                                                                                                               <  2:1321409/147‑1 (MQ=255)
        ctttgtctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                          <  2:522287/149‑1 (MQ=255)
          ttgtctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                        <  2:448795/149‑1 (MQ=255)
            gtctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                      <  1:694585/149‑1 (MQ=255)
            gtctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                      <  2:408788/149‑1 (MQ=255)
            ggctttgtTGACCGGGACGACGCTGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGGGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                      <  2:628434/147‑1 (MQ=255)
             tctttgtTGACCGTGACGACGCATGCAAAGGTACGCTGAAGCTCACGAGTAAGAAATCATTCTGGGGGTATAAAAAAGTTATTTCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                     <  1:1439863/149‑1 (MQ=38)
               tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                   >  2:806221/1‑149 (MQ=255)
                ttgttgACCGTGACGACGCCGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATAAGGCGTTCACGCCGCATCCGGcactg                                                                                                                                  <  2:270850/149‑5 (MQ=255)
                ttgttgACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                  >  1:1236393/1‑149 (MQ=255)
                tagttgACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCTTTCTGGTGGTAAAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGCTAACAGTGCCTGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                  <  1:1018011/147‑1 (MQ=255)
                 tgtcgAACGTGACGCCGCACGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGGAATGGGGAAGATTTCGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa                                                                                                                                 <  2:1067163/149‑1 (MQ=39)
                  gttgACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAt                                                                                                                                >  2:336690/1‑149 (MQ=255)
                     gACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTc                                                                                                                             <  2:1525725/149‑1 (MQ=255)
                      aCCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCg                                                                                                                            >  1:50918/1‑149 (MQ=255)
                        cGTGCCGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCAGGTGGTATAAAAAAGGGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCa                                                                                                                          <  1:1018651/149‑1 (MQ=255)
                        cGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCa                                                                                                                          <  1:791867/149‑1 (MQ=255)
                           gacgacGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGa                                                                                                                       >  1:1113562/1‑149 (MQ=255)
                            actaagaaGGCAAGGGAACTCTTACGCGCACGCTTAAGAAATCATTCTGGTGGTATAAAAAAGTTATTGCCAGTAATGGGGAAGATTTAGAGTATGTATCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGaa                                                                                                                      <  1:1279497/142‑1 (MQ=37)
                              gccgCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAACGTGATTGCCAGTAATGGGGAAGATTTAGAGGAGGTAACAGTGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACATGTCGCCCGAAcc                                                                                                                    <  1:377960/147‑1 (MQ=255)
                              gacgCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAAcc                                                                                                                    >  1:529712/1‑149 (MQ=255)
                                cgCAGGCAACGGTACGCTTACGCGCACTCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCAt                                                                                                                  <  2:132761/149‑1 (MQ=255)
                                 gCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATg                                                                                                                 >  1:642428/1‑149 (MQ=255)
                                     gCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAgg                                                                                                             >  2:370926/1‑149 (MQ=255)
                                      cAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGc                                                                                                            >  1:1181815/1‑149 (MQ=255)
                                       aaCGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGcc                                                                                                           <  2:1383051/149‑1 (MQ=255)
                                          ggTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGTGGAATATTTAGAGTAGGTAACAGTGCAGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGa                                                                                                        <  2:225214/149‑1 (MQ=255)
                                          ggTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGa                                                                                                        <  1:464908/149‑1 (MQ=255)
                                            tACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATa                                                                                                      >  2:1432079/1‑149 (MQ=255)
                                             aCGCGGACGCGCACGCGTAAGACATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTACTGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATaa                                                                                                     <  1:702084/149‑1 (MQ=255)
                                               gCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagc                                                                                                   >  1:816638/1‑149 (MQ=255)
                                               gCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagc                                                                                                   >  2:750553/1‑149 (MQ=255)
                                                cTGACGCGCACGAGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagcg                                                                                                  <  1:1398527/149‑1 (MQ=255)
                                                 tGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagcgc                                                                                                 <  1:574146/149‑1 (MQ=255)
                                                       gcACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAgcgcagcgca                                                                                           >  2:1293217/1‑149 (MQ=255)
                                                           gcgTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCgg                                                                                       >  2:32590/1‑149 (MQ=255)
                                                            cgTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCggg                                                                                      <  1:1318202/149‑1 (MQ=255)
                                                              tAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGATCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAGTCGCACGAACCATAAAAGCCCGATAACCGCCACGCAACgggaa                                                                                    >  2:1011892/1‑147 (MQ=38)
                                                                 taactCATTCTGGCGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAAtt                                                                                 <  1:1547756/145‑1 (MQ=255)
                                                                 gAAATCATTCTGGTGGTATAAAAAAGTGATTGTCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAAtt                                                                                 >  2:1252070/1‑149 (MQ=255)
                                                                     tCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGc                                                                             >  2:489739/1‑149 (MQ=255)
                                                                        ttCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAGCCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCAtt                                                                          <  1:854213/149‑1 (MQ=255)
                                                                            ggtggtATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGc                                                                      <  1:701670/149‑1 (MQ=255)
                                                                            ggtggtATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGc                                                                      <  2:900893/149‑1 (MQ=255)
                                                                            ggtggaCTAAAAAAGTGATTGCCAGTACTGGGGAAGATTTAGAGTAGGTACCAGGGCCGGATGCGGCGTGCACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGc                                                                      <  1:1133133/149‑1 (MQ=255)
                                                                             gtggtATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCa                                                                     >  2:821601/1‑149 (MQ=255)
                                                                             gtggtATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCa                                                                     >  1:263531/1‑149 (MQ=255)
                                                                                gtATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCa                                                                  <  2:642428/149‑1 (MQ=255)
                                                                                     aaaaaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCt                                                             >  1:1334605/1‑149 (MQ=255)
                                                                                             aTTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCATCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCAcc                                                     <  1:1204200/149‑1 (MQ=255)
                                                                                             aTTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCAcc                                                     >  1:1541474/1‑149 (MQ=255)
                                                                                                tccAGTAATGGTGTAGATTTAGAGTAGGTAACAGTGCCGGATCCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATATGGCAATTTAGCATTTGGCATCATAGCTGGTTCACCCCa                                                  <  1:1052859/148‑1 (MQ=38)
                                                                                                 ccAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCag                                                 >  2:752535/1‑148 (MQ=255)
                                                                                                     tAATGGGGAAGTTTTAGAGTAGGTAACAGTTCGGGATTCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGAGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAg                                             <  2:944141/149‑1 (MQ=39)
                                                                                                     tAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAg                                             >  2:1119979/1‑149 (MQ=255)
                                                                                                     tAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAg                                             <  2:728098/149‑1 (MQ=255)
                                                                                                       aTGTGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGaa                                           <  1:1144490/149‑1 (MQ=255)
                                                                                                        tGGGGAATATTTAGAGAAGGTAACAGTGCCGGTTGCGGCGGGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGCAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAAc                                          <  1:853110/149‑1 (MQ=39)
                                                                                                          gggAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGACCCATGAAGGCCCGATAAGCGCAGCGCAGCGGGCATTTTAGCAGTTGGCATCAGAGCTGGTTCACCCCCTCAAGAACAt                                        >  2:1132259/1‑149 (MQ=38)
                                                                                                          gggAAGATTTAGAGTACGTAACAGTGCCTGATGCGGCGTCAACGCCTTATACTGCCTACATGTCGCACGAAGCATGAAGGCCCTATAAGAGCAGCGCATCGGGCAAATTAGCATATGGCATCAGAGCTGGTTCACCCCAACAAGAACAt                                        >  2:382147/1‑149 (MQ=255)
                                                                                                            gAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATcc                                      <  2:214535/149‑1 (MQ=255)
                                                                                                               gATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTg                                   >  2:169117/1‑149 (MQ=255)
                                                                                                                 tttAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAAtt                                                                                 <  2:1492081/101‑1 (MQ=255)
                                                                                                                 tttAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAAtt                                                                                 >  1:1492081/1‑101 (MQ=255)
                                                                                                                 tttAGAGTAGGTAACAGTGCCGGATGAGGCGTGAACGCCTTATCCGGCCTACTTGTCGCACGAACCATGAATGCCCGATAAGCGCAGCGTATCGGGCAATTTAGCAGTTGGCATCAGATCAGGTTCACTCCATCAAGAACAACCCAGTc                                 >  2:855133/1‑149 (MQ=38)
                                                                                                                  ttAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCACCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTcc                                >  2:1488170/1‑149 (MQ=255)
                                                                                                                  ttAGAGTAGGTAACAGAGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGACCCATCAAGGCCCGATAAGCGCATCGCATCGGGCAATTTAGCATTTGGAATCAGAGCTGGTTCACACCATCAAGAACATCCCTGtgc                                >  2:327146/1‑147 (MQ=38)
                                                                                                                    agagTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTg                              >  2:776524/1‑149 (MQ=255)
                                                                                                                     gagTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTgt                             >  2:470053/1‑148 (MQ=255)
                                                                                                                       gTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGAtt                           <  2:1544665/149‑1 (MQ=255)
                                                                                                                        tAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCGTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCAGCAAGCACATCCCTGTCCTGATTc                          >  1:394647/1‑149 (MQ=255)
                                                                                                                           gTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCtt                       >  2:154621/1‑149 (MQ=255)
                                                                                                                             aaCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTaa                     <  1:598826/149‑1 (MQ=255)
                                                                                                                             aaCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTaa                     >  1:664744/1‑149 (MQ=255)
                                                                                                                              aCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAAt                    <  1:750553/149‑1 (MQ=255)
                                                                                                                              aCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAAt                    <  1:450342/149‑1 (MQ=255)
                                                                                                                              aCAGGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAAt                    <  2:433141/149‑1 (MQ=255)
                                                                                                                               cAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATg                   <  2:500951/149‑1 (MQ=255)
                                                                                                                                taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACCTGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAGCATCCCTGTGCTGATTCCTTAATGa                  >  1:1130918/4‑149 (MQ=255)
                                                                                                                                taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACATGTCGCTCGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGa                  >  1:994267/4‑149 (MQ=255)
                                                                                                                                aGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGa                  <  1:389432/149‑1 (MQ=255)
                                                                                                                                aGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACAGCCCTGTCCTGATTCCTTAATGa                  >  1:366764/1‑149 (MQ=255)
                                                                                                                                 gTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGaa                 <  1:776524/149‑1 (MQ=255)
                                                                                                                                  tGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGGCCTGATTCCTTAATGaaa                >  2:1568546/1‑149 (MQ=255)
                                                                                                                                  tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGaaa                <  1:763348/149‑1 (MQ=255)
                                                                                                                                  tGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGACGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGaaa                <  1:1553486/149‑1 (MQ=255)
                                                                                                                                    ccGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGCTTCCTTAATGaaaaa              >  2:260543/1‑149 (MQ=255)
                                                                                                                                    ccGGATAAGGCGTTCACGCCGCATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGaaaaa              >  2:408152/1‑149 (MQ=255)
                                                                                                                                     cGGATTCTGCGTGAACGCCTTTTCCGGCCTACATGTGGCACTAACCTTGAATGCCCTATAAGCGGAGCGCATGGGGCAATTTAGCATTTGGCATCATAGATGGTTCACCCCATCAATAACATCCCTGTCCTGATTCCTTAATGAAAAac             <  2:932318/149‑1 (MQ=255)
                                                                                                                                     cGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGAACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAac             <  1:1011892/149‑1 (MQ=255)
                                                                                                                                       gataaggcgTTCACGCCGCATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAAGTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGAGTCCTTAATGAAAAacta           >  1:1121000/6‑149 (MQ=255)
                                                                                                                                       gATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAacta           <  1:781295/149‑1 (MQ=255)
                                                                                                                                          gcggcgTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACACCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAactacta        >  1:820319/1‑148 (MQ=255)
                                                                                                                                          ccggcgTGAACGCCTTATACGCCCTACATGACGCACGAACCATGAAGGCACGATAAGCGCAGCGCAACGCGCAATATAGCATTAGGCATCAGAGCTGGTTCACCCCAACAAGAACATCCCTGTCCAGATTCCTTAATGAAAAactacta        >  2:50171/2‑148 (MQ=255)
                                                                                                                                            ggcgTGAACGCCTTATCCGGCCTACAAGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAAGTTAGCATTTGGCATCAGCGCTGGTTAACCCCATCAAGAACATCCCTGTCCTGAGTCCTTAATGAAAAACTACgctt      >  1:300571/1‑149 (MQ=39)
                                                                                                                                               gTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTcttctt   >  2:831021/1‑149 (MQ=255)
                                                                                                                                                tGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCttcttc  <  2:1181815/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTC  >  NC_000913/2842292‑2842584

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: