Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,597,890 A→G 11.1% D464G (GAC→GGC) ‡ PP_4923 → outer membrane efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,597,8900AG11.1% 163.2 / 10.1 63D464G (GAC→GGC) ‡PP_4923outer membrane efflux protein
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base G (4/3);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

CAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAATA  >  NC_002947/5597750‑5598038
                                                                                                                                            |                                                                                                                                                     
caacaaTAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTggcgg                                                                                                                                               <  6:196540/148‑1 (MQ=255)
   caaTAGCCGTTACGACTACATTCTCGACACCCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGgcag                                                                                                                                            <  2:36196/149‑1 (MQ=255)
        gCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCagcagc                                                                                                                                        <  2:58585/148‑1 (MQ=255)
        gCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCACGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCg                                                                                                                                       <  6:18530/149‑1 (MQ=255)
         ccGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCgg                                                                                                                                      <  6:225487/149‑1 (MQ=255)
           gTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGc                                                                                                                                     >  2:270347/1‑148 (MQ=255)
            ttACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGcc                                                                                                                                    >  1:266177/1‑148 (MQ=255)
                  cTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCCGCTGCTGCCGCCaggtccg                                                                                                                             <  5:113480/149‑7 (MQ=255)
                    aCATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTc                                                                                                                           >  7:25532/1‑149 (MQ=255)
                     cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTc                                                                                                                           <  8:25532/148‑1 (MQ=255)
                     cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTc                                                                                                                           <  2:266177/148‑1 (MQ=255)
                     cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgcag                                                                                                                                            <  4:6410/131‑1 (MQ=255)
                     cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgcag                                                                                                                                            >  3:6410/1‑131 (MQ=255)
                        tctcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCCGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAgc                                                                                                                       >  4:127426/1‑149 (MQ=255)
                        tctcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgca                                                                                                                                             <  1:239439/127‑1 (MQ=255)
                        tctcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgca                                                                                                                                             >  2:239439/1‑127 (MQ=255)
                            gACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAg                                                                                                                   <  5:134370/149‑1 (MQ=255)
                                 ccTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGa                                                                                                                         >  6:108853/1‑138 (MQ=255)
                                 ccTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGa                                                                                                                         <  5:108853/138‑1 (MQ=255)
                                     aGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCtt                                                                                                              >  2:153067/1‑145 (MQ=255)
                                           aaGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACACCCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGAAAgggg                                                                                                    >  3:207729/1‑149 (MQ=255)
                                           aaGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCATGCCTTGAGATAgg                                                                                                      <  8:148856/147‑1 (MQ=255)
                                           aaGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAgg                                                                                                      <  2:262123/147‑1 (MQ=255)
                                            tcacggcagCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGc                                                                                                   >  4:95698/5‑149 (MQ=255)
                                                 gcagCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGAc                                                                                              >  8:71559/1‑149 (MQ=255)
                                                      cGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTggcgg                                                                                                                                               <  7:34865/95‑1 (MQ=255)
                                                      cGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTggcgg                                                                                                                                               >  8:34865/1‑95 (MQ=255)
                                                        gCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTcgc                                                                                       >  4:210738/1‑149 (MQ=255)
                                                            ccTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGAc                                                                                                                                                      <  6:276018/82‑1 (MQ=255)
                                                            ccTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGAc                                                                                                                                                      >  5:276018/1‑82 (MQ=255)
                                                            ccTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCgg                                                                                     <  4:27662/147‑1 (MQ=255)
                                                            ccTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCgg                                                                                     >  3:27662/1‑147 (MQ=255)
                                                            ccTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGc                                                                                   >  5:98798/1‑149 (MQ=255)
                                                                   ccGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCg                                                                            >  7:63629/1‑149 (MQ=255)
                                                                     gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCt                                                                         >  8:216040/1‑149 (MQ=255)
                                                                     gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgca                                                                                                                                             <  4:63041/82‑1 (MQ=255)
                                                                     gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgca                                                                                                                                             >  3:63041/1‑82 (MQ=255)
                                                                      cAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTGCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTc                                                                        <  6:214222/149‑1 (MQ=255)
                                                                               cAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACgg                                                                <  6:98798/148‑1 (MQ=255)
                                                                                  gACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGcc                                                            >  4:56969/1‑149 (MQ=255)
                                                                                     cTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAgg                                                         >  2:218322/1‑149 (MQ=255)
                                                                                      tGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCGTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGc                                                        >  6:171836/1‑149 (MQ=255)
                                                                                       gAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGcc                                                       >  2:266463/1‑149 (MQ=255)
                                                                                       gAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGcc                                                       >  6:173426/1‑149 (MQ=255)
                                                                                                ccTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGc                                                <  7:216040/147‑1 (MQ=255)
                                                                                                 cTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCCTCCCACACGGTGCCAGGCCTGAAGGCat                                              <  1:153067/149‑1 (MQ=255)
                                                                                                 cTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCat                                              >  7:147379/1‑148 (MQ=255)
                                                                                                       ccTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAc                                       >  8:68624/1‑149 (MQ=255)
                                                                                                        cTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAcc                                      >  4:119758/1‑149 (MQ=255)
                                                                                                        cTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGACCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAcc                                      >  2:167644/1‑149 (MQ=255)
                                                                                                           aCTACAACCCCGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCTCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCtgt                                   <  3:127426/149‑1 (MQ=255)
                                                                                                               cAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTga                                 <  3:56969/147‑1 (MQ=255)
                                                                                                                aaCCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGagag                              >  3:65641/1‑149 (MQ=255)
                                                                                                                   ccGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTg                           <  3:210738/149‑1 (MQ=255)
                                                                                                                       aCAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCGTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCtt                       <  5:171836/149‑1 (MQ=255)
                                                                                                                         aaGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGc                     >  4:109718/1‑149 (MQ=255)
                                                                                                                               ttCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGa               >  1:30903/1‑149 (MQ=255)
                                                                                                                                 ccTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGAt              <  1:266463/148‑1 (MQ=255)
                                                                                                                                  cTGCCGCCAGGTCCGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAg            >  6:185341/1‑149 (MQ=255)
                                                                                                                                   tGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTg                           <  3:221126/133‑1 (MQ=255)
                                                                                                                                   tGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTg                           >  4:221126/1‑133 (MQ=255)
                                                                                                                                         cgGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCa     <  4:65641/149‑1 (MQ=255)
                                                                                                                                          ggACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCaa    <  8:63629/149‑1 (MQ=255)
                                                                                                                                            aCCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAata  <  7:71559/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
CAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAATA  >  NC_002947/5597750‑5598038

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: