Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,535,978 G→T 33.6% D49Y (GAT→TAT)  crr → glucose‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,535,9780GT33.6% 523.7 / 311.6 315D49Y (GAT→TAT) crrglucose‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base G (102/107);  new base T (51/55);  total (153/162)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

GATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGA  >  NC_000913/2535830‑2536126
                                                                                                                                                    |                                                                                                                                                    
gaTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTg                                                                                                                                                      >  3:253535/1‑149 (MQ=255)
 aTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGa                                                                                                                                                     >  3:1747934/1‑149 (MQ=255)
 aTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGa                                                                                                                                                     >  4:795437/1‑149 (MQ=255)
      gggTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTa                                                                                                                                                >  2:2402787/1‑149 (MQ=255)
          ttGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGc                                                                                                                                            <  1:1840773/149‑1 (MQ=255)
            gTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTTTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTa                                                                                                                                          <  2:1814899/149‑1 (MQ=255)
            gTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTa                                                                                                                                          >  3:1183915/1‑149 (MQ=255)
            gTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTa                                                                                                                                          <  4:1053461/149‑1 (MQ=255)
                gATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCaa                                                                                                                                      >  4:629393/1‑149 (MQ=255)
                  tAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAAc                                                                                                                                    >  4:1291985/1‑149 (MQ=255)
                    aaCTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCa                                                                                                                                  >  4:1248225/1‑149 (MQ=255)
                      cTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAAc                                                                                                                                <  4:1163346/149‑1 (MQ=255)
                      cTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAAc                                                                                                                                <  3:1849366/149‑1 (MQ=255)
                      cTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAAc                                                                                                                                >  3:287483/1‑149 (MQ=255)
                        gAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAACCGTTGGTGCTGGTATTGCTATCAAACCAACgg                                                                                                                              <  4:1512064/149‑1 (MQ=255)
                          aaTCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGt                                                                                                                            <  1:329099/149‑1 (MQ=255)
                           aTCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTa                                                                                                                           >  1:2012583/1‑149 (MQ=255)
                           aTCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTa                                                                                                                           >  2:60876/1‑149 (MQ=255)
                            tctcTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTaa                                                                                                                          >  2:1447491/1‑149 (MQ=255)
                              tctGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACa                                                                                                                        >  2:1057857/1‑149 (MQ=255)
                                 ggTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAAACAACGGGTAACaaaa                                                                                                                     <  4:1834133/149‑1 (MQ=255)
                                  gTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAAt                                                                                                                    <  3:629393/149‑1 (MQ=255)
                                  gTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAAt                                                                                                                    <  1:511022/149‑1 (MQ=255)
                                      ccGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTc                                                                                                                <  1:2402787/149‑1 (MQ=255)
                                       cgacgaCAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACaaaa                                                                                                                     >  4:1421287/1‑143 (MQ=255)
                                             cAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAg                                                                                                         >  1:1940996/1‑149 (MQ=255)
                                              aagaagGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGt                                                                                                        >  4:168819/1‑149 (MQ=255)
                                               agaagGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTa                                                                                                       <  2:1940996/149‑1 (MQ=255)
                                               agaagGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTa                                                                                                       <  4:1127267/149‑1 (MQ=255)
                                                gaagGATACCGGAACTATTGAGATCATTACTGCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGCTTTTGCGGGAAAAATCGGTAGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAg                                                                                                      >  4:1523949/1‑149 (MQ=255)
                                                    gATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACgg                                                                                                  <  2:1526141/149‑1 (MQ=255)
                                                    gATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACgg                                                                                                  <  1:1615257/149‑1 (MQ=255)
                                                       aCCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTAGTGCTATCAAACCAACGGGTAACAGAATGGTCGCGCCAGTAGACGGCAc                                                                                               <  1:700561/149‑1 (MQ=255)
                                                        ccGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCAcc                                                                                              >  1:2473107/1‑149 (MQ=255)
                                                            aaCTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTg                                                                                          >  4:2385600/1‑149 (MQ=255)
                                                               tATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTa                                                                                       >  4:1395393/1‑149 (MQ=255)
                                                                aTTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTaa                                                                                      <  2:2098692/149‑1 (MQ=255)
                                                                 ttGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTaaa                                                                                     >  1:1001427/1‑149 (MQ=255)
                                                                 ttGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTaaa                                                                                     >  3:71830/1‑149 (MQ=255)
                                                                 ttGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTaaa                                                                                     >  3:1384929/1‑149 (MQ=255)
                                                                    agaTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACgg                                                                                                                              <  2:2208259/105‑1 (MQ=255)
                                                                    agaTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACgg                                                                                                                              >  1:2208259/1‑105 (MQ=255)
                                                                    agaTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATc                                                                                                                                        >  3:439930/1‑95 (MQ=255)
                                                                    agaTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATc                                                                                                                                        <  4:439930/95‑1 (MQ=255)
                                                                      aTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACg                                                                                                                               <  2:322001/102‑1 (MQ=255)
                                                                      aTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCAt                                                                                            <  1:1444332/137‑1 (MQ=255)
                                                                       tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                               <  1:63377/149‑1 (MQ=255)
                                                                       tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                               >  4:1751372/1‑149 (MQ=255)
                                                                       tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                               <  3:958419/149‑1 (MQ=255)
                                                                       tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                               <  4:2365013/149‑1 (MQ=255)
                                                                       tCATTGCTCCGCTCTCTGGCGAGATCGGCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                               <  1:2336924/149‑1 (MQ=255)
                                                                           tGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGaaa                                                                           <  2:135546/149‑1 (MQ=255)
                                                                             cTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAaacc                                                                         >  3:1114686/1‑149 (MQ=255)
                                                                              tCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAaacca                                                                        >  1:1292750/1‑149 (MQ=255)
                                                                                cGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAaaccaac                                                                      >  4:1858724/1‑149 (MQ=255)
                                                                                  ctctctGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAaccaacca                                                                    <  3:1248225/149‑1 (MQ=255)
                                                                                    ctctGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACg                                                                  >  3:1716435/1‑149 (MQ=255)
                                                                                      ctGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGaaa                                                                           >  1:788851/1‑138 (MQ=255)
                                                                                        ggCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCAtt                                                              >  1:1123404/1‑149 (MQ=255)
                                                                                        ggCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCAtt                                                              >  4:2308360/1‑149 (MQ=255)
                                                                                          cGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATtct                                                            <  3:795437/149‑1 (MQ=255)
                                                                                           gagaTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATtctc                                                           <  1:11715/149‑1 (MQ=255)
                                                                                           gagaTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATtctc                                                           >  1:2311327/1‑149 (MQ=255)
                                                                                             gaTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTa                                                         <  4:274038/149‑1 (MQ=255)
                                                                                             gaTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTAACAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTa                                                         >  1:953866/1‑149 (MQ=255)
                                                                                              aTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTAt                                                        <  4:1193363/149‑1 (MQ=255)
                                                                                               tCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATc                                                       <  2:228887/149‑1 (MQ=255)
                                                                                                 gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTACCGa                                                     >  3:1500893/1‑149 (MQ=255)
                                                                                                 gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGa                                                     >  2:1766268/1‑149 (MQ=255)
                                                                                                 gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGa                                                     >  2:1685586/1‑149 (MQ=255)
                                                                                                   cAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAAt                                                   <  2:1090196/149‑1 (MQ=255)
                                                                                                     atatCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCt                                                 >  2:1276153/1‑149 (MQ=255)
                                                                                                     atatCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCt                                                                                                                                           <  4:2004280/59‑1 (MQ=255)
                                                                                                     atatCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCt                                                                                                                                           >  3:2004280/1‑59 (MQ=255)
                                                                                                     atatCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCt                                                 <  1:390272/149‑1 (MQ=255)
                                                                                                         cGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATa                                             <  3:1417641/149‑1 (MQ=255)
                                                                                                            aGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgcg                                          >  3:1586743/1‑149 (MQ=255)
                                                                                                              aCGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgcggc                                        <  2:2473107/149‑1 (MQ=255)
                                                                                                                 tGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGtt                                     <  2:1225916/149‑1 (MQ=255)
                                                                                                                 tGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGtt                                     <  1:1775845/149‑1 (MQ=255)
                                                                                                                   ccGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                               <  2:2351779/105‑1 (MQ=255)
                                                                                                                   ccGGATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGa                                   >  2:1339138/1‑149 (MQ=255)
                                                                                                                   ccGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGa                                   <  4:211011/149‑1 (MQ=255)
                                                                                                                     ggacggcGTTTTTGCGGAAAAAATCGTGGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAAc                                 <  4:1603298/143‑1 (MQ=255)
                                                                                                                      gATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATtctc                                                           <  4:55400/122‑1 (MQ=255)
                                                                                                                      gATGTCGTTTTTGCGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATtctc                                                           >  3:55400/1‑122 (MQ=255)
                                                                                                                          tCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTc                            <  1:1057857/149‑1 (MQ=255)
                                                                                                                           cGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCg                           <  3:1802876/149‑1 (MQ=255)
                                                                                                                                ttGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCAc                      <  3:2354956/149‑1 (MQ=255)
                                                                                                                                 tGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTaa                                                                                      >  2:1010499/1‑84 (MQ=255)
                                                                                                                                   cGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTc                   >  3:2153456/1‑149 (MQ=255)
                                                                                                                                   cGGAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTc                   <  4:2470262/149‑1 (MQ=255)
                                                                                                                                   cGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTc                   <  2:1177645/149‑1 (MQ=255)
                                                                                                                                     gAAAAAATCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTCTCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCgg                 >  2:782162/1‑149 (MQ=255)
                                                                                                                                           aTCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGa           >  2:98596/1‑149 (MQ=255)
                                                                                                                                            tCGTTGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGac          >  4:916697/1‑149 (MQ=255)
                                                                                                                                              gTTGGTGATGGTATTGATATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTCGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGacac        >  4:739808/1‑149 (MQ=255)
                                                                                                                                               ttGGTTATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACAACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACAcc       >  2:68789/1‑149 (MQ=255)
                                                                                                                                                  gTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGtt    <  2:851341/149‑1 (MQ=255)
                                                                                                                                                    gATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGa  >  1:2356362/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGA  >  NC_000913/2535830‑2536126

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