Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,068,014 A→T 5.6% S421R (AGT→AGA)  yihQ ← putative sulpholipid alpha‑glucosidase; alpha‑glucosyl fluoride glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,068,0140AT 5.6% 224.5 / 3.3 72S421R (AGT→AGA) yihQputative sulpholipid alpha‑glucosidase; alpha‑glucosyl fluoride glucosidase
Reads supporting (aligned to +/- strand):  ref base A (30/38);  new base T (2/2);  total (32/40)
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.21e-01

GCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCG  >  NC_000913/4067868‑4068161
                                                                                                                                                  |                                                                                                                                                   
gcccACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACt                                                                                                                                                   <  5:351040/149‑1 (MQ=255)
    acaTCATGGTGGAGTATTTCTGGCTACCGGTAGCACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGATTCTTCACGGGCGTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACg                                                                                                                                               >  5:104498/1‑149 (MQ=255)
     catcatGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGGTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGc                                                                                                                                              <  6:171233/149‑1 (MQ=255)
      atcatGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGcc                                                                                                                                             <  1:3300/149‑1 (MQ=255)
         atGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGt                                                                                                                                           <  7:14197/148‑1 (MQ=255)
            gtggAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTAt                                                                                                                                        <  4:200021/148‑1 (MQ=255)
             tggAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATg                                                                                                                                       >  7:250942/1‑148 (MQ=255)
             tggAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGc                                                                                                                                      <  6:104093/149‑1 (MQ=255)
               gAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGGCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCaa                                                                                                                                    <  3:36336/149‑1 (MQ=255)
                 gTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGt                                                                                                                                  <  7:329919/149‑1 (MQ=255)
                          gCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACgtgt                                                                                                                            <  4:191509/146‑1 (MQ=255)
                              ccGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCCGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTggg                                                                                                                     >  1:301377/1‑149 (MQ=255)
                              ccGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTg                                                                                                                       <  6:277647/147‑1 (MQ=255)
                              ccGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGAACGTGTCGGTggg                                                                                                                     >  3:51945/1‑149 (MQ=255)
                               cGGTAGAACCGGCGCGCATAAAGTATAGTATCTCGCCGAGGTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCTCAGCGCATGCCAGGCGTTATGCATAATTTGGGCACTGACGCGGTTATGCAAGTACGTGTCGGTGGGc                                                                                                                    <  6:104498/149‑1 (MQ=255)
                                 gTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTggg                                                                                                                     >  2:221687/1‑146 (MQ=255)
                                 gTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAg                                                                                                                  <  3:295013/149‑1 (MQ=255)
                                    gAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATa                                                                                                               <  8:250942/149‑1 (MQ=255)
                                          gcgcgcATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCACACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGcc                                                                                                         <  8:336161/149‑1 (MQ=255)
                                               cATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGt                                                                                                    >  8:76853/1‑149 (MQ=255)
                                                 tAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGt                                                                                                    >  5:101514/1‑147 (MQ=255)
                                                     gAAAAGGATCTCGCCGAGCTTGCCCGTATCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCa                                                                                              >  2:288268/1‑149 (MQ=255)
                                                         aGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCa                                                                                          <  3:291226/149‑1 (MQ=255)
                                                                    gAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGc                                                                               >  3:103556/1‑149 (MQ=255)
                                                                     aGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGcc                                                                              <  7:76853/149‑1 (MQ=255)
                                                                      gCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGcc                                                                              >  7:146892/1‑148 (MQ=255)
                                                                            ccGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTAc                                                                                                                                <  2:161908/92‑1 (MQ=255)
                                                                            ccGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTAc                                                                                                                                >  1:161908/1‑92 (MQ=255)
                                                                                   ttcAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGt                                                                >  4:103250/1‑149 (MQ=255)
                                                                                      aaGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGAACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTc                                                              <  4:53262/148‑1 (MQ=255)
                                                                                        gggCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCttt                                                           >  2:383423/1‑149 (MQ=255)
                                                                                         ggCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCtttt                                                          <  3:139523/149‑1 (MQ=255)
                                                                                            ttCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCTCTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGa                                                       >  6:319988/1‑149 (MQ=255)
                                                                                              cGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTTTGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGAt                                                      <  1:288268/148‑1 (MQ=255)
                                                                                               gTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATc                                                     <  6:101514/148‑1 (MQ=255)
                                                                                                  gTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGACAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACtt                                                 >  3:277191/1‑149 (MQ=255)
                                                                                                   ttACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTc                                                >  6:298846/1‑149 (MQ=255)
                                                                                                     acacTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTcc                                               <  1:383423/148‑1 (MQ=255)
                                                                                                       acTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCGATcc                                                                                           <  1:362459/102‑1 (MQ=255)
                                                                                                       acTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATcc                                                                                           >  2:362459/1‑102 (MQ=255)
                                                                                                         ttCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTg                                                                                                                       >  4:94492/1‑72 (MQ=255)
                                                                                                         ttCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTg                                                                                                                       <  3:94492/72‑1 (MQ=255)
                                                                                                         ttCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTATGATCACTTCCTTGaa                                          >  5:205526/1‑149 (MQ=255)
                                                                                                             cccACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAg                                      >  2:28572/1‑149 (MQ=255)
                                                                                                             cccACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAg                                      <  1:221687/149‑1 (MQ=255)
                                                                                                             cccACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGACAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAg                                      <  4:277191/149‑1 (MQ=255)
                                                                                                              ccACAGCGCAGGCCAGGCGTAATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAgg                                     >  2:79089/1‑149 (MQ=255)
                                                                                                               cacaGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATATTCTTTTTGATCACTTCCTTg                                            >  1:148912/1‑141 (MQ=255)
                                                                                                               cacaGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATATTCTTTTTGATCACTTCCTTg                                            <  2:148912/141‑1 (MQ=255)
                                                                                                                acaGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACt                                                  >  1:203019/1‑134 (MQ=255)
                                                                                                                acaGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACt                                                  <  2:203019/134‑1 (MQ=255)
                                                                                                                 caGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGAACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGt                                  <  4:51945/149‑1 (MQ=255)
                                                                                                                         gccaggcGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCt                          <  1:28572/149‑1 (MQ=255)
                                                                                                                           caggcGTTATGCATAATTTCGGCACTGACGCCGTTATGCAATTACGTGTCGGTGGGCAGATACTCGCCGAAGACAGCAATACAGCCGCCGAAGCCGAGTTCAATCATGTTCATTTTCATCACTTCCTTGACCCCGGCCTAGGCtaatg                         >  3:82090/1‑144 (MQ=255)
                                                                                                                                 ttATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGt                  <  5:298846/149‑1 (MQ=255)
                                                                                                                                  ttttcATAATGTCGGCACTGACGCCGTTATTCAAGTACGTGTGGGTGGGCAGATACTCGCCGAAGTCAGCCATCCTGCCGGCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGt                  <  4:82090/144‑1 (MQ=255)
                                                                                                                                    tgaaTAATTTCGGCTCTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCtttt                                                          <  1:170888/103‑1 (MQ=255)
                                                                                                                                    tGCATAATTTCGGCTCTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCtttt                                                          >  2:170888/1‑106 (MQ=255)
                                                                                                                                    tGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCt                          >  4:283159/1‑138 (MQ=255)
                                                                                                                                    tGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCt                          <  3:283159/138‑1 (MQ=255)
                                                                                                                                    tGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGGCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTga                <  2:301377/148‑1 (MQ=255)
                                                                                                                                     gCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGt                  >  3:98416/1‑145 (MQ=255)
                                                                                                                                      cATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTa                    <  8:232593/142‑1 (MQ=255)
                                                                                                                                      cATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTa                    >  7:232593/1‑142 (MQ=255)
                                                                                                                                       aTAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATc                                                                                            <  1:10176/69‑1 (MQ=255)
                                                                                                                                       aTAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATc                                                                                            >  2:10176/1‑69 (MQ=255)
                                                                                                                                       aTAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTg                                            <  7:16409/117‑1 (MQ=255)
                                                                                                                                       aTAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTg                                            >  8:16409/1‑117 (MQ=255)
                                                                                                                                          aTTTCGGCTCTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGAc         <  5:319988/149‑1 (MQ=255)
                                                                                                                                           tttCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGAc         >  4:334308/1‑148 (MQ=255)
                                                                                                                                             tCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCTACaa       <  1:79089/148‑1 (MQ=255)
                                                                                                                                                 cACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAAcgccg  >  1:169615/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGGCCAGGCGTTATGCATAATTTCGGCACTGACGCCGTTATGCAAGTACGTGTCGGTGGGCAGATACTCGCCGAAGTCAGCCATCCAGCCGCCGCAGCCGAGTTCAATCATGTTCTTTTTGATCACTTCCTTGAACCAGGCGTAGGCTTCTGGATTAGTGAGATCGACAACGCCG  >  NC_000913/4067868‑4068161

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