Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,552,838 Δ1 bp 6.2% intergenic (+17/‑64) 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,8380T. 6.2% 196.4 / 7.5 65intergenic (+17/‑64)uxuA/uxuBmannonate hydrolase/D‑mannonate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base T (32/28);  new base . (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

CACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATTGTGCA  >  NC_000913/4552690‑4552981
                                                                                                                                                    |                                                                                                                                                
cACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCt                                                                                                                                                  >  1:306458/1‑149 (MQ=255)
  cGGTCATCAGATGCTGGACGACCTGAAGAAGAAACCCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTgc                                                                                                                                                <  1:393027/149‑1 (MQ=255)
           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                                                                                                                  >  1:502911/1‑149 (MQ=255)
                               agaaAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCg                                                                                                                  >  2:781209/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                                                               <  2:502911/149‑1 (MQ=255)
                                                                                  aaGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACa                                                               <  1:781209/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)
                                                                                                                                                    |                                                                                                                                                
CACGGTCATCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATTGTGCA  >  NC_000913/4552690‑4552981

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