Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 325,159 C→G 5.2% P149A (CCA→GCA)  oprQ → outer‑membrane porin D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947325,1590CG 5.2% 246.0 / 4.0 77P149A (CCA→GCA) oprQouter‑membrane porin D
Reads supporting (aligned to +/- strand):  ref base C (38/34);  new base G (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

TGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGT  >  NC_002947/325011‑325305
                                                                                                                                                    |                                                                                                                                                  
tggCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGc                                                                                                                                                    >  4:249066/1‑149 (MQ=255)
tggCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGc                                                                                                                                                    <  7:362414/149‑1 (MQ=255)
       gTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtgct                                                                                                                                             <  3:249066/149‑1 (MQ=255)
        tCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtgctg                                                                                                                                            <  2:154536/149‑1 (MQ=255)
          gACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtgctg                                                                                                                                            >  8:243251/1‑147 (MQ=255)
           aCTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGc                                                                                                                                                    >  1:11116/1‑138 (MQ=255)
           aCTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGc                                                                                                                                                    <  2:11116/138‑1 (MQ=255)
              tcttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctggc                                                                                                                                          >  5:311476/1‑143 (MQ=255)
              tcttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctggc                                                                                                                                          <  6:311476/145‑3 (MQ=255)
                ttcAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTacga                                                                                                                                    >  4:339431/1‑149 (MQ=255)
                   aaGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtgctg                                                                                                                                            >  1:224493/1‑138 (MQ=255)
                   aaGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtgctg                                                                                                                                            <  2:224493/138‑1 (MQ=255)
                     gCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATg                                                                                                                               <  5:382437/149‑1 (MQ=255)
                      ccGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctggc                                                                                                                                          <  4:218334/137‑3 (MQ=255)
                      ccGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGGCAGtgctggc                                                                                                                                          >  3:218334/1‑135 (MQ=255)
                       cGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGCGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCa                                                                                                                             >  5:264814/1‑149 (MQ=255)
                          aagAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACg                                                                                                                           >  7:365037/3‑148 (MQ=255)
                           cGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTc                                                                                                                         <  7:243463/149‑1 (MQ=255)
                            gAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTct                                                                                                                        >  2:334346/1‑149 (MQ=255)
                               accaccAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGGTGTGAAGTTCCGCATTATCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctgct                                                                                                                     <  1:137258/149‑1 (MQ=255)
                               accaccAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtgctg                                                                                                                                            <  8:17970/126‑1 (MQ=255)
                               accaccAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtgctg                                                                                                                                            >  7:17970/1‑126 (MQ=255)
                               accaccAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGCCGATGCACGTctgct                                                                                                                     >  5:351231/1‑149 (MQ=255)
                               accaccAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctgct                                                                                                                     <  3:339431/149‑1 (MQ=255)
                                ccaccaATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTctgctg                                                                                                                    >  1:26558/1‑149 (MQ=255)
                                 caccaATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCtgctgc                                                                                                                   >  5:70421/1‑149 (MQ=255)
                                   ccaATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCa                                                                                                                 <  3:350944/149‑1 (MQ=255)
                                    caATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCtgctgc                                                                                                                   >  7:15244/1‑146 (MQ=255)
                                    caATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCtgctgc                                                                                                                   >  6:408889/1‑146 (MQ=255)
                                    caATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCtgctgc                                                                                                                   >  5:185291/1‑146 (MQ=255)
                                     aaTGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAg                                                                                                                >  5:65536/1‑148 (MQ=255)
                                         ccGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGc                                                                                                           >  2:190574/1‑149 (MQ=255)
                                            aaGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTc                                                                                                        <  6:185291/149‑1 (MQ=255)
                                                         cTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGa                                                                                           <  7:243251/149‑1 (MQ=255)
                                                           cgccgcACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGAt                                                                                          >  6:257178/1‑148 (MQ=255)
                                                            gccgcACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATc                                                                                         <  8:365037/148‑1 (MQ=255)
                                                             ccgcACAACCTGGCGCGTGGTGGTGCTGCTGTGAACTTCCGCATTTCCAACACCGTCCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCAGATAAc                                                                                       >  7:88669/1‑149 (MQ=255)
                                                               gcACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCt                                                                                     <  5:233279/149‑1 (MQ=255)
                                                               gcACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCt                                                                                     >  7:181952/1‑149 (MQ=255)
                                                                      cTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGCGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAaga                                                                              <  6:264814/149‑1 (MQ=255)
                                                                          cgcgTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATc                                                                          <  8:15244/149‑1 (MQ=255)
                                                                             gtggtggtgCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAg                                                                       <  1:334346/149‑1 (MQ=255)
                                                                              tggtggtgCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAgg                                                                      <  6:65536/149‑1 (MQ=255)
                                                                               ggtggtgATGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGGAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAggg                                                                     <  8:88669/149‑1 (MQ=255)
                                                                                  ggtgCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCt                                                                  >  8:285367/1‑149 (MQ=255)
                                                                                      ctgctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGaa                                                              <  5:87945/149‑1 (MQ=255)
                                                                                      ctgctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGaa                                                              >  6:281569/1‑149 (MQ=255)
                                                                                            tgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCTAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAAc                                                        >  2:132259/1‑149 (MQ=255)
                                                                                            tgAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGa                                                          <  4:7431/147‑1 (MQ=255)
                                                                                                 ttCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACgccgg                                                   >  4:86076/1‑149 (MQ=255)
                                                                                                     gCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCCCCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGc                                               >  8:113890/1‑149 (MQ=255)
                                                                                                     gCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGaa                                                              >  1:318819/1‑134 (MQ=255)
                                                                                                     gCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGaa                                                              <  2:318819/134‑1 (MQ=255)
                                                                                                        tttCCAACACCGTGCTGAAGGACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTc                                            <  3:144078/149‑1 (MQ=255)
                                                                                                           tggtcacCGTACTTCAGCACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAcc                                         >  2:324348/5‑149 (MQ=255)
                                                                                                           ccAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAc                                          <  6:354611/148‑1 (MQ=255)
                                                                                                             aaCACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCCAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCa                                       <  7:210693/149‑1 (MQ=255)
                                                                                                             aaCACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAcc                                         >  2:322964/1‑147 (MQ=255)
                                                                                                              tcacCGTACTTCAGCACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAg                                      >  5:83240/2‑149 (MQ=255)
                                                                                                              acacCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAg                                      >  3:39983/1‑149 (MQ=255)
                                                                                                              acacCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAg                                      <  4:39983/149‑1 (MQ=255)
                                                                                                              acacCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAg                                      <  8:181952/149‑1 (MQ=255)
                                                                                                               cacCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTTCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAg                                      >  2:82787/1‑148 (MQ=255)
                                                                                                                acCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAgg                                     >  8:249355/1‑148 (MQ=255)
                                                                                                                 ccGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGaa                                   <  5:408889/149‑1 (MQ=255)
                                                                                                                   gTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgc                                 <  5:281569/149‑1 (MQ=255)
                                                                                                                         aaGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCa                            >  4:145029/1‑148 (MQ=255)
                                                                                                                          aGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgcgcgc                             >  8:346685/1‑146 (MQ=255)
                                                                                                                          aGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAag                          <  7:285367/149‑1 (MQ=255)
                                                                                                                           gTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAaga                         <  7:275616/149‑1 (MQ=255)
                                                                                                                              cGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagc                       <  6:70421/148‑1 (MQ=255)
                                                                                                                                    ccAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCACAAgagcgccgagc                >  1:418298/1‑149 (MQ=255)
                                                                                                                                    ccAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACACATGAAGCGCGCAAgagcgactagc                >  6:81878/1‑149 (MQ=255)
                                                                                                                                    ccAGATGCCAGCACTGACAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGGCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccgag                 <  6:353366/148‑1 (MQ=255)
                                                                                                                                     cAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgacgagcg               >  8:83365/1‑149 (MQ=255)
                                                                                                                                            cAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAACGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAt        <  6:433683/149‑1 (MQ=255)
                                                                                                                                            cAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAt        <  7:346685/149‑1 (MQ=255)
                                                                                                                                            cAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAt        >  3:207090/1‑149 (MQ=255)
                                                                                                                                             aGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtg       >  5:292338/1‑149 (MQ=255)
                                                                                                                                             aGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAtg       <  3:86076/149‑1 (MQ=255)
                                                                                                                                                  tGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGGGAtggtggt  >  7:170535/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGCGGCGTCGACTTCTTCAAGCCGCACGAAACCACCAATGCCGAAGGCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGT  >  NC_002947/325011‑325305

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