Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,873,370 A→T 6.7% C81* (TGT→TGA 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,3700AT 6.7% 228.8 / 4.1 80C81* (TGT→TGAdgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (27/45);  new base T (3/3);  total (31/49)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTC  >  NC_000913/3873223‑3873513
                                                                                                                                                   |                                                                                                                                               
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:24853/149‑1 (MQ=255)
                         gcgcCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGAc                                                                                                                       <  5:201466/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                                                                                                <  5:218278/149‑1 (MQ=255)
                                                ggCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAgtaccgtac                                                                                                <  1:340399/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                       >  1:41798/1‑149 (MQ=255)
                                                                                                                         tGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAAtt                       <  3:65454/149‑1 (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                                                                   >  6:150768/1‑92 (MQ=255)
                                                                                                                                      tAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTAc                                                                   <  5:150768/92‑1 (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)
                                                                                                                                               gCTGACAGCCAATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTc  >  3:102129/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTC  >  NC_000913/3873223‑3873513

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: