Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,842,438 A→C 5.3% 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.3% 394.8 / 10.8 113intergenic (+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 (56/51);  new base C (4/2);  total (60/53)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 115 total aligned reads displayed.

AACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCA  >  NC_000913/2842290‑2842586
                                                                                                                                                    |                                                                                                                                                    
aaCGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGa                                                                                                                                                      >  1:370292/1‑149 (MQ=255)
  cGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAAc                                                                                                                                                    <  3:163386/149‑1 (MQ=255)
    cTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGc                                                                                                                                                  >  1:732575/1‑149 (MQ=255)
     tACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACTGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGcc                                                                                                                                                 >  4:1001178/1‑149 (MQ=255)
     tACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGcc                                                                                                                                                 >  4:797379/1‑149 (MQ=255)
       cGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATAAGGCGTTCACGccgc                                                                                                                                               <  1:1089513/149‑3 (MQ=255)
         gCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                             >  1:16085/1‑149 (MQ=255)
            ttgtctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                          <  2:729978/149‑1 (MQ=255)
            ttgtctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                          <  4:596529/149‑1 (MQ=255)
               tctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                       >  3:162205/1‑149 (MQ=255)
               tctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                       <  4:307568/149‑1 (MQ=255)
               tctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                       >  4:1059259/1‑149 (MQ=255)
                 tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                          <  1:666207/144‑1 (MQ=255)
                 tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                          >  2:666207/1‑144 (MQ=255)
                 tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                     >  1:216702/1‑149 (MQ=255)
                 tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                     >  1:280550/1‑149 (MQ=255)
                 tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                     >  3:49761/1‑149 (MQ=255)
                 tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                         <  2:145416/145‑1 (MQ=255)
                 tttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                         >  1:145416/1‑145 (MQ=255)
                     ttgACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATg                                                                                                                                 >  4:735450/1‑149 (MQ=255)
                     ttgACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATg                                                                                                                                 >  1:287259/1‑149 (MQ=255)
                                gacgCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGTTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAAcc                                                                                                                      <  2:909592/149‑1 (MQ=255)
                                gacgCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                       >  4:404035/1‑132 (MQ=255)
                                gacgCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                       <  3:404035/132‑1 (MQ=255)
                                  cgCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCAt                                                                                                                    >  3:452295/1‑149 (MQ=255)
                                  cgCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCAt                                                                                                                    <  4:902764/149‑1 (MQ=255)
                                  cgCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCAt                                                                                                                    <  4:890936/149‑1 (MQ=255)
                                   gCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATg                                                                                                                   >  2:547305/1‑149 (MQ=255)
                                    cAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGa                                                                                                                  <  3:797379/149‑1 (MQ=255)
                                      ggCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAg                                                                                                                >  2:111415/1‑149 (MQ=255)
                                      ggCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAg                                                                                                                <  2:501015/149‑1 (MQ=255)
                                       gCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAgg                                                                                                               >  3:562668/1‑149 (MQ=255)
                                         aaCGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGcc                                                                                                             >  4:953786/1‑149 (MQ=255)
                                          aCGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGccc                                                                                                            <  4:630945/149‑1 (MQ=255)
                                            ggTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGa                                                                                                          >  4:589983/1‑149 (MQ=255)
                                             gTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGAt                                                                                                         <  1:623054/149‑1 (MQ=255)
                                               aCGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATaa                                                                                                       <  3:1049073/149‑1 (MQ=255)
                                                 gCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagc                                                                                                     >  2:845257/1‑149 (MQ=255)
                                                     aCGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAagcgcag                                                                                                 <  3:589983/149‑1 (MQ=255)
                                                           aCGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATc                                                                                           >  4:468883/1‑149 (MQ=255)
                                                            cgcgTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCg                                                                                          <  2:314963/149‑1 (MQ=255)
                                                             gcgTAAGAAATCATTCTGGTTGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCgg                                                                                         <  1:599033/149‑1 (MQ=255)
                                                               gTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGc                                                                                       >  4:407827/1‑149 (MQ=255)
                                                                tAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCa                                                                                      >  1:827138/1‑149 (MQ=255)
                                                                   gAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGTAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAAtt                                                                                   >  3:944435/1‑149 (MQ=255)
                                                                   gAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGGAGCGCATCGGGCAAtt                                                                                   >  3:922815/1‑149 (MQ=255)
                                                                   gAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGACCGATAAGCGCAGCGCATCGGACAAtt                                                                                   >  4:935028/1‑149 (MQ=255)
                                                                     aaTCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTa                                                                                 >  4:102009/1‑149 (MQ=255)
                                                                     aaTCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTa                                                                                 <  3:450558/149‑1 (MQ=255)
                                                                      aTCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAg                                                                                >  3:869572/1‑149 (MQ=255)
                                                                       tCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGc                                                                               >  2:243100/1‑149 (MQ=255)
                                                                         aTTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCAt                                                                             >  2:325491/1‑149 (MQ=255)
                                                                          ttCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCa                                                                              <  1:165618/147‑1 (MQ=255)
                                                                          ttCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCa                                                                              >  2:165618/1‑147 (MQ=255)
                                                                            cTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTg                                                                          <  2:287259/149‑1 (MQ=255)
                                                                               gtggtATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCa                                                                       >  4:778147/1‑149 (MQ=255)
                                                                                   tataAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCag                                                                   >  4:833289/1‑149 (MQ=255)
                                                                                   tataAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCag                                                                   >  4:816080/1‑149 (MQ=255)
                                                                                    ataAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCaga                                                                  <  4:601764/149‑1 (MQ=255)
                                                                                     taAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATaa                                                                                                       <  3:99971/111‑1 (MQ=255)
                                                                                      aaaaaaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGATCCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGc                                                                >  1:488712/1‑149 (MQ=255)
                                                                                      aaaaaaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGc                                                                <  3:740703/149‑1 (MQ=255)
                                                                                       aaaaaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCt                                                               >  4:743276/1‑149 (MQ=255)
                                                                                        aaaaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTg                                                              >  1:781788/1‑149 (MQ=255)
                                                                                        aaaaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTg                                                              <  3:724646/149‑1 (MQ=255)
                                                                                        aaaaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTg                                                              >  4:639250/1‑149 (MQ=255)
                                                                                        aaaaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTg                                                              >  4:454113/1‑149 (MQ=255)
                                                                                        aaaaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTg                                                              >  4:517244/1‑149 (MQ=255)
                                                                                          aaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGt                                                            >  4:26679/1‑149 (MQ=255)
                                                                                          aaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGt                                                            <  2:216702/149‑1 (MQ=255)
                                                                                            gTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTc                                                          <  4:562668/149‑1 (MQ=255)
                                                                                             tGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCa                                                         <  1:369628/149‑1 (MQ=255)
                                                                                               aTTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCAcc                                                       <  3:517244/149‑1 (MQ=255)
                                                                                                ttGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCAccc                                                      >  3:799069/1‑149 (MQ=255)
                                                                                                  gCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCa                                                    >  3:422647/1‑149 (MQ=255)
                                                                                                      gTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCaa                                                <  2:370292/149‑1 (MQ=255)
                                                                                                            gggTAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCCAGAACAt                                          >  4:896284/1‑149 (MQ=255)
                                                                                                             ggAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAATGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATc                                         >  3:704948/1‑149 (MQ=255)
                                                                                                              gAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATcc                                        >  4:197103/1‑149 (MQ=255)
                                                                                                              gAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATcc                                        <  2:103536/149‑1 (MQ=255)
                                                                                                                aGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCt                                      <  2:346183/149‑1 (MQ=255)
                                                                                                                   tttAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCATGTc                                   >  4:396720/1‑149 (MQ=255)
                                                                                                                     tAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCt                                 <  4:944435/149‑1 (MQ=255)
                                                                                                                     tAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCt                                 <  4:704948/149‑1 (MQ=255)
                                                                                                                     tAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCt                                 <  3:223230/149‑1 (MQ=255)
                                                                                                                         gTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGAtt                             <  1:420867/149‑1 (MQ=255)
                                                                                                                           aGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCAGATTcc                           <  2:781788/149‑1 (MQ=255)
                                                                                                                             gTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCAcccc                                                     >  4:683889/1‑121 (MQ=255)
                                                                                                                             gTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCAcccc                                                     <  3:683889/121‑1 (MQ=255)
                                                                                                                             gTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATcc                                        >  2:472215/1‑134 (MQ=255)
                                                                                                                             gTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATcc                                        <  1:472215/134‑1 (MQ=255)
                                                                                                                               aaCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTaa                       <  3:1059259/149‑1 (MQ=255)
                                                                                                                               aaCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTaa                       <  2:600360/149‑1 (MQ=255)
                                                                                                                                    ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGaaa                  >  1:790384/2‑149 (MQ=255)
                                                                                                                                    ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGaaa                  >  3:574497/2‑149 (MQ=255)
                                                                                                                                         gataaggcgTTCACGCCGCATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAacta             >  2:1074165/6‑149 (MQ=255)
                                                                                                                                         gATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAacta             >  2:397896/1‑149 (MQ=255)
                                                                                                                                                  tGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTaa                       >  2:522360/1‑130 (MQ=255)
                                                                                                                                                  tGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTaa                       <  1:522360/130‑1 (MQ=255)
                                                                                                                                                    aaCGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCa  >  3:879775/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCA  >  NC_000913/2842290‑2842586

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: