Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,597,893 T→C 11.1% L465P (CTG→CCG)  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,8930TC11.1% 164.4 / 10.6 63L465P (CTG→CCG) PP_4923outer membrane efflux protein
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base C (4/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAATAT  >  NC_002947/5597750‑5598039
                                                                                                                                               |                                                                                                                                                   
caacaaTAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACC‑GAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTggcgg                                                                                                                                                <  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                                                                                                                                             >  3:6410/1‑131 (MQ=255)
                     cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgcag                                                                                                                                             <  4:6410/131‑1 (MQ=255)
                        tctcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCCGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAgc                                                                                                                        >  4:127426/1‑149 (MQ=255)
                        tctcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgca                                                                                                                                              >  2:239439/1‑127 (MQ=255)
                        tctcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgca                                                                                                                                              <  1:239439/127‑1 (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)
                                                            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)
                                                                     gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑Ct                                                                          >  8:216040/1‑149 (MQ=255)
                                                                     gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgca                                                                                                                                              <  4:63041/82‑1 (MQ=255)
                                                                     gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgca                                                                                                                                              >  3:63041/1‑82 (MQ=255)
                                                                      cAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTGCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTc                                                                         <  6:214222/149‑1 (MQ=255)
                                                                               cAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACgg                                                                 <  6:98798/148‑1 (MQ=255)
                                                                                  gACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGcc                                                             >  4:56969/1‑149 (MQ=255)
                                                                                     cTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAgg                                                          >  2:218322/1‑149 (MQ=255)
                                                                                      tGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCGTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGc                                                         >  6:171836/1‑149 (MQ=255)
                                                                                       gAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGcc                                                        >  6:173426/1‑149 (MQ=255)
                                                                                       gAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGcc                                                        >  2:266463/1‑149 (MQ=255)
                                                                                                ccTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGc                                                 <  7:216040/147‑1 (MQ=255)
                                                                                                 cTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCCTCCCACACGGTGCCAGGCCTGAAGGCat                                               <  1:153067/149‑1 (MQ=255)
                                                                                                 cTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCat                                               >  7:147379/1‑148 (MQ=255)
                                                                                                       ccTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAc                                        >  8:68624/1‑149 (MQ=255)
                                                                                                        cTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAcc                                       >  4:119758/1‑149 (MQ=255)
                                                                                                        cTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGACCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAcc                                       >  2:167644/1‑149 (MQ=255)
                                                                                                           aCTACAACCCCGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCTCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCtgt                                    <  3:127426/149‑1 (MQ=255)
                                                                                                               cAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTga                                  <  3:56969/147‑1 (MQ=255)
                                                                                                                aaCCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGagag                               >  3:65641/1‑149 (MQ=255)
                                                                                                                   ccGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTg                            <  3:210738/149‑1 (MQ=255)
                                                                                                                       aCAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCGTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCtt                        <  5:171836/149‑1 (MQ=255)
                                                                                                                         aaGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGc                      >  4:109718/1‑149 (MQ=255)
                                                                                                                               ttCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGa                >  1:30903/1‑149 (MQ=255)
                                                                                                                                 ccTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGAt               <  1:266463/148‑1 (MQ=255)
                                                                                                                                  cTGCCGCCAGGTCCGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAg             >  6:185341/1‑149 (MQ=255)
                                                                                                                                   tGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTg                            <  3:221126/133‑1 (MQ=255)
                                                                                                                                   tGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTg                            >  4:221126/1‑133 (MQ=255)
                                                                                                                                        aagttcttGGCCGCTGCTGCCGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCtt        >  3:96385/8‑148 (MQ=255)
                                                                                                                                         cgGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCa      <  4:65641/149‑1 (MQ=255)
                                                                                                                                          ggACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCaa     <  8:63629/149‑1 (MQ=255)
                                                                                                                                            aCCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAata   <  7:71559/149‑1 (MQ=255)
                                                                                                                                             ccTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAatat  <  3:109718/149‑1 (MQ=255)
                                                                                                                                               tGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATa              <  7:254360/135‑1 (MQ=255)
                                                                                                                                               tGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATa              >  8:254360/1‑135 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAATAT  >  NC_002947/5597750‑5598039

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: ‑