Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,802,992 A→T 8.5% K459M (AAG→ATG) ‡ rbsA‑I → ribose ABC transporter ‑ ATP‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,802,9920AT 8.5% 173.3 / 5.4 59K459M (AAG→ATG) ‡rbsA‑Iribose ABC transporter ‑ ATP‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base T (3/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGT  >  NC_002947/2802846‑2803137
                                                                                                                                                   |                                                                                                                                                 
tGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCcaaggcaagg                                                                                                                                                 >  2:282405/1‑149 (MQ=255)
    ggTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCCGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCta                                                                                                                                             <  3:237884/149‑2 (MQ=255)
        aCCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCACCAAGGCCTTGCCTtggc                                                                                                                                          <  8:159293/148‑2 (MQ=255)
           agcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                      <  2:207163/149‑1 (MQ=255)
           agcagAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGCGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                      <  5:98269/149‑1 (MQ=255)
           agcagAAGGTGGTGATCGGCCGCTAGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTtggtggtg                                                                                                                                      <  3:383712/149‑1 (MQ=255)
                  ggtggtGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCa                                                                                                                               <  1:282405/149‑1 (MQ=255)
                     ggtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTg                                                                                                                            <  4:202661/149‑1 (MQ=255)
                      gtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGa                                                                                                                           <  7:281782/149‑1 (MQ=255)
                      gtgATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGa                                                                                                                           >  4:165612/1‑149 (MQ=255)
                           cGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTg                                                                                                                       >  1:304388/1‑149 (MQ=255)
                                gctggctgGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGaa                                                                                                                 <  1:186635/149‑1 (MQ=255)
                                     ctggAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGAt                                                                                                            <  6:309490/149‑1 (MQ=255)
                                              gATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTg                                                                                                   >  3:380777/1‑149 (MQ=255)
                                              gATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTg                                                                                                   >  7:98680/1‑149 (MQ=255)
                                                     aGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCg                                                                                             >  2:25096/1‑148 (MQ=255)
                                                      ggTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGTCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCa                                                                                           >  5:365585/1‑149 (MQ=255)
                                                       gTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCAt                                                                                          >  5:321955/1‑149 (MQ=255)
                                                        tgctgcTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATc                                                                                         <  4:380777/149‑1 (MQ=255)
                                                         gctgctGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCg                                                                                        >  3:89019/1‑149 (MQ=255)
                                                         gctgctGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCACCAAGGCCTTGCCTtggc                                                                                                                                          <  4:269275/99‑2 (MQ=255)
                                                         gctgctGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCACCAAGGCCTTGCCTtggc                                                                                                                                          >  3:269275/1‑98 (MQ=255)
                                                          ctgctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCg                                                                                        <  3:165612/148‑1 (MQ=255)
                                                          ctgctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGc                                                                                       >  3:312261/1‑149 (MQ=255)
                                                           tgctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGcc                                                                                      >  7:215787/1‑149 (MQ=255)
                                                            gctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCg                                                                                     >  5:297862/1‑149 (MQ=255)
                                                             ctgTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGcc                                                                                      <  4:89019/147‑1 (MQ=255)
                                                                  cgacgaACCCACGCGTGGCATCTATGTGGGCGCCAAGTGCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTg                                                                                <  2:294649/148‑1 (MQ=255)
                                                                       aaCCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCt                                                                          <  6:321955/149‑1 (MQ=255)
                                                                        aCCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGATGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTg                                                                         <  3:303680/149‑1 (MQ=255)
                                                                          ccACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGAttt                                                                                                    <  2:281603/120‑1 (MQ=255)
                                                                          ccACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGAttt                                                                                                    >  1:281603/1‑120 (MQ=255)
                                                                          ccACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGgc                                                                       >  1:33280/1‑149 (MQ=255)
                                                                               gcgTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGgccgcc                                                                   <  1:305814/149‑1 (MQ=255)
                                                                                 gTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTg                                                                 <  1:25096/149‑1 (MQ=255)
                                                                                  tGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGa                                                                >  2:278621/1‑149 (MQ=255)
                                                                                    gCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGAt                                                               <  2:304388/148‑1 (MQ=255)
                                                                                    gCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTCATc                                                              >  8:304614/1‑149 (MQ=255)
                                                                                     cATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCg                                                             <  8:215787/149‑1 (MQ=255)
                                                                                         gATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGacac                                                         <  2:33280/149‑1 (MQ=255)
                                                                                                 cgcCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTcgcgc                                                 >  2:321104/1‑149 (MQ=255)
                                                                                                  gcCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGac                                                           >  8:17277/1‑138 (MQ=255)
                                                                                                  gcCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGac                                                           <  7:17277/138‑1 (MQ=255)
                                                                                                   ccAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCAGCTGGCCGCCTGATCGACACCTTCGCGCGt                                               >  5:316746/1‑149 (MQ=255)
                                                                                                       gTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATc                                           <  6:143058/149‑1 (MQ=255)
                                                                                                        ttCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCa                                          >  2:124455/1‑149 (MQ=255)
                                                                                                         tCGACATTTATGGCCTGCAGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGt                                               <  8:163557/143‑1 (MQ=255)
                                                                                                         tCGACATTTATGGCCTGCAGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGt                                               >  7:163557/1‑143 (MQ=255)
                                                                                                           gACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTg                                       <  1:321104/149‑1 (MQ=255)
                                                                                                                      ccTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTc                                                     <  1:20082/124‑1 (MQ=255)
                                                                                                                      ccTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTc                                                     >  2:20082/1‑124 (MQ=255)
                                                                                                                       ctgctgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGAc                             >  1:247593/1‑147 (MQ=255)
                                                                                                                         gctgGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGt                         >  8:313305/1‑149 (MQ=255)
                                                                                                                             gCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGttg                       <  7:313305/147‑1 (MQ=255)
                                                                                                                               cGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTgg                                                                        <  2:182431/96‑1 (MQ=255)
                                                                                                                               cGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTgg                                                                        >  1:182431/1‑96 (MQ=255)
                                                                                                                                     ccaccaccAAGGCCTTGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTgccgccg             >  6:28279/7‑149 (MQ=255)
                                                                                                                                      caccaccAAGGCCTTGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTgccgccgc            >  1:50140/6‑149 (MQ=255)
                                                                                                                                        gcgcCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTgccgccgc            >  1:146344/1‑147 (MQ=255)
                                                                                                                                        gcgcCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCt          >  5:132847/1‑149 (MQ=255)
                                                                                                                                              aggcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCg    <  6:297862/149‑1 (MQ=255)
                                                                                                                                                gcaaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGt  >  7:556/1‑149 (MQ=255)
                                                                                                                                                 caaggcCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGt  >  5:52687/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGGCGGTAACCAGCAGAAGGTGGTGATCGGCCGCTGGCTGGAGCGCGATTGCCAGGTGCTGCTGTTCGACGAACCCACGCGTGGCATCGATGTGGGCGCCAAGTTCGACATTTATGGCCTGCTGGCCGAGCTGGCGCGCCAAGGCAAGGCCTTGGTGGTGGTGTCCAGTGACCTGCGCGAACTGATGCTGATTTGCGACCGCATCGCCGTGCTGTCCGCTGGCCGCCTGATCGACACCTTCGCGCGTGATCATTGGAGCCAGGACCAGTTGCTTGCCGCCGCCTTCGCCGGT  >  NC_002947/2802846‑2803137

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: