Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,754,924 G→T 7.1% intergenic (‑496/+49) yiaW ← / ← aldB DUF3302 family inner membrane protein/aldehyde dehydrogenase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,754,9240GT 7.1% 440.5 / 17.8 126intergenic (‑496/+49)yiaW/aldBDUF3302 family inner membrane protein/aldehyde dehydrogenase B
Reads supporting (aligned to +/- strand):  ref base G (48/69);  new base T (4/5);  total (52/74)
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.63e-01

Only 100 of 129 total aligned reads displayed.

AAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCATCATCTTGTGGGTTTCGCGACCGATAC  >  NC_000913/3754778‑3755070
                                                                                                                                                  |                                                                                                                                                  
aaaCGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGt                                                                                                                                                   <  3:108773/148‑1 (MQ=255)
   cGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAg                                                                                                                                                <  3:111829/148‑1 (MQ=255)
   cGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGc                                                                                                                                               <  2:375268/149‑1 (MQ=255)
     tACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCg                                                                                                                                              >  3:186432/1‑148 (MQ=255)
      aCTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCTTACGCGt                                                                                                                                             <  3:376441/148‑1 (MQ=255)
      aCTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTc                                                                                                                                            <  2:21257/149‑1 (MQ=255)
       cTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCTTACGCGTc                                                                                                                                            <  4:37086/148‑1 (MQ=255)
        ttCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGc                                                                                                                                          <  1:42146/149‑1 (MQ=255)
            gTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAGACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATc                                                                                                                                       <  3:150182/148‑1 (MQ=255)
            gTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATc                                                                                                                                       >  4:61412/1‑148 (MQ=255)
                  gTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAAt                                                                                                                                <  1:181846/149‑1 (MQ=255)
                  gTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAAt                                                                                                                                <  2:6546/149‑1 (MQ=255)
                  gTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAAt                                                                                                                                <  1:330211/149‑1 (MQ=255)
                   tAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATGCGACGCTTACGCGTctta                                                                                                                                         <  4:235145/139‑4 (MQ=38)
                   tAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATGCGACGCTTACGCGTctta                                                                                                                                         >  3:235145/1‑136 (MQ=38)
                       cTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATa                                                                                                                           <  2:196150/149‑1 (MQ=255)
                        tGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATa                                                                                                                           <  4:232733/148‑1 (MQ=255)
                        tGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATa                                                                                                                           <  4:312271/148‑1 (MQ=255)
                        tGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATAc                                                                                                                          >  2:277080/1‑149 (MQ=255)
                                  cgcACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGAc                                                                                                                <  1:69453/149‑1 (MQ=255)
                                    cACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCa                                                                                                                                  <  3:113505/129‑1 (MQ=255)
                                    cACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATAcc                                                                                                                         >  3:357131/1‑138 (MQ=255)
                                     aCGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCa                                                                                                              <  4:259250/148‑1 (MQ=255)
                                         aGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGcttc                                                                                                         <  2:30196/149‑1 (MQ=255)
                                         aGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGcttc                                                                                                         >  2:149650/1‑149 (MQ=255)
                                         aGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGctt                                                                                                          <  3:61412/148‑1 (MQ=255)
                                         aGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGctt                                                                                                          <  4:186432/148‑1 (MQ=255)
                                          ggTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGcttct                                                                                                        <  2:263640/149‑1 (MQ=255)
                                            tttAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGt                                                                                                                                             >  4:268098/1‑110 (MQ=255)
                                            tttAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGt                                                                                                                                             <  3:268098/110‑1 (MQ=255)
                                            tttAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGTCAATGAATACCCAATGCGACCAGCTTCTta                                                                                                      <  1:326141/149‑1 (MQ=255)
                                            tttAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGcttctt                                                                                                       <  4:371992/148‑1 (MQ=255)
                                            tttAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTta                                                                                                      >  2:133819/1‑149 (MQ=255)
                                             ttAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTtat                                                                                                     <  2:21187/149‑1 (MQ=255)
                                             ttAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTta                                                                                                      >  4:130899/1‑148 (MQ=255)
                                               aaaaCGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTatat                                                                                                   <  2:302394/149‑1 (MQ=255)
                                               aaaaCGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTatat                                                                                                   <  1:326729/149‑1 (MQ=255)
                                                aaaCGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTatat                                                                                                   <  4:221134/148‑1 (MQ=255)
                                                aaaCGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTatat                                                                                                   >  4:137876/1‑148 (MQ=255)
                                                    gTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGaa                                                                                              >  1:17729/1‑149 (MQ=255)
                                                     tAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGaa                                                                                              >  3:336787/1‑148 (MQ=255)
                                                     tAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCAAGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGaa                                                                                              >  4:229395/1‑148 (MQ=255)
                                                      aaaTTGTAACAACGTGCGTTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAAc                                                                                             <  4:180088/148‑1 (MQ=255)
                                                         ttGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGcc                                                                                         >  2:352420/1‑149 (MQ=255)
                                                           gTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGcccc                                                                                       >  1:352630/1‑149 (MQ=255)
                                                                aacGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACgg                                                                                  >  2:114779/1‑149 (MQ=255)
                                                                 acGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACgg                                                                                  >  3:218202/1‑148 (MQ=255)
                                                                    tGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGttt                                                                               <  3:26743/148‑1 (MQ=255)
                                                                       gcTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATcc                                                                           <  2:77227/149‑1 (MQ=255)
                                                                         tttGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGa                                                                         <  2:197209/149‑1 (MQ=255)
                                                                         tttGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGa                                                                         <  1:266645/149‑1 (MQ=255)
                                                                          ttGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGa                                                                         >  3:371482/1‑148 (MQ=255)
                                                                          ttGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGa                                                                         >  3:145355/1‑148 (MQ=255)
                                                                              ttATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGa                                                                         >  3:76115/1‑144 (MQ=255)
                                                                              ttATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGa                                                                         <  4:76115/144‑1 (MQ=255)
                                                                                aTGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTc                                                                  >  2:125669/1‑149 (MQ=255)
                                                                                     ggATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACca                                                              <  4:119624/148‑1 (MQ=255)
                                                                                       aTGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAg                                                                                                                                                <  4:281408/64‑1 (MQ=255)
                                                                                       aTGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAg                                                                                                                                                >  3:281408/1‑64 (MQ=255)
                                                                                        taaggcgTACACGCCGCATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACcagca                                                           >  3:187565/4‑148 (MQ=39)
                                                                                        tGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACcagcag                                                          <  1:173080/149‑1 (MQ=255)
                                                                                            gcgTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCa                                                       <  4:329428/148‑1 (MQ=255)
                                                                                                 tACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGt                                                 <  1:125286/149‑1 (MQ=255)
                                                                                                  aCGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATTTCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGt                                                 <  3:34898/148‑1 (MQ=255)
                                                                                                    gCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTg                                              >  1:155898/1‑149 (MQ=255)
                                                                                                      cTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCa                                                                                                 >  2:48480/1‑96 (MQ=255)
                                                                                                      cTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCa                                                                                                 <  1:48480/96‑1 (MQ=255)
                                                                                                      cTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGc                                             <  4:336683/148‑1 (MQ=255)
                                                                                                      cTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCt                                            <  1:114779/149‑1 (MQ=255)
                                                                                                      cTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCt                                            <  1:91842/149‑1 (MQ=255)
                                                                                                       ttATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATAc                                                                                                                          <  4:21407/70‑1 (MQ=255)
                                                                                                       ttATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATAc                                                                                                                          >  3:21407/1‑70 (MQ=255)
                                                                                                       ttATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCAACAGCAGGCACTTGGTTTGCTg                                           >  1:170725/1‑149 (MQ=255)
                                                                                                        tATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTgg                                          <  2:118648/149‑1 (MQ=255)
                                                                                                         aTCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGt                                         <  1:358243/149‑1 (MQ=255)
                                                                                                          tCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGt                                         <  4:371482/148‑1 (MQ=255)
                                                                                                           ccGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTaa                                       <  2:17729/149‑1 (MQ=255)
                                                                                                             ggCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACtt                                                    <  4:239947/134‑1 (MQ=255)
                                                                                                             ggCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAAt                                      <  4:236686/148‑1 (MQ=255)
                                                                                                                  aCGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCa                                <  1:131437/149‑1 (MQ=255)
                                                                                                                     gTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGc                                    <  1:50567/142‑1 (MQ=255)
                                                                                                                     gTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGc                                    >  2:50567/1‑142 (MQ=255)
                                                                                                                       cAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGcatc                           <  1:149650/149‑1 (MQ=255)
                                                                                                                          gTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCatcatc                        >  1:247732/1‑149 (MQ=255)
                                                                                                                          gTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCatcatc                        >  1:296914/1‑149 (MQ=255)
                                                                                                                          gTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCatcatc                        >  1:242082/1‑149 (MQ=255)
                                                                                                                               cGTATGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATc                                                                                                  <  4:351394/70‑1 (MQ=255)
                                                                                                                               cGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCATCATCTtgtg                   >  2:135829/1‑149 (MQ=255)
                                                                                                                                gTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCATCATCTTGTgg                  <  1:133819/149‑1 (MQ=255)
                                                                                                                                 tAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGt                                                                                 <  1:145357/85‑1 (MQ=255)
                                                                                                                                 tAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCATCATCTTGTgg                  >  4:71020/1‑148 (MQ=255)
                                                                                                                                     cATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCATCATCTTGTGGGttt              <  3:77836/148‑1 (MQ=255)
                                                                                                                                       tGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCATCATCTTGTGGGTTTcgc           <  1:32920/149‑1 (MQ=255)
                                                                                                                                       tGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCATCATCTTGTGGGTTTcg            >  3:54561/1‑148 (MQ=255)
                                                                                                                                        gATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCATCATCTTGTGGGTTTcgcg          >  1:317533/1‑149 (MQ=255)
                                                                                                                                          tgcgACGCTTACGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGcatcat                         >  4:247589/4‑132 (MQ=255)
                                                                                                                                          tgcgACGCTTACGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGcatcat                         <  3:247589/129‑1 (MQ=255)
                                                                                                                                          tAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCATCATCTTGTGGGTTTCGCGa         >  4:190285/1‑148 (MQ=255)
                                                                                                                                               cgcgTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCATCATCTTGTGGGTTTCGCGACCGAt    >  3:340008/1‑148 (MQ=255)
                                                                                                                                               cgcgTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCATCATCTTGTGGGTTTCGCGACCGATa   <  2:352630/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAACGTACTTCTGTTGGCGTAGTCTGGGTTATTGCGCACGTAGGTTTAAAACGTAAATTGTAACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGGTCAGGTTCCCGTAGGCATGATAAGACGCGTAAGCGTCGCATCAGGCAATGAATACCCAATGCGACCAGCTTCTTATATCAGAACAGCCCCAACGGTTTATCCGAGTAGCTCACCAGCAGGCACTTGGTTTGCTGGTAATGCTCCAGCATCATCTTGTGGGTTTCGCGACCGATAC  >  NC_000913/3754778‑3755070

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: