Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,873,365 A→G 6.0% L83P (CTC→CCC)  dgoA ← 2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,873,3650AG 6.0% 251.1 / 5.5 84L83P (CTC→CCC) dgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (28/49);  new base G (3/2);  total (33/51)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

GCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGT  >  NC_000913/3873218‑3873512
                                                                                                                                                   |                                                                                                                                                   
gCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAg                                                                                                                                                    <  3:58016/149‑1 (MQ=255)
     tGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGa                                                                                                                                                >  1:241479/1‑148 (MQ=255)
     ccgatctaaTATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGa                                                                                                                                                <  7:145013/141‑1 (MQ=255)
         ggAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGa                                                                                                                                                >  2:68507/1‑144 (MQ=255)
         ggAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGa                                                                                                                                                <  1:68507/144‑1 (MQ=255)
             atatTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGa                                                                                                                                                >  8:145013/1‑140 (MQ=255)
                ttttCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTg                                                                                                                                    <  7:48101/149‑1 (MQ=255)
                  ttCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGc                                                                                                                                  >  6:24853/1‑149 (MQ=255)
                   tCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGc                                                                                                                                  <  5:311426/148‑1 (MQ=255)
                      gcgcCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAgc                                                                                                                              <  4:156345/149‑1 (MQ=255)
                           tGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATc                                                                                                                         >  2:340399/1‑149 (MQ=255)
                              gcgcCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGAc                                                                                                                      <  5:201466/149‑1 (MQ=255)
                              gcgcCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGAc                                                                                                                      <  5:24853/149‑1 (MQ=255)
                                 cccGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTg                                                                                                                   <  7:364965/149‑1 (MQ=255)
                                    gCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGtt                                                                                                                 >  3:185680/1‑148 (MQ=255)
                                    gCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGtt                                                                                                                 >  2:294735/1‑148 (MQ=255)
                                    gCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTc                                                                                                                <  4:1183/149‑1 (MQ=255)
                                        cGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAgg                                                                                                             <  3:147257/148‑1 (MQ=255)
                                         gAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGt                                                                                                            <  8:248065/148‑1 (MQ=255)
                                            cgcgGGAAAGGCGTAGGTTGGAGTCGGTGAGCGGGTGCAGTCGTTGATTCTGTAGCCCACGGTACGGCGGATCACTTCACAATGGATCTTGGGCTGTAGGATGAGGTGTCAGTTCATGCTGTAGAGCGCTTCGACCTGTTCAGGTTTCa                                                                                                        <  2:221240/149‑1 (MQ=255)
                                              cgGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAgt                                                                                                      <  1:293523/149‑1 (MQ=255)
                                               ggTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAgta                                                                                                     <  1:348695/149‑1 (MQ=255)
                                                gTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAgtac                                                                                                    <  3:145411/149‑1 (MQ=255)
                                                 tAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAgtacc                                                                                                   <  6:243485/149‑1 (MQ=255)
                                                  aaaGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAgtaccg                                                                                                  >  7:167205/1‑149 (MQ=255)
                                                     ggCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGCGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAgtaccgtac                                                                                               >  2:131691/1‑149 (MQ=255)
                                                     ggCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAgtaccgtac                                                                                               <  1:340399/149‑1 (MQ=255)
                                                     ggCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAgtaccgtac                                                                                               <  5:218278/149‑1 (MQ=255)
                                                      gCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAgtaccgtacc                                                                                              <  2:241479/149‑1 (MQ=255)
                                                       cTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCt                                                                                             <  8:263406/149‑1 (MQ=255)
                                                           ggTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTgcgc                                                                                         <  2:127547/149‑1 (MQ=255)
                                                            gTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTgcgc                                                                                         >  8:180762/1‑148 (MQ=255)
                                                            gTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCAATCCTGGCGAGCGCATCGACCTGTGCAGGTTTCAGTACCGTACCTgcgc                                                                                         >  4:334027/1‑148 (MQ=255)
                                                            gTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACAGCACGGCGCATCACTTCACTATGGATATTGGGCGTAACGATGCGCTGACAGCCCCTCCTGGCGAGCGCAACGACCTGAGCAGGGTTCCGTACCGTACCTgcgc                                                                                         >  2:264454/1‑148 (MQ=255)
                                                               gCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAAt                                                                                     <  5:340901/149‑1 (MQ=255)
                                                                  gTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCaa                                                                                  <  7:180762/149‑1 (MQ=255)
                                                                   tCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAAc                                                                                 <  5:214278/149‑1 (MQ=255)
                                                                       gcAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAAc                                                                                 <  2:23710/145‑1 (MQ=255)
                                                                       gcAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAAc                                                                                 >  1:23710/1‑145 (MQ=255)
                                                                        cAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCAACCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCtt                                                                            >  6:43233/1‑149 (MQ=255)
                                                                                gCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGt                                                                    <  7:145037/149‑1 (MQ=255)
                                                                                          aTGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAAc                                                          >  6:281871/1‑149 (MQ=255)
                                                                                             ccGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGAt                                                       >  8:291966/1‑149 (MQ=255)
                                                                                             ccGTAGCCCACCGCACGGCGGATAACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCAGGCGAGCGCATCGACCGGTTCAGGTTTCAGTACCGTACCAGCGCCAAACAACCCCTTGTCGCCGTACGCAAAAACGAt                                                       >  6:348790/1‑149 (MQ=255)
                                                                                              cGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATg                                                      <  4:174876/149‑1 (MQ=255)
                                                                                                   cccACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCgg                                                  >  8:346732/1‑148 (MQ=255)
                                                                                                          cACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCtt                                          <  1:44910/149‑1 (MQ=255)
                                                                                                          cACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCtt                                          <  8:167205/149‑1 (MQ=255)
                                                                                                             ggcggATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGc                                       <  6:170827/149‑1 (MQ=255)
                                                                                                             ggcggAACACTTCACTATCGATCTTGGGCGTCACGAGGAGCCGCCACCCCATCCTGGCGAGCGCATCGACCTGGTCAGGTGTCAGTACCGTACCTGCGCCCATCACCGCGTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGc                                       <  6:9899/149‑1 (MQ=255)
                                                                                                                ggATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTcc                                    >  2:25587/1‑149 (MQ=255)
                                                                                                                  aTCACTTCACTCTGGATATTGGGCGTCACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGGGCCAATCAACGCCTTGTCGCCGTACGCATCCACGATGGCGGGAATGCTTTGCTccct                                  <  3:334027/149‑2 (MQ=255)
                                                                                                                  aTCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCa                                  <  1:110856/149‑1 (MQ=255)
                                                                                                                  aTCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCa                                  <  5:348790/149‑1 (MQ=255)
                                                                                                                      cttcactATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATtgt                              <  1:294735/149‑1 (MQ=255)
                                                                                                                       ttcactATGGATATTGGGCGTAACGATGAGCTGCCAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATtgtg                             >  4:294214/1‑149 (MQ=255)
                                                                                                                         cactATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATtgtg                             <  4:185680/147‑1 (MQ=255)
                                                                                                                         cactATGGATATTGGGCGTAACGAGGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTgg                            <  4:173970/148‑1 (MQ=255)
                                                                                                                          actATGGATATTGGGCTTAACGATCAGCTGACATCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGAACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTggg                           <  6:125144/148‑1 (MQ=255)
                                                                                                                            tATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCa                                  <  4:1364/139‑1 (MQ=255)
                                                                                                                            tATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCa                                  >  3:1364/1‑139 (MQ=255)
                                                                                                                             aTGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAAt                       >  4:65454/1‑149 (MQ=255)
                                                                                                                              tGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAAtt                      <  3:65454/149‑1 (MQ=255)
                                                                                                                              tGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAAtt                      >  1:41798/1‑149 (MQ=255)
                                                                                                                                gATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCa                    <  7:291966/149‑1 (MQ=255)
                                                                                                                                   atTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCg                 >  8:279893/1‑149 (MQ=255)
                                                                                                                                    ttGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTCCCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCgg                <  3:227117/149‑1 (MQ=255)
                                                                                                                                    ttGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCgg                <  3:294214/149‑1 (MQ=255)
                                                                                                                                    ttGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCgg                >  7:51077/1‑149 (MQ=255)
                                                                                                                                    ttGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCg                 <  5:281871/148‑1 (MQ=255)
                                                                                                                                          gTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGAttt           >  2:158971/1‑148 (MQ=255)
                                                                                                                                          gTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTc          <  1:25587/149‑1 (MQ=255)
                                                                                                                                           tAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTAc                                                                  <  5:150768/92‑1 (MQ=255)
                                                                                                                                           tAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTAc                                                                  >  6:150768/1‑92 (MQ=255)
                                                                                                                                            caggATGGGCTGTCAGCTCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCaa        >  2:235893/4‑149 (MQ=255)
                                                                                                                                            aaCGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCa                                                                                                        <  1:352164/53‑1 (MQ=255)
                                                                                                                                            aaCGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCa                                                                                                        >  2:352164/1‑53 (MQ=255)
                                                                                                                                             aggATGGGCTGTCAGCTCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGAtt            <  1:93941/142‑1 (MQ=255)
                                                                                                                                             aggATGGGCTGTCAGCTCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGAtt            <  7:314106/142‑1 (MQ=255)
                                                                                                                                             aggATGGGCTGTCAGCTCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGAtt            >  2:93941/3‑144 (MQ=255)
                                                                                                                                             aggATGGGCTGTCAGCTCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGAtt            >  8:314106/3‑144 (MQ=255)
                                                                                                                                               gATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCAtt                                >  2:23785/1‑122 (MQ=255)
                                                                                                                                               gATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCAtt                                <  1:23785/122‑1 (MQ=255)
                                                                                                                                                aTGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACcgc    >  5:96683/1‑149 (MQ=255)
                                                                                                                                                   aGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGt  <  8:96811/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGT  >  NC_000913/3873218‑3873512

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