Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,835,576 G→A 5.5% R21H (CGT→CAT)  ynjA → carboxymuconolactone decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,835,5760GA 5.5% 230.7 / 2.3 73R21H (CGT→CAT) ynjAcarboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base G (38/31);  new base A (2/2);  total (40/33)
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.96e-01

ACCGGACTGTTTATTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGC  >  NC_000913/1835429‑1835719
                                                                                                                                                   |                                                                                                                                               
aCCGGACTGTTTATTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCAtcgt                                                                                                                                                >  2:42235/1‑149 (MQ=255)
       tGTTTATTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATAAGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCAGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCAt                                                                                                                                          >  1:362303/1‑148 (MQ=255)
            aTTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGc                                                                                                                                    >  1:131045/1‑149 (MQ=255)
            aTTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGc                                                                                                                                    >  6:3273/1‑149 (MQ=255)
             ttCTGTTTGCACTGACCGTGGTGATTTTTATGGAGAAGAAAATATGGCTTTAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGAATCGTCGCCATTACgtc                                                                                                                                    >  2:38125/1‑146 (MQ=255)
                tGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCATCGCCATTACGGCGAgg                                                                                                                                >  5:423361/1‑149 (MQ=255)
                  tttGCCCTGACCGTGTGGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAg                                                                                                                                 <  2:362303/146‑1 (MQ=255)
                  tttGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAgg                                                                                                                                <  4:83736/147‑1 (MQ=255)
                    tGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGACCGGCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCGGAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGc                                                                                                                             >  8:117273/1‑148 (MQ=255)
                      cccTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTaa                                                                                                                          <  1:42235/149‑1 (MQ=255)
                      cccAGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCt                                                                                                                            >  7:235936/1‑147 (MQ=255)
                       ccTGACCGTTGTGATTTTTATGGCGAAGAAAATATGGCTTGCACGCCAGAAGAGGAATGCCTGATCGGTTTACCGCCGCTTAGCAAAATTACTTTAATTTTACGTCCACAGTCGTGGCTGCATCGTCGCCATTACGGAGAGCTGCTAAg                                                                                                                         >  4:22919/1‑149 (MQ=255)
                       ccTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAg                                                                                                                         >  2:16244/1‑149 (MQ=255)
                         tGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGcc                                                                                                                       >  1:321248/1‑149 (MQ=255)
                           aCCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGccc                                                                                                                      >  2:222692/1‑148 (MQ=255)
                               tggtgATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTc                                                                                                                 >  3:262024/1‑149 (MQ=255)
                                   gATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGc                                                                                                                        <  2:33671/138‑1 (MQ=255)
                                   gATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGc                                                                                                                        >  1:33671/1‑138 (MQ=255)
                                       tttATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCtggtg                                                                                                          <  2:74182/148‑1 (MQ=255)
                                          aTGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTggggg                                                                                                      >  8:258589/1‑149 (MQ=255)
                                               gaagaaAATATGGATTGAACGCCAGAATAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTGCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCCTTCGCTCGTGGGGGCggct                                                                                                 >  2:306533/1‑147 (MQ=255)
                                                aagaaAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATc                                                                                                >  3:137653/1‑149 (MQ=255)
                                                aagaaAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATc                                                                                                <  8:417658/149‑1 (MQ=255)
                                                 agaaAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATcc                                                                                               <  4:337536/149‑1 (MQ=255)
                                                 agaaAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATcc                                                                                               >  4:176484/1‑149 (MQ=255)
                                                   aaaaTATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCg                                                                                             <  7:258589/149‑1 (MQ=255)
                                                     aaTATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtt                                                                                           <  2:321248/149‑1 (MQ=255)
                                                      atatGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtt                                                                                           >  2:111484/1‑148 (MQ=255)
                                                         tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatc                                                                                       <  4:262024/149‑1 (MQ=255)
                                                          ggCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatct                                                                                      >  6:407110/1‑149 (MQ=255)
                                                           gCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatctt                                                                                     <  2:369127/149‑1 (MQ=255)
                                                                aaCGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCATCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatcttttat                                                                                 <  6:423361/148‑1 (MQ=255)
                                                                  cGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTg                                                                              <  6:400170/149‑1 (MQ=255)
                                                                          gagGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAt                                                                      <  2:131045/149‑1 (MQ=255)
                                                                          gacgactGCCTGATGGGGTTACCGCCGCTTAGCAAAATTCCTTTAATTATACGTCCACAGGCGTGGCTGCATCGTCGCCATTACCGCGAGGTTCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAt                                                                      <  2:445313/143‑1 (MQ=255)
                                                                            ggAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgt                                                                    <  3:22919/149‑1 (MQ=255)
                                                                             gAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtt                                                                   >  5:73407/1‑149 (MQ=255)
                                                                             gAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAt                                                                      >  3:179221/1‑146 (MQ=255)
                                                                             gAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAt                                                                      <  4:179221/146‑1 (MQ=255)
                                                                                   cTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTtgg                                                             <  1:16244/149‑1 (MQ=255)
                                                                                      aTGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttg                                                                 <  1:433484/142‑1 (MQ=255)
                                                                                      aTGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttg                                                                 >  2:433484/1‑142 (MQ=255)
                                                                                         ggTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCGACAGGCGTGGCTGCATCGTCGCCATTAAGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGTGCGGATCCCGTTTATATTTTATCTGTTGTCGATGTTTGTtggctggc                                                        <  3:241343/148‑1 (MQ=255)
                                                                                           tttACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGgctggctg                                                      >  4:368759/1‑148 (MQ=255)
                                                                                              aCCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAg                                                   <  6:73407/148‑1 (MQ=255)
                                                                                               ccgccgCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAgcg                                                 >  1:37537/1‑149 (MQ=255)
                                                                                                 gccgcTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGT‑GGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCa                                               >  4:131705/1‑148 (MQ=255)
                                                                                                  ccgcTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCaa                                              <  5:3273/149‑1 (MQ=255)
                                                                                                        aGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTtgg                                                             >  4:237202/1‑128 (MQ=255)
                                                                                                        aGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTtgg                                                             <  3:237202/128‑1 (MQ=255)
                                                                                                         gCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCAt                                                                                                                   <  6:327004/73‑1 (MQ=255)
                                                                                                         gCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCAt                                                                                                                   >  5:327004/1‑73 (MQ=255)
                                                                                                         gCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCa                                       <  8:235936/149‑1 (MQ=255)
                                                                                                           aaaaTTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCAcc                                     >  4:271871/1‑149 (MQ=255)
                                                                                                           aaaaTTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCAcc                                     >  3:16056/1‑149 (MQ=255)
                                                                                                                   tttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatcttttatct                                                                               <  5:156536/99‑1 (MQ=255)
                                                                                                                   tttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatcttttatct                                                                               >  6:156536/1‑99 (MQ=255)
                                                                                                                   tttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTg                                                                              <  5:91871/100‑1 (MQ=255)
                                                                                                                   tttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTg                                                                              >  6:91871/1‑100 (MQ=255)
                                                                                                                   tttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATc                             >  1:73490/1‑149 (MQ=255)
                                                                                                                   tttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATc                             >  4:431611/1‑149 (MQ=255)
                                                                                                                      aaTTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTc                                 <  7:186611/142‑1 (MQ=255)
                                                                                                                      aaTTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTc                                 >  8:186611/1‑142 (MQ=255)
                                                                                                                       aTTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCgg                          <  3:368759/148‑1 (MQ=255)
                                                                                                                        ttttACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGATGTCGATGTTTGTTGGCTGGCTGGAGCTCAAACGCTCACCGCTCGATCCGgta                        <  1:222692/149‑1 (MQ=255)
                                                                                                                                      gCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGt          >  3:274050/1‑149 (MQ=255)
                                                                                                                                        gTGGCTGCAGCGTCGCCATTACGGCGATGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCa        <  1:306533/149‑1 (MQ=255)
                                                                                                                                           gCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTTGAGAGCAAACGCTCACCTCTCGATCCGGTAGTACTATCGCTTGTCAgct     >  6:25067/1‑148 (MQ=255)
                                                                                                                                             tGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAgcgcg   >  6:391899/1‑149 (MQ=255)
                                                                                                                                             tGCATCATCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATc                             >  2:4271/1‑123 (MQ=255)
                                                                                                                                             tGCATCATCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATc                             <  1:4271/123‑1 (MQ=255)
                                                                                                                                              gCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAgcgcgc  >  7:267414/1‑149 (MQ=255)
                                                                                                                                                aTCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGgtagta                     <  4:18766/128‑1 (MQ=255)
                                                                                                                                                aTCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGgtagta                     >  3:18766/1‑128 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACCGGACTGTTTATTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGC  >  NC_000913/1835429‑1835719

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