Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,797 T→G 5.0% N581H (AAC→CAC)  yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,7970TG 5.0% 1964.5 / 78.4 559N581H (AAC→CAC) yicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (285/246);  new base G (12/16);  total (297/262)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-02

Only 100 of 578 total aligned reads displayed.

GGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAA  >  NC_000913/3832650‑3832945
                                                                                                                                                   |                                                                                                                                                    
ggTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGtt                                                                                                                                                      <  3:215742/148‑1 (MQ=255)
ggTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTc                                                                                                                                                     >  2:1011842/1‑149 (MQ=255)
 gTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCg                                                                                                                                                    >  2:145125/1‑149 (MQ=255)
  tAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGtt                                                                                                                                                      >  4:234803/1‑146 (MQ=255)
  tAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCg                                                                                                                                                    >  3:4453/1‑148 (MQ=255)
   aGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCa                                                                                                                                                  >  1:812324/1‑149 (MQ=255)
      aCTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcgc                                                                                                                                               >  1:31754/1‑149 (MQ=255)
         gCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACgcgc                                                                                                                                             >  3:718003/1‑148 (MQ=255)
          cacaTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcgg                                                                                                                                           <  1:552560/149‑3 (MQ=255)
           acaTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcgg                                                                                                                                           <  3:860990/148‑3 (MQ=255)
               cGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAgggtac                                                                                                                                      <  1:510267/149‑6 (MQ=255)
                   cGCTTCAGTGGACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGa                                                                                                                                   >  3:832389/1‑148 (MQ=255)
                   cGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGa                                                                                                                                   >  4:610925/1‑148 (MQ=255)
                    gCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATa                                                                                                                                 >  2:42693/1‑149 (MQ=255)
                     cTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATAc                                                                                                                                <  2:31754/149‑1 (MQ=255)
                          gTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGAt                                                                                                                            >  4:855864/1‑148 (MQ=255)
                           tGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATa                                                                                                                           <  3:781986/148‑1 (MQ=255)
                           tGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATa                                                                                                                           >  4:1052229/1‑148 (MQ=255)
                           tGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATAc                                                                                                                          <  1:875491/149‑1 (MQ=255)
                            gAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCa                                                                                                                                      >  3:745368/1‑136 (MQ=255)
                                acCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGcat                                                                                                                     >  2:569113/1‑149 (MQ=255)
                                 ccGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcggg                                                                                                                                          <  1:894367/127‑4 (MQ=255)
                                    gcgcAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATg                                                                                                                 >  1:602027/1‑149 (MQ=255)
                                        aaCCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGc                                                                                                             <  1:650674/149‑1 (MQ=255)
                                            aTCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCAttt                                                                                                         >  1:223248/1‑149 (MQ=255)
                                              cACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTTCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGCTCGCACGCGCAGCTTCACGTTACAGATACGGCATCATGCGGCATTTc                                                                                                        >  4:694966/1‑148 (MQ=255)
                                              cACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTc                                                                                                        <  3:54636/148‑1 (MQ=255)
                                              cACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCa                                                                                                       <  1:1000915/149‑1 (MQ=255)
                                                   tGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGc                                                                                                  <  2:82306/149‑1 (MQ=255)
                                                     tCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGt                                                                                                <  1:347850/149‑1 (MQ=255)
                                                        ccTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGtt                                                                                                                                                      >  3:690680/1‑92 (MQ=255)
                                                           aaCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaagaa                                                                                          <  1:707675/149‑1 (MQ=255)
                                                           aaCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaagaa                                                                                          <  2:891425/149‑1 (MQ=255)
                                                           aaCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaagaa                                                                                          >  2:541416/1‑149 (MQ=255)
                                                            aCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaagaa                                                                                          <  3:636175/148‑1 (MQ=255)
                                                                 tATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGcac                                                                                    <  2:930586/149‑1 (MQ=255)
                                                                 tATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGca                                                                                     <  4:832389/148‑1 (MQ=255)
                                                                   ttGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGcacc                                                                                   >  4:464069/1‑148 (MQ=255)
                                                                     gACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCTGTTGCGTGAAGAAGCGCACcaca                                                                                >  1:170232/1‑149 (MQ=255)
                                                                      aCGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACcaca                                                                                >  4:318244/1‑148 (MQ=255)
                                                                      aCGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCACACcaca                                                                                <  3:178953/148‑1 (MQ=255)
                                                                         gTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcac                                                                            >  2:734881/1‑149 (MQ=255)
                                                                          tCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcaca                                                                           <  1:200532/149‑1 (MQ=255)
                                                                          tCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcac                                                                            >  4:520033/1‑148 (MQ=255)
                                                                              ggTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGc                                                                                                  >  2:355669/1‑122 (MQ=255)
                                                                              ggTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGa                                                                        <  4:1018442/148‑1 (MQ=255)
                                                                               gTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACt                                                                      >  2:180159/1‑149 (MQ=255)
                                                                                tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtc                                                                     <  1:408932/149‑1 (MQ=255)
                                                                                 aaTCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGcatc                                                                                                                    <  2:405669/101‑1 (MQ=255)
                                                                                 aaTCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtc                                                                     >  4:4900/1‑148 (MQ=255)
                                                                                 aaTCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTAGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtc                                                                     <  4:368400/148‑1 (MQ=255)
                                                                                   tCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGAGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcatc                                                                  <  1:569113/149‑1 (MQ=255)
                                                                                    cacaAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcat                                                                   >  1:731619/1‑147 (MQ=255)
                                                                                     acaAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcatca                                                                 >  3:14442/1‑148 (MQ=255)
                                                                                     acaAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcatca                                                                 >  3:336331/1‑148 (MQ=255)
                                                                                      caAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTc                                                                                                        <  4:383169/108‑1 (MQ=255)
                                                                                      caAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTcatcat                                                                >  3:245723/1‑148 (MQ=255)
                                                                                      caAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTcatcat                                                                >  3:884883/1‑148 (MQ=255)
                                                                                      caAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCatcgtc                                                               <  2:95828/149‑1 (MQ=255)
                                                                                       aaGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCg                                                              <  2:791761/149‑1 (MQ=255)
                                                                                       aaGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCg                                                              >  2:203236/1‑149 (MQ=255)
                                                                                        aGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt                                                             >  2:650874/1‑149 (MQ=255)
                                                                                         gCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa                                                            <  2:1057434/149‑1 (MQ=255)
                                                                                          ccGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcac                                                                            >  4:562743/1‑132 (MQ=255)
                                                                                           cGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAgg                                                          >  2:265272/1‑149 (MQ=255)
                                                                                             gATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGcc                                                        >  1:990674/1‑149 (MQ=255)
                                                                                             gATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGc                                                         >  4:151902/1‑148 (MQ=255)
                                                                                                   ccGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGTCCACgg                                                   >  3:1008712/1‑148 (MQ=255)
                                                                                                   ccGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACgg                                                   >  4:279215/1‑148 (MQ=255)
                                                                                                            ccATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATa                                          >  3:400149/1‑148 (MQ=255)
                                                                                                            ccATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATa                                          >  4:756986/1‑148 (MQ=255)
                                                                                                                catcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGAt                                      >  3:671273/1‑148 (MQ=255)
                                                                                                                 atcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGAttt                                    >  1:937488/1‑149 (MQ=255)
                                                                                                                 atcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGAtt                                     <  4:206185/148‑1 (MQ=255)
                                                                                                                  tcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTg                                   >  2:713966/1‑149 (MQ=255)
                                                                                                                   catGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGc                                  <  1:265272/149‑1 (MQ=255)
                                                                                                                     tGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCt                                 <  4:797833/148‑1 (MQ=255)
                                                                                                                         ccGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGtt                                                                                                    <  1:761397/77‑1 (MQ=255)
                                                                                                                          cGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCgtg                           <  1:1065689/149‑1 (MQ=255)
                                                                                                                            catcatCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTa                         >  2:114941/1‑149 (MQ=255)
                                                                                                                              tcatcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTaaa                       <  1:369652/149‑1 (MQ=255)
                                                                                                                               catcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAAc                      >  1:468223/1‑149 (MQ=255)
                                                                                                                               catcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAAc                      <  2:484030/149‑1 (MQ=255)
                                                                                                                                 tcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACgg                    >  1:667599/1‑149 (MQ=255)
                                                                                                                                 tcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACgg                    >  2:840412/1‑149 (MQ=255)
                                                                                                                                   ggCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGc                   <  3:1041317/148‑1 (MQ=255)
                                                                                                                                   ggCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGc                   <  3:375637/148‑1 (MQ=255)
                                                                                                                                       tACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATg               <  3:505851/148‑1 (MQ=255)
                                                                                                                                         ccccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCt            >  1:1032153/1‑149 (MQ=255)
                                                                                                                                           ctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTcatcat                                                                >  4:966890/3‑95 (MQ=255)
                                                                                                                                           ctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTg           >  4:715257/3‑148 (MQ=255)
                                                                                                                                            tgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCt            <  3:927638/145‑1 (MQ=255)
                                                                                                                                            tgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCt            >  4:927638/2‑146 (MQ=255)
                                                                                                                                            cgcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTgg          <  3:973442/148‑1 (MQ=255)
                                                                                                                                            cgcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTgg          >  3:28814/1‑148 (MQ=255)
                                                                                                                                            cgcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGga         <  2:627053/149‑1 (MQ=255)
                                                                                                                                              cgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGgaga       <  2:763526/149‑1 (MQ=255)
                                                                                                                                              cgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGgag        >  4:668039/1‑148 (MQ=255)
                                                                                                                                              cgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGgag        <  4:89185/148‑1 (MQ=255)
                                                                                                                                                   tCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGcaaa  <  2:672471/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAA  >  NC_000913/3832650‑3832945

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: