Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,552,842:1 +A 6.2% intergenic (+21/‑60) uxuA → / → uxuB mannonate hydrolase/D‑mannonate oxidoreductase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,552,8421.A 6.2% 256.5 / 11.4 64intergenic (+21/‑60)uxuA/uxuBmannonate hydrolase/D‑mannonate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base . (32/28);  new base A (2/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

GATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATTGTGCA  >  NC_000913/4552701‑4552981
                                                                                                                                              |                                                                                                                                           
gATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGc                                                                                                                                      >  2:613711/1‑149 (MQ=255)
gATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGc                                                                                                                                      <  2:639765/149‑1 (MQ=255)
   gCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGATGCGC‑GCCGTGCCgg                                                                                                                                   >  2:638103/1‑149 (MQ=255)
      gggcgacCCGAAGAATAAAAGCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTc                                                                                                                                <  1:2495/146‑1 (MQ=255)
      ggACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCACGGC‑GCGCAGCCATGCCGGtgg                                                                                                                                <  1:452613/149‑3 (MQ=255)
      gaccgacCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGCCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTc                                                                                                                                <  2:517935/146‑1 (MQ=255)
       gacgacCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATGTCCACCGGCCTGGCTCCGC‑GCCGTGCCGGGTcc                                                                                                                               >  2:55577/1‑149 (MQ=255)
           acCTGAAGAAGAAATCCAACCCAGGTTACTCCGCAATTGGTCGTCGGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCcttct                                                                                                                           <  2:396356/149‑1 (MQ=255)
           acCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAAGTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCACGGC‑GCGCAGCCATGCCGGtggatatt                                                                                                                           <  2:200267/149‑8 (MQ=255)
                   aagaaAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTAAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGcc                                                                                                                   >  2:430399/1‑149 (MQ=255)
                    agaaAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCg                                                                                                                  >  2:781209/1‑149 (MQ=255)
                    agaaAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCg                                                                                                                  >  1:502911/1‑149 (MQ=255)
                    acaaaaCCATCCCAGGTGACTCCGCAAGTGGTCGTCGGAAAGGCCTGGCCGAAGTCCGCGGTGTCTAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGCCATGGCTGCGC‑TCCGTGCCGTTTCCTTCTGCCTTGCCg                                                                                                                  <  2:354427/147‑1 (MQ=255)
                       aaaCCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCa                                                                                                               <  1:560153/149‑1 (MQ=255)
                                cAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTATAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGaa                                                                                                      <  2:580539/149‑1 (MQ=255)
                                   gTTACTCCGCAATTGGCCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGAc                                                                                                   <  2:735245/149‑1 (MQ=255)
                                       ccccGCAATGGGTCGTCGGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGAt                                                                                               <  1:171186/147‑1 (MQ=255)
                                        tCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGAtt                                                                                              >  2:686143/1‑149 (MQ=255)
                                          cGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCt                                                                                            >  2:706568/1‑149 (MQ=255)
                                          cGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGAACCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCt                                                                                            >  1:25979/1‑149 (MQ=255)
                                            cAATTGGTCGTCTGAAAGGCCTAGCGTAAGTTCGCGGTGTCGATCTGGCGATCCAGCGCGCTTTCTTTAGCCGTTACTATCCACCGCCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTGCCTTGCCGTCACTCTTTGAAGACGGATTCTgg                                                                                          <  1:7401/149‑1 (MQ=255)
                                             aactGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGa                                                                                         <  1:72919/146‑1 (MQ=255)
                                             aaTTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGa                                                                                         <  2:629873/149‑1 (MQ=255)
                                                tGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGACGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGtt                                                                                      >  2:472271/1‑149 (MQ=255)
                                                   tcgtcTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTGGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTAc                                                                                   <  1:558108/149‑1 (MQ=255)
                                                   tcgtcTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGGGTTTAc                                                                                   >  1:276267/1‑149 (MQ=255)
                                                         gAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATggc                                                                             >  2:314245/1‑147 (MQ=255)
                                                                cTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTAtt                                                                      <  1:618812/149‑1 (MQ=255)
                                                                 tGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTAttg                                                                     >  2:60016/1‑149 (MQ=255)
                                                                    ccGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTAttgttg                                                                  >  2:81683/1‑149 (MQ=255)
                                                                     cGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGCTGCCTACTAttgttgg                                                                 >  1:444552/1‑148 (MQ=255)
                                                                     cGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGa                                                                 >  1:663496/1‑149 (MQ=255)
                                                                       aaGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACa                                                               <  1:781209/149‑1 (MQ=255)
                                                                       aaGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACa                                                               <  2:502911/149‑1 (MQ=255)
                                                                          ttCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCa                                                            <  1:602692/149‑1 (MQ=255)
                                                                                  gtCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGACGACTACTATTGTTGACAGCAATCTGCCg                                                    <  1:617540/149‑1 (MQ=255)
                                                                                              aTCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGcccgcccg                                        >  2:88827/1‑149 (MQ=255)
                                                                                              aTCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGAACGCAATCTCCCGGTTgcccgcccc                                        >  2:693384/1‑148 (MQ=255)
                                                                                              aTCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGGTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGaccgcccg                                        >  1:93162/1‑149 (MQ=255)
                                                                                                ccAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTc                                      <  1:643162/149‑1 (MQ=255)
                                                                                                 cAGCGGGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCa                                     <  2:739820/149‑1 (MQ=255)
                                                                                                 cAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCa                                     <  1:28635/149‑1 (MQ=255)
                                                                                                 cAGCGCGCTTTCTTGAGCCGTTAACATCCACCGGCATGGCCGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGCCGGCTTCTTGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCa                                     <  2:285121/149‑1 (MQ=255)
                                                                                                   gcgcgcTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGGAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGCCATg                                   >  1:522736/1‑149 (MQ=255)
                                                                                                   gcgcgcTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCATTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATg                                   >  1:522835/1‑149 (MQ=255)
                                                                                                   gcgcgcTGTCTTTGGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATg                                   <  2:236741/149‑1 (MQ=255)
                                                                                                   gcgcgcATTCTTTAGCCGTTAATACCCACCGGCACGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACTTATCCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCAGTCATg                                   <  1:23372/149‑1 (MQ=255)
                                                                                                    cgcgcTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATgg                                  <  1:367357/149‑1 (MQ=255)
                                                                                                      cgcTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGa                                <  1:81683/149‑1 (MQ=255)
                                                                                                        ctttctttAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATc                              >  2:346060/1‑149 (MQ=255)
                                                                                                                  ccGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCt                    >  2:351187/1‑149 (MQ=255)
                                                                                                                    gTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTgg                  >  2:176421/1‑149 (MQ=255)
                                                                                                                    gTTAATATCCACCGGCATGGCTGAGC‑GCCGTGCCGGTTACTTCTTCATTGCCGTCACTCTTTGAAGACGGATGCTGGAGTTTACGATGACTACTATTGTTGACAGCGATCTACCGGTTGCCAGACAAACAGGGGATCATTCTCGTCggg                  >  2:585348/1‑149 (MQ=255)
                                                                                                                     ttAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGa                 >  1:670109/1‑149 (MQ=255)
                                                                                                                             cACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACTGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGGTCATTCTCGGCTGGAATCCCGCa         >  1:520572/1‑149 (MQ=255)
                                                                                                                             cACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCTTTGGGAATcgcggc         >  2:415907/1‑147 (MQ=255)
                                                                                                                             cACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCa         >  1:655774/1‑149 (MQ=255)
                                                                                                                             cACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCa         >  1:45717/1‑149 (MQ=255)
                                                                                                                             cACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCa         >  2:156629/1‑149 (MQ=255)
                                                                                                                             cACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACCACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCGTGGGATCATTCTCGTCTGGAATCACGCa         <  1:472271/149‑1 (MQ=255)
                                                                                                                              aCCGGCACGGC‑GCGCAGCCATGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCAt        >  1:348968/1‑149 (MQ=255)
                                                                                                                              aCCGGCACGGC‑GCGCAGCCATGCCGCTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGCACTCACGCAt        >  1:391450/1‑149 (MQ=255)
                                                                                                                                 ggCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATtgt     <  1:60016/149‑1 (MQ=255)
                                                                                                                                    aTGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATTGTGCa  >  1:237702/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                           
GATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATTGTGCA  >  NC_000913/4552701‑4552981

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