Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,397,825 C→A 5.9% V266V (GTG→GTT uspE ← stress‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,397,8250CA 5.9% 219.7 / 2.6 68V266V (GTG→GTTuspEstress‑induced protein
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base A (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

CTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTA  >  NC_000913/1397680‑1397972
                                                                                                                                                 |                                                                                                                                                   
cttcttCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACcac                                                                                                                                                   >  4:105513/1‑148 (MQ=255)
   cttcGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATc                                                                                                                                               >  1:267130/1‑149 (MQ=255)
      cGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACca                                                                                                                                                    >  3:44502/1‑141 (MQ=255)
      cGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACca                                                                                                                                                    <  4:44502/141‑1 (MQ=255)
        tCATCCAGTTCAACGGGTGGCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGc                                                                                                                                          <  6:230370/149‑1 (MQ=255)
           tCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGAGCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTg                                                                                                                                       <  3:7600/149‑1 (MQ=255)
              aGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGt                                                                                                                                      >  2:241672/1‑147 (MQ=255)
               gTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGAGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAg                                                                                                                                   >  3:241784/1‑149 (MQ=255)
                    aCGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTc                                                                                                                              >  3:86845/1‑149 (MQ=255)
                     cGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTcc                                                                                                                             <  4:115437/149‑1 (MQ=255)
                        gtgtCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGcc                                                                                                                          <  2:267130/149‑1 (MQ=255)
                         tgtCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCa                                                                                                                         <  7:313013/149‑1 (MQ=255)
                           tCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCaaa                                                                                                                       >  8:213371/1‑149 (MQ=255)
                             tGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAAt                                                                                                                      >  5:277235/1‑148 (MQ=255)
                              gATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCaaagcc                                                                                                                    >  6:137066/1‑146 (MQ=255)
                                 aCTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGTACAACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAgg                                                                                                                  >  5:87357/1‑148 (MQ=255)
                                 aCTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACAACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAgg                                                                                                                  >  3:225567/1‑148 (MQ=255)
                                  cTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGa                                                                                                                 <  6:302012/148‑1 (MQ=255)
                                      tCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCa                                                                                                             >  7:202085/1‑148 (MQ=255)
                                          gTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAg                                                                                                                   <  4:200761/138‑1 (MQ=255)
                                          gTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAg                                                                                                                   >  3:200761/1‑138 (MQ=255)
                                           tttAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCtctt                                                                                                       >  2:18767/1‑149 (MQ=255)
                                              aaTGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCtcttct                                                                                                     <  4:71454/148‑1 (MQ=255)
                                               aTGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGTACAACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTgg                                                                                                   <  6:87357/149‑1 (MQ=255)
                                               aTGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTgg                                                                                                   <  7:145752/149‑1 (MQ=255)
                                               aTGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACAACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTgg                                                                                                   <  4:225567/149‑1 (MQ=255)
                                                tGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTgg                                                                                                   >  6:318975/1‑148 (MQ=255)
                                                tGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGCTCCGCCGTGTCGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGAGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGc                                                                                                  <  5:148661/149‑1 (MQ=255)
                                                  aCCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCa                                                                                                 <  3:105513/148‑1 (MQ=255)
                                                      gcagGTCGCAGCGAAGATGATCAATCCCCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCt                                                                                            >  4:113422/1‑149 (MQ=255)
                                                          gTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCttttt                                                                                        <  6:277235/149‑1 (MQ=255)
                                                          gTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCttttt                                                                                        >  8:198336/1‑149 (MQ=255)
                                                               aGCGAAGATTAGCAATCACCTGTTCCGCCGTGTTGGCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAAGCACCTCTTCTGGCAGACCTTTTTCt                                                                                      <  4:313758/146‑1 (MQ=255)
                                                                gCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGt                                                                                  <  3:158530/149‑1 (MQ=255)
                                                                        gatCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTTCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCagg                                                                          >  2:227584/1‑147 (MQ=255)
                                                                        gatCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATg                                                                          >  4:227289/1‑149 (MQ=255)
                                                                        gatCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATg                                                                          >  3:231197/1‑149 (MQ=255)
                                                                         atcaatcaCCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTACGCCAAATCAGGAATCACCTCTTCTGGCATACCTTTTTCTACGTGCGTCATg                                                                          <  1:242551/148‑1 (MQ=255)
                                                                              tcaCCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGtttt                                                                      <  2:31922/147‑1 (MQ=255)
                                                                              tcaCCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGtttt                                                                      >  1:31922/1‑147 (MQ=255)
                                                                              tcaCCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCa                                                                    <  5:318975/149‑1 (MQ=255)
                                                                              tcaCCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCa                                                                    <  7:213371/149‑1 (MQ=255)
                                                                                aCCTGTTCCGCCGTGTTGCCGAGGACTGCTGCTGAAATACCGGTGCGGCCTACCGTGCCCAGACCCACAATCCCCGCCTGTAACTGCGCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCAtt                                                                  <  8:243832/149‑1 (MQ=255)
                                                                                  cTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTaa                                                                <  3:167785/149‑1 (MQ=255)
                                                                                        cgccgTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCg                                                           >  2:72018/1‑148 (MQ=255)
                                                                                         gccgTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGaa                                                         >  4:31264/1‑149 (MQ=255)
                                                                                          ccgTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAAt                                                        <  8:202085/149‑1 (MQ=255)
                                                                                            gtgtTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAAttt                                                      <  5:113353/149‑1 (MQ=255)
                                                                                               ttGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTgg                                                                                                   <  1:60912/101‑1 (MQ=255)
                                                                                               ttGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTgg                                                                                                   >  2:60912/1‑101 (MQ=255)
                                                                                               ttGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTg                                                   <  4:86845/149‑1 (MQ=255)
                                                                                                    gAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCg                                                >  1:199311/1‑147 (MQ=255)
                                                                                                          tgctgcTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAAGCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCTTTAATGCCGAATTTCTGCCGCAGGGCtt                                        >  6:155586/1‑149 (MQ=255)
                                                                                                               ctgAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCa                                     <  1:18767/147‑1 (MQ=255)
                                                                                                                tgAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTg                                  >  6:174023/1‑149 (MQ=255)
                                                                                                                 gAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGc                                                                                                                           <  5:333102/59‑1 (MQ=255)
                                                                                                                 gAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGc                                                                                                                           >  6:333102/1‑59 (MQ=255)
                                                                                                                   aaTACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCa                               <  7:198336/149‑1 (MQ=255)
                                                                                                                     tACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAg                              >  3:289302/1‑148 (MQ=255)
                                                                                                                        cGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCaaa                          >  7:143791/1‑149 (MQ=255)
                                                                                                                        cGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCaaa                          >  2:18017/1‑149 (MQ=255)
                                                                                                                                          ccAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGCATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTATCATTAATGCCGAATTTCTGCCGCAGGGCTTACATAGCCCGCAAATCTTGCCCACGAATGGCa        >  8:293630/1‑149 (MQ=255)
                                                                                                                                          ccAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCa        >  6:142810/1‑149 (MQ=255)
                                                                                                                                           cAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCAt       >  2:46107/1‑149 (MQ=255)
                                                                                                                                           cAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCAt       >  3:46209/1‑149 (MQ=255)
                                                                                                                                            aGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAATTCAGGAATCACCTCTTCTGGCAGACCTTTTGCTACTTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATc      <  5:137066/149‑1 (MQ=255)
                                                                                                                                            aGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATg                        <  1:191165/131‑1 (MQ=255)
                                                                                                                                            aGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATg                        >  2:191165/1‑131 (MQ=255)
                                                                                                                                            aGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATc      <  5:174023/149‑1 (MQ=255)
                                                                                                                                            aGAACCACAATCCCCGCCTGTAAGTGCGCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTAATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATc      <  4:241784/149‑1 (MQ=255)
                                                                                                                                             gAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCg     <  4:46209/149‑1 (MQ=255)
                                                                                                                                                 cacaATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTta  >  3:323781/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTTAATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCCGCAGGGCTTTCATTGCCAGCAAATGTTGCCCACGAATGGCATCGTTA  >  NC_000913/1397680‑1397972

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