Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,835,815 G→A 6.9% A101T (GCT→ACT)  ynjA → carboxymuconolactone decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,835,8150GA 6.9% 188.5 / 4.0 58A101T (GCT→ACT) ynjAcarboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base G (32/22);  new base A (2/2);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

GGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGCTCAGGCGCTCACCGAACT  >  NC_000913/1835667‑1835963
                                                                                                                                                    |                                                                                                                                                    
ggAGCGCAAACGATCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTgg                                                                                                                                                      <  1:705629/149‑1 (MQ=255)
   gcgcAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTg                                                                                                                                                   <  2:687720/149‑1 (MQ=255)
       aaaCGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTg                                                                                                                                               >  4:393888/1‑149 (MQ=255)
       aaaCGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTg                                                                                                                                               >  1:208222/1‑149 (MQ=255)
       aaaCGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTg                                                                                                                                               >  2:25131/1‑149 (MQ=255)
           gCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTggcgg                                                                                                                                           >  1:62075/1‑149 (MQ=255)
               aCCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCaaa                                                                                                                                       >  1:108082/1‑149 (MQ=255)
               aCCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCaaa                                                                                                                                       <  4:6944/149‑1 (MQ=255)
               aCCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCaaa                                                                                                                                       >  1:535204/1‑149 (MQ=255)
                       aTCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGctc                                                                                                                               <  2:340191/149‑1 (MQ=255)
                        tCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGctct                                                                                                                              >  1:423386/1‑149 (MQ=255)
                        tCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGACTGACTGGCGGCAAAGCCCGctct                                                                                                                              >  4:375054/1‑149 (MQ=255)
                        tCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGACTGACTGGCGGCAAAGCCCGctct                                                                                                                              >  2:100192/1‑149 (MQ=255)
                          cGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCttt                                                                                                                            <  2:423386/149‑1 (MQ=255)
                            gtagtaCGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAg                                                                                                                          >  4:698007/1‑149 (MQ=255)
                              agtaCGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCg                                                                                                                        >  1:9139/1‑149 (MQ=255)
                                taCGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGAt                                                                                                                      <  4:620699/149‑1 (MQ=255)
                                taCGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGAt                                                                                                                      <  4:714945/149‑1 (MQ=255)
                                  cGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATga                                                                                                                    <  1:321848/149‑1 (MQ=255)
                                  cGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATga                                                                                                                    >  4:202757/1‑149 (MQ=255)
                                         ttGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCg                                                                                                                                  >  3:649291/1‑128 (MQ=255)
                                         ttGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCg                                                                                                                                  <  4:649291/128‑1 (MQ=255)
                                         ttGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACgg                                                                                                             >  2:362541/1‑149 (MQ=255)
                                          tGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGc                                                                                                            <  1:291001/149‑1 (MQ=255)
                                                 gcgcgcATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCt                                                                                                     >  1:748465/1‑149 (MQ=255)
                                                    cgcATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGa                                                                                                  <  3:391025/149‑1 (MQ=255)
                                                    cgcATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGa                                                                                                  <  1:803085/149‑1 (MQ=255)
                                                       aTTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCGGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCTCTGGAGTa                                                                                               >  2:662722/1‑149 (MQ=255)
                                                                     ccTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCaa                                                                                 <  1:574260/149‑1 (MQ=255)
                                                                      cTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAg                                                                                <  2:9139/149‑1 (MQ=255)
                                                                        gtgtgAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGaa                                                                                        >  3:681590/1‑139 (MQ=255)
                                                                        gtgtgAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGaa                                                                                        <  4:681590/139‑1 (MQ=255)
                                                                         tgtgAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCAGGCGCTGGAGTACGCCGAAGCCGCAAGCGt                                                                             >  3:478320/1‑149 (MQ=255)
                                                                           tgAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTaa                                                                           >  4:591222/1‑149 (MQ=255)
                                                                                tttGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAAcgccg                                                                      <  3:591222/149‑1 (MQ=255)
                                                                                              aCCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGc                                                                               <  2:532728/126‑1 (MQ=255)
                                                                                              aCCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGc                                                                               >  1:532728/1‑126 (MQ=255)
                                                                                              aCCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCgatga                                                        >  2:799607/1‑149 (MQ=255)
                                                                                                  gTATGAGAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGcc                                                    <  2:535204/149‑1 (MQ=255)
                                                                                                  gTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGcc                                                    >  4:467595/1‑149 (MQ=255)
                                                                                                      gAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGa                                                                                                  >  4:194271/1‑99 (MQ=255)
                                                                                                      gAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGa                                                                                                  <  3:194271/99‑1 (MQ=255)
                                                                                                            cGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTAccc                                          >  1:164307/1‑149 (MQ=255)
                                                                                                                gAGCGCACCGGCAGCAGCGATAAACTGCTGGCGGTGGCTGACCGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGACTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACt                                      <  1:662722/149‑1 (MQ=255)
                                                                                                                   cgcACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGc                                   <  2:108082/149‑1 (MQ=255)
                                                                                                                      aCCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTgc                                >  2:27643/1‑149 (MQ=255)
                                                                                                                      aCCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTgc                                <  2:62075/149‑1 (MQ=255)
                                                                                                                          gcagcaGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCAt                            >  3:250756/1‑149 (MQ=255)
                                                                                                                                gcGATAAACTGCTGGCAGTGACTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGAc                      <  1:100192/149‑1 (MQ=255)
                                                                                                                                gcGATAAACTGCTGGCAGTGACTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGAc                      <  3:375054/149‑1 (MQ=255)
                                                                                                                                  gATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGc                    >  2:91534/1‑149 (MQ=255)
                                                                                                                                  gATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGCACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGc                    >  4:656762/1‑149 (MQ=255)
                                                                                                                                  gATAAACTGCTAGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGc                    >  3:516815/1‑149 (MQ=255)
                                                                                                                                   aTAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGCt                   >  4:800872/1‑149 (MQ=255)
                                                                                                                                    tAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGCTc                  >  3:222492/1‑149 (MQ=255)
                                                                                                                                         tgctgGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGCTCAGgcg             >  1:303283/1‑149 (MQ=255)
                                                                                                                                            tgGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGCTCAGGCGCTc          >  3:328193/1‑149 (MQ=255)
                                                                                                                                                  tGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGCTCAGGCGCTCACCGaa    >  4:512788/1‑149 (MQ=255)
                                                                                                                                                  tGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGCTCAGGCGCTCACCGaa    <  1:362541/149‑1 (MQ=255)
                                                                                                                                                   ggCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGCTCAGGCGCTCACCGAAc   <  3:698007/149‑1 (MQ=255)
                                                                                                                                                    gCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGCTCAGGCGCTCACCGAACt  >  1:648480/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAATGTGCCTGTGTGAGTTTTGTGTGGATATCACCAGTATGAAAGTCGCCGAGCGCACCGGCAGCAGCGATAAACTGCTGGCAGTGGCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGCTCAGGCGCTCACCGAACT  >  NC_000913/1835667‑1835963

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑