Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 696,413 A→G 5.5% intergenic (+137/+17) ubiF → / ← glnX 2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol oxygenase/tRNA‑Gln

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913696,4130AG 5.5% 227.9 / 2.6 72intergenic (+137/+17)ubiF/glnX2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol oxygenase/tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base A (28/40);  new base G (2/2);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

GTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGA  >  NC_000913/696273‑696555
                                                                                                                                            |                                                                                                                                              
gTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCt                                                                                                                                         <  7:221124/148‑1 (MQ=255)
 tAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACTCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGAAATTGAATGAACGCAGAAAAGCAAAAAGCTCGCAGAAGAGAGCttt                                                                                                                                       >  8:121533/1‑149 (MQ=255)
 tAGCCTTACAACATTGCCGGATGCGTGCCAACCGGAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGACCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCttt                                                                                                                                       <  1:65692/149‑1 (MQ=255)
   gCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCttttt                                                                                                                                     <  8:83802/149‑1 (MQ=255)
    ccTTACAACATTGCCGGATGCGTGCCAACCGTCGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCtttttt                                                                                                                                    <  6:290121/149‑1 (MQ=255)
    ccTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGCCGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCtttttt                                                                                                                                    <  7:210782/149‑1 (MQ=255)
    ccTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCTTCGGCGAGCtttttg                                                                                                                                    <  1:190374/149‑2 (MQ=255)
    ccTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCttttt                                                                                                                                     >  4:239288/1‑148 (MQ=255)
    ccTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCttttt                                                                                                                                     >  1:326981/1‑148 (MQ=255)
        acaacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGAAAAGACGGGTCAACATCGCATTCGACATTGAATGAACGCAGAAAACCAAAAAGCTCGCCGAAGCGAGCTTTTTTaa                                                                                                                                  <  1:111520/147‑1 (MQ=255)
          aacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAAtgtg                                                                                                                              >  7:26953/1‑149 (MQ=255)
          aacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAAtgtg                                                                                                                              >  2:247772/1‑149 (MQ=255)
                 ccGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAAt                                                                                                                                 >  1:139736/1‑139 (MQ=255)
                 ccGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAAt                                                                                                                                 <  2:139736/139‑1 (MQ=255)
                 ccGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTgggg                                                                                                                       >  2:17467/1‑149 (MQ=255)
                  cGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGt                                                                                                                      <  5:211702/149‑1 (MQ=255)
                         gTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGa                                                                                                               >  1:88667/1‑149 (MQ=255)
                         gTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGa                                                                                                               >  1:238779/1‑149 (MQ=255)
                            ccAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGAtt                                                                                                             <  4:170331/148‑1 (MQ=255)
                               aCCGTAGGTCGGATAAGACCCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGGGTAAACATCGCATTAGACATTGAATGAACGCAGAAAAGCAAAAAGCTAGCCGAAGCGAGCTTTTTTAATGTGGCTTGGGTACGAGGATTCGaa                                                                                                         >  2:4625/1‑149 (MQ=255)
                                   tAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTc                                                                                                     <  1:306595/149‑1 (MQ=255)
                                   tAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTc                                                                                                     <  8:261955/149‑1 (MQ=255)
                                   tAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTc                                                                                                     <  8:26953/149‑1 (MQ=255)
                                     ggTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCg                                                                                                    >  8:324850/1‑148 (MQ=255)
                                      gTCGGATAAGACGCGGTAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCgg                                                                                                   >  6:101636/1‑148 (MQ=255)
                                       tCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGaa                                                                                                 <  1:17467/149‑1 (MQ=255)
                                       tCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGaa                                                                                                 <  7:121533/149‑1 (MQ=255)
                                         ggATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGCTTCGCACCTCGGAATg                                                                                               >  1:228007/1‑149 (MQ=255)
                                          gATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGc                                                                                              <  5:183396/149‑1 (MQ=255)
                                           aTAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGcc                                                                                             <  3:111215/149‑1 (MQ=255)
                                            tAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCg                                                                                            <  1:247772/149‑1 (MQ=255)
                                               gACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAAt                                                                                                                                 >  4:131200/1‑109 (MQ=255)
                                                 cgcgGCAGCGCCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATc                                                                                       <  2:212660/149‑1 (MQ=255)
                                                  gcgGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCa                                                                                      <  1:7548/149‑1 (MQ=255)
                                                    ggCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGa                                                                                    <  2:88667/149‑1 (MQ=255)
                                                       aGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATc                                                                                 <  2:326981/149‑1 (MQ=255)
                                                         cGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCg                                                                               <  3:239288/149‑1 (MQ=255)
                                                           tCGCATCCGACATTGAAGGATAAGACGTGTCAGCATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCAGGAGTCTGCATCCGGt                                                                             >  3:74952/1‑149 (MQ=39)
                                                           tCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGCTTCGAACCTCGGAATGCCGGAATCAGACTCCGGt                                                                             >  6:287157/1‑149 (MQ=255)
                                                             gCATCCGACATTGAAGGATAAGCCGTGTCCACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGc                                                                           >  1:44277/1‑149 (MQ=255)
                                                              cATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGcc                                                                          >  3:65528/1‑149 (MQ=255)
                                                                  cGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAc                                                                      <  1:61840/149‑1 (MQ=255)
                                                                   gACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAGTGCCGGAATCAGAATCAGGTGCCTTagc                                                                     >  4:12694/1‑147 (MQ=255)
                                                                   gACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAcc                                                                     >  5:214800/1‑149 (MQ=255)
                                                                      aTTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                                                  >  1:152965/1‑149 (MQ=255)
                                                                       ttGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCtt                                                                 >  7:299964/1‑149 (MQ=255)
                                                                         gAACGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTg                                                                <  2:44544/148‑1 (MQ=255)
                                                                          aaGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAACGCAAAACGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                              <  1:86930/149‑1 (MQ=255)
                                                                                   aCGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATAcccc                                                      >  5:2763/1‑148 (MQ=255)
                                                                                     gtgtCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAAc                                                   <  1:45115/149‑1 (MQ=255)
                                                                                        tCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGTTACGAGGATTCGAACCTCGGAATGCCGTAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTgg                                                <  6:137179/149‑1 (MQ=255)
                                                                                        tCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGCACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTgg                                                >  5:72568/1‑149 (MQ=255)
                                                                                          aaCATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTggg                                               <  6:2763/148‑1 (MQ=255)
                                                                                             aTCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTgg                                                >  1:68263/1‑144 (MQ=255)
                                                                                             aTCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTgg                                                <  2:68263/144‑1 (MQ=255)
                                                                                                gCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCg                                                                                                    <  4:77648/89‑1 (MQ=255)
                                                                                                gCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCg                                                                                                    >  3:77648/1‑89 (MQ=255)
                                                                                                     cGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACaa                                    <  4:65528/148‑1 (MQ=255)
                                                                                                             aatgaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGc                                                                           <  1:141326/101‑1 (MQ=255)
                                                                                                             aatgaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGc                                                                           >  2:141326/1‑101 (MQ=255)
                                                                                                                 aaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCtt                        <  6:72568/148‑1 (MQ=255)
                                                                                                                 aaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCtt                        <  7:324850/148‑1 (MQ=255)
                                                                                                                      aGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATAcccc                                                      >  6:252744/1‑113 (MQ=255)
                                                                                                                      aGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATAcccc                                                      <  5:252744/113‑1 (MQ=255)
                                                                                                                              aaaaaGCTCGCTTCGGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTg                                                                >  2:82961/1‑95 (MQ=17)
                                                                                                                              aaaaaGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTg          <  3:24214/149‑1 (MQ=255)
                                                                                                                               aaaaGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTgg         <  8:299964/149‑1 (MQ=255)
                                                                                                                                aaaGCTCGCTTCGGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCtt                                                                 >  3:226244/1‑92 (MQ=17)
                                                                                                                                aaaGCTCGCTTCGGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCtt                                                                 <  4:226244/92‑1 (MQ=17)
                                                                                                                                       gCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGTGGATTCGAACCTCGGAGTGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCTTCTTGTATAAATTGGCTGGGGTACGa  <  3:12694/148‑1 (MQ=255)
                                                                                                                                       gCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTAAAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGa  <  7:99962/148‑1 (MQ=255)
                                                                                                                                       gCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCAGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGa  <  5:287157/148‑1 (MQ=255)
                                                                                                                                           aaGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGcc                                                                          <  4:120664/72‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
GTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGA  >  NC_000913/696273‑696555

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: ‑