Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,952,587 T→A 5.8% W108R (TGG→AGG)  yecN → MAPEG family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,952,5870TA 5.8% 538.6 / 14.3 152W108R (TGG→AGG) yecNMAPEG family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (67/76);  new base A (4/5);  total (71/81)
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 153 total aligned reads displayed.

TGCTGATGCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTA  >  NC_000913/1952441‑1952735
                                                                                                                                                   |                                                                                                                                                    
tgCTGATGCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTgg                                                                                                                                                    <  1:601928/149‑1 (MQ=255)
  cTGATGCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGgtg                                                                                                                                                  >  3:114503/1‑149 (MQ=255)
   tGATGCTGATTATGGAAATGAATGATGTAGAAGCGTGGATGGTGCATATTTGCGGCATCGTTTTGCTAACTGGTTGACTGATGCATTATTCCGGTTTGCCTCACCGTCTGTTCCGCTTGCGGCGTTACGGCAGGAGCGCCTGCTGgtgt                                                                                                                                                 <  4:1382593/149‑1 (MQ=255)
    gATGCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGtgtg                                                                                                                                                >  1:118029/1‑149 (MQ=255)
     aTGCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTgc                                                                                                                                               <  1:1437680/149‑1 (MQ=255)
       gCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCACACCAGgtggcgc                                                                                                                                              <  4:1226513/149‑5 (MQ=255)
        cTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCt                                                                                                                                            >  4:506800/1‑149 (MQ=255)
           tttATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGtt                                                                                                                                         <  1:885540/149‑1 (MQ=255)
             tATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCACACCAGGTGGCGCtcatgc                                                                                                                                       <  4:120805/149‑6 (MQ=255)
             tATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCACACCAGGTGGCGCtcatgc                                                                                                                                       <  4:233678/149‑6 (MQ=255)
                 gAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCACACCAGgtggcgc                                                                                                                                              <  3:1443175/139‑5 (MQ=255)
                 gAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCACACCAGgtggcgc                                                                                                                                              >  4:1443175/1‑135 (MQ=255)
                  aaaTGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATg                                                                                                                                  <  2:213514/149‑1 (MQ=255)
                  aaaTGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTGTTGCTTGCTGGTCGTCTGAGGCATTATTACGGTTTTGATCACCGTCTGTTCCGCTGGCGACGGTCCGGCATGAGCGCCACCTGGCGTGCGCTGTTGCTGATg                                                                                                                                  >  4:473026/1‑149 (MQ=255)
                   aatgaatgGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATgg                                                                                                                                 <  2:887829/149‑1 (MQ=255)
                     agaatgGCGCAGAAACCCGGATGGTACATATTGGCGGCATCGTGTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTg                                                                                                                               >  4:32755/2‑149 (MQ=255)
                          ggCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGc                                                                                                                          <  1:342647/149‑1 (MQ=255)
                             gcAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGaa                                                                                                                       >  4:479099/1‑149 (MQ=255)
                              cAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAAt                                                                                                                      <  2:304019/149‑1 (MQ=255)
                               aGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATc                                                                                                                     >  2:239264/1‑149 (MQ=255)
                                gAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCt                                                                                                                    >  3:496971/1‑149 (MQ=255)
                                 aaaCCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCtt                                                                                                                   <  1:66753/149‑1 (MQ=255)
                                  aaCCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCttt                                                                                                                  <  1:5902/149‑1 (MQ=255)
                                  aaCCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCttt                                                                                                                  <  2:1296146/149‑1 (MQ=255)
                                  aaCCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCttt                                                                                                                  <  3:25195/149‑1 (MQ=255)
                                   aCCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTg                                                                                                                 >  3:1055447/1‑149 (MQ=255)
                                   aCCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTg                                                                                                                 <  1:990165/149‑1 (MQ=255)
                                    ccTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTgg                                                                                                                <  3:1152921/149‑1 (MQ=255)
                                           ggTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGc                                                                                                         <  1:770895/149‑1 (MQ=255)
                                             tGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGccc                                                                                                       <  4:1205590/149‑1 (MQ=255)
                                              gCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGa                                                                                                                                    >  2:1415408/1‑119 (MQ=255)
                                              gCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCt                                                                                                      >  3:1084142/1‑149 (MQ=255)
                                               cATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTg                                                                                                     <  3:180819/149‑1 (MQ=255)
                                                atatTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTgg                                                                                                    <  4:68257/149‑1 (MQ=255)
                                                atatTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTgg                                                                                                    >  1:279183/1‑149 (MQ=255)
                                                 tatTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTggg                                                                                                   >  2:1014971/1‑149 (MQ=255)
                                                  atTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGa                                                                                                  <  1:974722/149‑1 (MQ=255)
                                                  atTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGa                                                                                                  <  1:231440/149‑1 (MQ=255)
                                                   tttGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAg                                                                                                 >  1:895194/1‑149 (MQ=255)
                                                   tttGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAg                                                                                                 >  4:53608/1‑149 (MQ=255)
                                                     tGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGtt                                                                                               >  2:124127/1‑149 (MQ=255)
                                                      gCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTg                                                                                              <  1:1014971/149‑1 (MQ=255)
                                                      gCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTCCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTg                                                                                              <  2:619559/149‑1 (MQ=255)
                                                        ggCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCACACCAGgtggcgc                                                                                                                                              >  2:531125/1‑96 (MQ=255)
                                                        ggCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGt                                                                                            <  1:10324/149‑1 (MQ=255)
                                                         gCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGtt                                                                                           <  1:1444747/149‑1 (MQ=255)
                                                         gCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGtt                                                                                           >  1:374992/1‑149 (MQ=255)
                                                            tCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTgc                                                                                                                                               >  2:614935/1‑94 (MQ=255)
                                                             cGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTtct                                                                                       <  3:32755/149‑1 (MQ=255)
                                                              gTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTtctc                                                                                      >  4:616641/1‑149 (MQ=255)
                                                               ttttGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTcc                                                                                     <  3:224184/149‑1 (MQ=255)
                                                                   gcttgctGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGc                                                                                 <  3:566795/149‑1 (MQ=255)
                                                                     ttgctGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGt                                                                               >  2:406727/1‑149 (MQ=255)
                                                                         tGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAgc                                                                           >  1:1258094/1‑149 (MQ=255)
                                                                         tGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAgc                                                                           >  1:334325/1‑149 (MQ=255)
                                                                           gtcgtcTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGt                                                                                            <  2:468352/130‑1 (MQ=255)
                                                                           gtcgtcTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGt                                                                                            >  1:468352/1‑130 (MQ=255)
                                                                           gtcgtcTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAgcgc                                                                         >  1:395756/1‑149 (MQ=255)
                                                                            tcgtcTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGca                                                                        <  2:334325/149‑1 (MQ=255)
                                                                                    tGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCTTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGc                                                                <  3:66526/149‑1 (MQ=255)
                                                                                    tGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCAGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGGCCTGGGAGCTGGTTTTCTCCCTGCGTTAGCGCACAATACGc                                                                >  4:1191181/1‑149 (MQ=255)
                                                                                      cATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCa                                                              <  1:29337/149‑1 (MQ=255)
                                                                                       attattACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTg                                                                                                                 <  3:337115/97‑1 (MQ=255)
                                                                                        ttattaCGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCAct                                                            >  1:1372842/1‑149 (MQ=255)
                                                                                            taCGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCActttct                                                        >  3:125781/1‑149 (MQ=255)
                                                                                            taCGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCActttct                                                        >  4:377089/1‑149 (MQ=255)
                                                                                              cGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCActttcttt                                                      <  1:1062043/149‑1 (MQ=255)
                                                                                               ggTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGt                                                                                            >  2:494911/1‑110 (MQ=255)
                                                                                               ggTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGt                                                                                            <  1:494911/110‑1 (MQ=255)
                                                                                               ggTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCtttt                                                     <  3:1263452/149‑1 (MQ=255)
                                                                                                gTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCttttt                                                    <  4:835751/149‑1 (MQ=255)
                                                                                                 ttttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTc                                                   >  3:475868/1‑149 (MQ=255)
                                                                                                 ttttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTc                                                   >  3:1385066/1‑149 (MQ=255)
                                                                                                 ttttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTc                                                   >  2:336511/1‑149 (MQ=255)
                                                                                                      atcaCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGa                                              <  2:395756/149‑1 (MQ=255)
                                                                                                      atcaCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGa                                              <  2:485185/149‑1 (MQ=255)
                                                                                                       tcaCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGAt                                             >  4:1338250/1‑149 (MQ=255)
                                                                                                        caCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGAGGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGAtt                                            >  4:1127728/1‑149 (MQ=255)
                                                                                                             tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTAc                                       <  1:239264/149‑1 (MQ=255)
                                                                                                             tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTAc                                       >  3:576979/1‑149 (MQ=255)
                                                                                                                gTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGtt                                    <  3:1338250/149‑1 (MQ=255)
                                                                                                                 ttCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTGTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGACTTTTACGTTa                                   >  3:576621/1‑149 (MQ=255)
                                                                                                                 ttCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTa                                   <  1:435889/149‑1 (MQ=255)
                                                                                                                 ttCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTa                                   >  4:459823/1‑149 (MQ=255)
                                                                                                                  tCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTAt                                  <  4:1084142/149‑1 (MQ=255)
                                                                                                                   ccGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATg                                 <  4:571916/149‑1 (MQ=255)
                                                                                                                           gACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGGGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATg                                 <  2:813979/141‑1 (MQ=255)
                                                                                                                              gTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGa                      <  2:279183/149‑1 (MQ=255)
                                                                                                                               ttCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGac                     >  1:65091/1‑149 (MQ=255)
                                                                                                                                 ccGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGATTCTTTGGTATATTCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGTACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGacac                   >  4:1069364/1‑149 (MQ=255)
                                                                                                                                   ggCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGc                 >  2:365662/1‑149 (MQ=255)
                                                                                                                                    gCATGAGCGCCACCTGGTGTGCACTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTTTGTCTCACCGCGACACGCt                >  2:1127546/1‑149 (MQ=255)
                                                                                                                                       tGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTAtt             >  4:1028214/1‑149 (MQ=255)
                                                                                                                                        cagcgcacACCAGGTGGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGTCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTAttt            >  3:777538/8‑149 (MQ=255)
                                                                                                                                           cgcCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCt         <  2:271019/149‑1 (MQ=255)
                                                                                                                                             ccACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGc       >  3:945018/1‑149 (MQ=255)
                                                                                                                                             ccACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGc       >  1:768771/1‑149 (MQ=255)
                                                                                                                                               cACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGcc      <  4:907166/149‑1 (MQ=255)
                                                                                                                                                  cTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTAc                                       >  2:760024/1‑113 (MQ=255)
                                                                                                                                                   tGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTa  >  4:1107632/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGCTGATGCTGTTTATGGAAATGAATGGCGCAGAAACCTGGATGGTGCATATTTGCGGCATCGTTTTGCTTGCTGGTCGTCTGATGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTA  >  NC_000913/1952441‑1952735

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