Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,595,022 G→A 5.1% A277V (GCC→GTC)  lgoR ← putative transcriptional activator for L‑galactonate catabolism

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,595,0220GA 5.1% 243.0 / 2.1 78A277V (GCC→GTC) lgoRputative transcriptional activator for L‑galactonate catabolism
Reads supporting (aligned to +/- strand):  ref base G (37/36);  new base A (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.43e-01

TATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGAAAGATGACTGAGATGATCTCAAGCGATTGATCAAAAA  >  NC_000913/4594879‑4595169
                                                                                                                                               |                                                                                                                                                   
tATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATTAGCATAACGCTTGTTCTCATTAATGGATAGAATCATTGATTGATTAGCTTAATTCAAATGGTTACGGAGCTCAAGGGTGGCag                                                                                                                                                  <  8:131884/147‑2 (MQ=255)
  tttATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAg                                                                                                                                              >  2:314594/1‑149 (MQ=255)
  tgtATCTTTTGATACAACCGTGCCCCCGCTAATAGTTTGCGCTTTTTAATCCGTCCTTACTGAGCATAACGCGTGTTCTCATTAATGGATCGACTCATTGATTGTTTAGCGTAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAg                                                                                                                                              <  2:218351/147‑1 (MQ=255)
    tATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAgg                                                                                                                                             >  7:142773/1‑148 (MQ=255)
          ttGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGGTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCa                                                                                                                                          <  1:69173/145‑1 (MQ=255)
           tGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCg                                                                                                                                     >  8:179176/1‑149 (MQ=255)
             aTAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCgg                                                                                                                                    >  2:90061/1‑148 (MQ=255)
              tAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGc                                                                                                                                      >  6:312618/1‑145 (MQ=255)
              tAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGc                                                                                                                                      <  5:312618/145‑1 (MQ=255)
              tAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGc                                                                                                                                   >  7:253952/1‑148 (MQ=255)
                      tGCACCCGCTAATAGTTTGAGGTTTTTAATCAGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTTTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCagag                                                                                                                          <  2:223804/149‑1 (MQ=255)
                         aCCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGGGGCATCAAGGTCACTGCGGCAAATCAGAgcgc                                                                                                                       <  6:223506/149‑1 (MQ=255)
                         aCCCGCTAATAGTTTGCGGGTTTTAATCCGTGCTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGGTTAGCTGAGTTCATATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAgcgc                                                                                                                       >  6:86858/1‑149 (MQ=255)
                           ccGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCt                                                                                                                      >  7:104496/1‑148 (MQ=255)
                              cTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTgag                                                                                                                   >  5:309379/1‑148 (MQ=255)
                                 aTAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATg                                                                                                               <  8:137078/149‑1 (MQ=255)
                                  tAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGAAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATgg                                                                                                              >  8:167472/1‑149 (MQ=255)
                                   aGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGt                                                                                                             <  1:239948/149‑1 (MQ=255)
                                    gTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGt                                                                                                             >  4:8967/1‑148 (MQ=255)
                                    gTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGt                                                                                                             >  4:87716/1‑148 (MQ=255)
                                    gTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTc                                                                                                            <  3:144009/149‑1 (MQ=255)
                                     tttGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCa                                                                                                           <  2:98774/149‑1 (MQ=255)
                                     tttGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCa                                                                                                           <  6:309379/149‑1 (MQ=255)
                                        gCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCata                                                                                                         <  1:314594/148‑1 (MQ=255)
                                            tttttAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCt                                                                                                     >  2:97461/1‑148 (MQ=255)
                                            tttttAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTc                                                                                                    <  8:104496/149‑1 (MQ=255)
                                            tttttAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTc                                                                                                    <  7:35657/149‑1 (MQ=255)
                                             ttttAATCCGTACTTAATGAGCATACCGCGTGATCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCa                                                                                                   >  4:309929/1‑149 (MQ=255)
                                                    ccGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACgg                                                                                            >  2:178997/1‑149 (MQ=255)
                                                      gTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGc                                                                                           >  2:302035/1‑148 (MQ=255)
                                                        aCTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAAt                                                                                        >  2:8617/1‑149 (MQ=255)
                                                          ttAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGa                                                                                      >  5:261138/1‑149 (MQ=255)
                                                               gAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGGCATATGCTCATCGACGGCAATGATATTg                                                                                 >  1:138007/1‑149 (MQ=255)
                                                                   aTAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGACATCAAGGTCACTCCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCt                                                                             >  5:39146/1‑149 (MQ=255)
                                                                    tAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCata                                                                                                         <  2:78476/120‑1 (MQ=255)
                                                                    tAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCata                                                                                                         >  1:78476/1‑120 (MQ=255)
                                                                      aCGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTTCGGCAAATAAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGAGCTGtt                                                                          >  7:174728/1‑149 (MQ=255)
                                                                      aCGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGACATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGtt                                                                          >  5:94946/1‑149 (MQ=255)
                                                                        gcgTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGttt                                                                         >  7:274405/1‑148 (MQ=255)
                                                                           tgtTCTCTTTAATGGATCGAATGATTGATTGTTTAGCTGAATTCAATTGGTTACGGATCGCAAGGGTGTCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCGCATCGACGGCAATGATATTGCGCTGTTTAAGa                                                                     <  5:196259/149‑1 (MQ=255)
                                                                           tgtTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGATCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGa                                                                     <  1:249466/149‑1 (MQ=255)
                                                                           tgtTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGa                                                                     <  8:253952/149‑1 (MQ=255)
                                                                           tgtTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGACATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGa                                                                     <  6:94946/149‑1 (MQ=255)
                                                                               ctcATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAAAGGTTACGGAGCGCAAGCGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATCGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTaaaaaac                                                                  >  6:196259/1‑143 (MQ=255)
                                                                                tcATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCt                                                                >  5:146015/1‑149 (MQ=255)
                                                                                 cATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGAGCTGAGAATGGTCATATGCTCATCGACGCGAATGATATTGCGCTGTTTAAGATCGCtt                                                               >  1:95336/1‑149 (MQ=255)
                                                                                  aTTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCGAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCt                                                                <  2:321481/147‑1 (MQ=255)
                                                                                  aTTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCtt                                                               <  6:261138/148‑1 (MQ=255)
                                                                                      aTGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTc                                                          >  7:26374/1‑149 (MQ=255)
                                                                                      aTGGATCGAATAATTGATTGTTTAGCTGAATTCAAATGGTTACAGAGCGCTAGGTTGGCATCAAGTTAACTGCGGCAAATCAGAGCGCTGCGAATGGTCATATGCTCCTCGACTGCAATGATATTGCGCTGTTTAAGATGGCTTTCGTc                                                          >  4:70291/1‑149 (MQ=255)
                                                                                          aTCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCAt                                                      <  6:245951/149‑1 (MQ=255)
                                                                                           tCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACTGAGCGCAAGGGTGGCATCTAGGTCACTGCGTCAAATCATTGCTCTGTGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCAtt                                                     <  3:309929/149‑1 (MQ=255)
                                                                                             gAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGa                                                   <  1:90061/149‑1 (MQ=255)
                                                                                             gAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGa                                                   <  2:162737/149‑1 (MQ=255)
                                                                                               aTCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGa                                                   >  2:136875/1‑147 (MQ=255)
                                                                                               aTCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATa                                                 <  1:178997/149‑1 (MQ=255)
                                                                                               aTCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGAGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATa                                                 >  7:178169/1‑149 (MQ=255)
                                                                                               aTCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGACATCAAGGTCACTCCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATa                                                 <  6:39146/149‑1 (MQ=255)
                                                                                                   ttgattgTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTg                                              >  7:313019/1‑148 (MQ=255)
                                                                                                     gagtgTTTAGCTGAATTCAAATGGTTACGGAGCCCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGATTAACATCGCTTTCGTCCCAGTGATAGGGaa                                            >  1:141206/4‑148 (MQ=255)
                                                                                                      attgTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGaaaa                                          <  6:146015/149‑1 (MQ=255)
                                                                                                      attgTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGaaaa                                          <  1:97461/149‑1 (MQ=255)
                                                                                                       ttgTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAAt                                         <  3:87716/149‑1 (MQ=255)
                                                                                                          tttAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGaa                                      <  1:8617/149‑1 (MQ=255)
                                                                                                            tAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGaaa                                     <  1:136875/148‑1 (MQ=255)
                                                                                                                    ttCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGAAAGATGACTga                            <  3:8967/149‑1 (MQ=255)
                                                                                                                    ttCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGAGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGAAAGATGACt                              >  5:162386/1‑147 (MQ=255)
                                                                                                                      cAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGAAAGATGACt                              >  5:281543/1‑145 (MQ=255)
                                                                                                                       aaaTGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTc                                                                                                            <  7:212898/66‑1 (MQ=255)
                                                                                                                       aaaTGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTc                                                                                                            >  8:212898/1‑66 (MQ=255)
                                                                                                                       aaaTGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGa                                             <  7:341760/129‑1 (MQ=255)
                                                                                                                       aaaTGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGa                                             >  8:341760/1‑129 (MQ=255)
                                                                                                                            gTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGAAAGATGACTGAGATGAtc                     <  8:142773/148‑1 (MQ=255)
                                                                                                                                cGGAGCGCAAGTGTGGCATCAATGTCACTGCGGCAAATAAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCTATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGAAAGATGACTGAGATGATCTCaa                 <  8:174728/148‑1 (MQ=255)
                                                                                                                                 ggAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGATTTCGTCCCATTGATAGTGAAAATGAAAGATGACTGAGATGATCTCTAGc               >  5:269477/1‑149 (MQ=255)
                                                                                                                                     cgcAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCTTCCCATTGATAGTGAAAATGAAAGATGACTGAGATGATCTCAAGCGAtt           >  2:324013/1‑149 (MQ=255)
                                                                                                                                       cAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGAAAGATGACTGAGATGATCTCAAGCGATTGa         <  5:86858/149‑1 (MQ=255)
                                                                                                                                       cAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGAAAGATGACTGAGATGATCTCAAGCGATTGa         <  4:204377/149‑1 (MQ=255)
                                                                                                                                              ggCATCAACGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGAAAGATGACTGAGATGATCTCAAGCGATTGATCaaaaa  <  2:151356/149‑1 (MQ=255)
                                                                                                                                               gCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGAAAGATGACTGAGATGATCTCAAGCGATTGATCaaaa   >  6:256634/1‑147 (MQ=255)
                                                                                                                                               gCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGAAAGATGACTGAGATGATCTCAAGCGATTGATCaaaa   >  1:166884/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATATTGCGCTGTTTAAGATCGCTTTCGTCCCATTGATAGTGAAAATGAAAGATGACTGAGATGATCTCAAGCGATTGATCAAAAA  >  NC_000913/4594879‑4595169

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