Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,893,484 G→C 6.9% G283A (GGC→GCC)  rarD → chloramphenicol‑sensitive protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,893,4840GC 6.9% 281.6 / 8.0 87G283A (GGC→GCC) rarDchloramphenicol‑sensitive protein
Reads supporting (aligned to +/- strand):  ref base G (41/40);  new base C (4/2);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

TGCCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAG  >  NC_002947/3893337‑3893625
                                                                                                                                                     |                                                                                                                                              
tgcCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGc                                                                                                                                               >  2:393440/1‑149 (MQ=255)
tgcCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGc                                                                                                                                               >  8:110977/1‑149 (MQ=255)
       cAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAg                                                                                                                                         <  2:115584/148‑1 (MQ=255)
        aGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGc                                                                                                                                        <  2:303693/148‑1 (MQ=255)
         gCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAc                                                                                                                                      >  4:362427/1‑149 (MQ=255)
         gCCTGTTCGGCATTCTCATCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAc                                                                                                                                      <  7:185545/149‑1 (MQ=255)
          ccTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAAGCCCTCAAGCAcc                                                                                                                                     <  1:426678/149‑1 (MQ=255)
            tGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTg                                                                                                                                   >  1:99564/1‑149 (MQ=255)
              ttCGGGCTGCTCAGCTACGTCGAACCCGTTCTTCTGGTTTGCGTGGCGTTGCTTCTGGGCGAAAACATCGGGCCGGATCAGTGGTCTTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCt                                                                                                                                    <  1:348361/148‑1 (MQ=255)
              ttCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTg                                                                                                                                   <  5:85644/147‑1 (MQ=255)
                   ccTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGt                                                                                                                            >  3:249595/1‑148 (MQ=255)
                          aGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCg                                                                                                                    >  8:166003/1‑149 (MQ=255)
                            cTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCCCCTGCTACGTCAGCGTCGcc                                                                                                                  >  6:303821/1‑149 (MQ=255)
                            cTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGcc                                                                                                                  >  8:139320/1‑149 (MQ=255)
                              aCGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGATGGGCTAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATATGGGCCGCAGTACTGCTGCTGTTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGt                                                                                                                >  4:12605/1‑149 (MQ=255)
                                gTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGTCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTc                                                                                                              <  6:23728/149‑1 (MQ=255)
                                 tCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCg                                                                                                             <  2:99564/149‑1 (MQ=255)
                                 tCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCg                                                                                                             >  6:30070/1‑149 (MQ=255)
                                     aCCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGgtg                                                                                                          >  1:236900/1‑148 (MQ=255)
                                      cccGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGgtgt                                                                                                         >  3:142072/1‑148 (MQ=255)
                                       ccGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTaa                                                                                                       <  1:393440/149‑1 (MQ=255)
                                         gTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGc                                                                                                     >  3:159538/1‑149 (MQ=255)
                                         gTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTACTACGTCAGCGTCGCCGTTCGGTGTAAGc                                                                                                     <  1:32423/149‑1 (MQ=255)
                                               ctgGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAAcc                                                                                               >  2:263315/1‑149 (MQ=255)
                                               ctgGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAAc                                                                                                <  1:220567/148‑1 (MQ=255)
                                                  gTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTgg                                                                                            <  1:317247/149‑1 (MQ=255)
                                                  gTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTg                                                                                             <  7:125983/148‑1 (MQ=255)
                                                   tAGGCGTTGCTTTGCTGCTGGGCGAAACCTTCGGCCCGTATCAGTTGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGTGCTTCAAGCACCTGCTACGTCAGCGTCGCGGTTCGGTGTAAGCCGAACCTGGc                                                                                           <  2:48840/149‑1 (MQ=255)
                                                             gtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCg                                                                                 <  7:139320/149‑1 (MQ=255)
                                                                 ctgctgGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGctct                                                                             >  4:5446/1‑149 (MQ=255)
                                                                  tgctgGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCtt                                                                            >  3:231118/1‑149 (MQ=255)
                                                                    ctgGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTc                                                                           >  6:152765/1‑148 (MQ=255)
                                                                    ctgGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCa                                                                          <  4:208340/149‑1 (MQ=255)
                                                                     tgGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTaa                                                                                                       <  5:232825/119‑1 (MQ=255)
                                                                     tgGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTaa                                                                                                       >  6:232825/1‑119 (MQ=255)
                                                                     tgGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCaa                                                                         >  8:52129/1‑149 (MQ=255)
                                                                      gggCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGtt                                                                                                               >  7:206073/1‑110 (MQ=255)
                                                                      gggCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGtt                                                                                                               <  8:206073/110‑1 (MQ=255)
                                                                      gggCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAAc                                                                        >  7:340645/1‑149 (MQ=255)
                                                                          gAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGTCGAACCTGGCGTACCCGTCGCTCTTCAACCGcc                                                                    >  4:340019/1‑149 (MQ=255)
                                                                               cATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGc                                                                                           <  7:44875/121‑1 (MQ=255)
                                                                               cATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGc                                                                                           >  8:44875/1‑121 (MQ=255)
                                                                                 tCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACgg                                                             >  3:423378/1‑149 (MQ=255)
                                                                                    gCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCg                                                                                 >  3:74258/1‑126 (MQ=255)
                                                                                    gCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCg                                                                                 <  4:74258/126‑1 (MQ=255)
                                                                                    gCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCa                                                          <  2:236900/149‑1 (MQ=255)
                                                                                          aTCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCtt                                                                            <  7:236951/125‑1 (MQ=255)
                                                                                          aTCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCtt                                                                            >  8:236951/1‑125 (MQ=255)
                                                                                           tCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATc                                                     >  2:278889/1‑147 (MQ=255)
                                                                                             aGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAg                                                  <  7:116576/148‑1 (MQ=255)
                                                                                             aGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAg                                                  <  5:303821/148‑1 (MQ=255)
                                                                                              gTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAg                                                  <  1:205465/147‑1 (MQ=255)
                                                                                              gTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAg                                                  >  2:205465/1‑147 (MQ=255)
                                                                                              gTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGc                                                 >  7:247289/1‑148 (MQ=255)
                                                                                                       acctacctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAg                                         <  7:52129/149‑1 (MQ=255)
                                                                                                        cctacctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAgg                                        <  4:159538/149‑1 (MQ=255)
                                                                                                         ctacctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGc                                       >  6:195757/1‑149 (MQ=255)
                                                                                                             ctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGc                                   <  5:195757/149‑1 (MQ=255)
                                                                                                               gCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCg                                 >  7:149886/1‑149 (MQ=255)
                                                                                                                ccAATCTGGGCCGCAGTACTGGTGGTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGGCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCg                                 <  8:2601/148‑1 (MQ=255)
                                                                                                                ccAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCg                                 <  4:142072/148‑1 (MQ=255)
                                                                                                                 cAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGt                               >  2:271133/1‑149 (MQ=255)
                                                                                                                  aaTCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGt                               >  7:382120/1‑148 (MQ=255)
                                                                                                                   aTCTGGGCCGCAGTACTGGTGCTGTTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGTCGTACCCGTCGctc                                                                              <  4:20968/100‑1 (MQ=255)
                                                                                                                   aTCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGctc                                                                              >  3:20968/1‑100 (MQ=255)
                                                                                                                        ggCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCAcc                           <  3:5446/146‑1 (MQ=255)
                                                                                                                        ggCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCg                          <  7:166003/147‑1 (MQ=255)
                                                                                                                        ggCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGa                         <  4:249595/148‑1 (MQ=255)
                                                                                                                            gCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATc                    >  2:36135/1‑149 (MQ=255)
                                                                                                                              aGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCgg                  <  3:362427/149‑1 (MQ=255)
                                                                                                                                   tggtgctggtgcTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGc                                                        <  3:195290/106‑1 (MQ=255)
                                                                                                                                   tggtgctggtgcTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGc                                                        >  4:195290/1‑106 (MQ=255)
                                                                                                                                       gcaggtgctTGAAGCCCTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGAcc                       >  2:186898/4‑135 (MQ=255)
                                                                                                                                       gcaggtgctTGAAGCCCTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGAcc                       <  1:186898/132‑1 (MQ=255)
                                                                                                                                       gcaggtgctTGAAGCCCTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTc         >  2:326058/4‑149 (MQ=255)
                                                                                                                                          ggtgctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAg      >  5:29936/1‑149 (MQ=255)
                                                                                                                                          ggtgctTGAAGCCCTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAg      >  2:286823/1‑149 (MQ=255)
                                                                                                                                          ggtgctTGAAGCACTCAAGCACCTGCTACGTCAGCGTCGCCGTTAGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCTACCAGCGGCAGGCGGTCCAg      >  8:48698/1‑149 (MQ=255)
                                                                                                                                           gtgctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCa       <  1:278889/147‑1 (MQ=255)
                                                                                                                                           gtgctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAg      <  5:152765/148‑1 (MQ=255)
                                                                                                                                             gctTGAGGGCTTCAAGCACCTGCTACCGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGc     >  1:225663/1‑148 (MQ=255)
                                                                                                                                              ctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGc     >  7:273205/1‑146 (MQ=255)
                                                                                                                                              ctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGGCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAg  <  5:30070/149‑1 (MQ=255)
                                                                                                                                               ttGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCg                                 >  7:331948/1‑117 (MQ=255)
                                                                                                                                               ttGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCg                                 <  8:331948/117‑1 (MQ=255)
                                                                                                                                                  aGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAg      >  2:7177/1‑141 (MQ=255)
                                                                                                                                                  aGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAg      <  1:7177/141‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TGCCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAG  >  NC_002947/3893337‑3893625

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: