Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,893,484 G→C 5.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 5.9% 221.6 / 3.7 68G283A (GGC→GCC) rarDchloramphenicol‑sensitive protein
Reads supporting (aligned to +/- strand):  ref base G (37/27);  new base C (2/2);  total (39/29)
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.51e-01

TGCCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCTCATC  >  NC_002947/3893337‑3893631
                                                                                                                                                    |                                                                                                                                                    
tgcCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGc                                                                                                                                                     >  7:341652/1‑149 (MQ=255)
  ccATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCTGTGGCGTTGCTGCTGGGCGAAACAATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCt                                                                                                                                                    >  1:283398/1‑149 (MQ=255)
  ccATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGc                                                                                                                                                     <  5:373473/147‑1 (MQ=255)
  ccATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGc                                                                                                                                                     <  2:34058/147‑1 (MQ=255)
    aTTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTc                                                                                                                                                  <  5:91234/148‑1 (MQ=255)
    aTTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTc                                                                                                                                                  >  6:91234/1‑148 (MQ=255)
        aGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGc                                                                                                                                              <  7:125867/148‑1 (MQ=255)
         gCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAc                                                                                                                                            >  1:9877/1‑149 (MQ=255)
         gCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAAGCCCTCAAGCAc                                                                                                                                            <  1:371793/149‑1 (MQ=255)
          ccTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                           >  7:193678/1‑149 (MQ=255)
          ccTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAAGCCCTCAAGCAcc                                                                                                                                           <  4:213530/149‑1 (MQ=255)
           cTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCt                                                                                                                                          <  2:283398/149‑1 (MQ=255)
            tGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTg                                                                                                                                         >  1:236127/1‑149 (MQ=255)
              ttCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGc                                                                                                                                        >  5:163280/1‑148 (MQ=255)
                  gCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGt                                                                                                                                   >  7:206513/1‑149 (MQ=255)
                  gCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGt                                                                                                                                   <  7:253307/149‑1 (MQ=255)
                    cTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCa                                                                                                                                                 >  2:271638/1‑133 (MQ=255)
                    cTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCa                                                                                                                                                 <  1:271638/133‑1 (MQ=255)
                     tGCTCAGCTACGTCGAACCCGTGCTTCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACATACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTAAGTccg                                                                                                                                >  2:260219/1‑147 (MQ=255)
                       cTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCg                                                                                                                              >  8:28780/1‑149 (MQ=255)
                                 tCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCg                                                                                                                   <  2:9877/149‑1 (MQ=255)
                                 tCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCg                                                                                                                   <  1:348786/149‑1 (MQ=255)
                                               ctgGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGCTTACACCGAAcg                                                                                                     <  2:236127/149‑2 (MQ=255)
                                                   tAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCGGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGc                                                                                                 <  4:42858/149‑1 (MQ=255)
                                                   tAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGc                                                                                                 <  8:193678/149‑1 (MQ=255)
                                                   tAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGc                                                                                                 <  7:28780/149‑1 (MQ=255)
                                                          tggcgtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTAccc                                                                                           >  7:32995/1‑149 (MQ=255)
                                                              gtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                           <  8:178399/98‑1 (MQ=255)
                                                              gtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                           >  7:178399/1‑98 (MQ=255)
                                                              gtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCg                                                                                       >  4:345507/1‑149 (MQ=255)
                                                              gtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCg                                                                                       >  2:336812/1‑149 (MQ=255)
                                                                tgctgcTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGct                                                                                     >  6:252798/1‑149 (MQ=255)
                                                                      tgGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCaa                                                                               >  2:329074/1‑149 (MQ=255)
                                                                        ggCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAAc                                                                              <  8:32995/148‑1 (MQ=255)
                                                                          cGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGc                                                                           >  1:173480/1‑149 (MQ=255)
                                                                           gAAACCATCGGCCCGGATCAGTGGCAGACCTACCTGCCAAACTGGGCCGCAGTACTGGAGCTGGTGATTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGcc                                                                          >  8:239591/1‑149 (MQ=255)
                                                                               ccATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCa                                                                      >  8:288801/1‑149 (MQ=255)
                                                                                 aTCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACg                                                                    >  7:15340/1‑149 (MQ=255)
                                                                                     gCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCa                                                                <  6:163280/149‑1 (MQ=255)
                                                                                     gCCCGGATCAGAGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCa                                                                <  8:341652/149‑1 (MQ=255)
                                                                                      cccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                           <  6:357688/74‑1 (MQ=255)
                                                                                      cccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                           >  5:357688/1‑74 (MQ=255)
                                                                                      cccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTgg                                                                                                  >  6:212446/1‑114 (MQ=255)
                                                                                      cccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTgg                                                                                                  <  5:212446/114‑1 (MQ=255)
                                                                                      cccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAg                                                               >  7:229813/1‑149 (MQ=255)
                                                                                          gATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATc                                                           <  2:173480/149‑1 (MQ=255)
                                                                                            tCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCa                                                         <  5:252798/149‑1 (MQ=255)
                                                                                                      acctacctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTcc                                                 <  8:15340/147‑1 (MQ=255)
                                                                                                        ctacctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGc                                             >  7:351579/1‑149 (MQ=255)
                                                                                                           cctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCCTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAg                                          >  4:90711/1‑149 (MQ=255)
                                                                                                             tGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGcc                                        >  4:103469/1‑149 (MQ=255)
                                                                                                              gCCAATCTGGGCCGCAGTACTGGTACTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTTTAAGCCGAACCTGGCGTACCCGTCGATATTCAACCGCCTCCACGGTCAGAAACAAGCCATACCAGGCAAGCCg                                       >  7:235557/1‑149 (MQ=255)
                                                                                                               ccAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCg                                       <  1:336812/148‑1 (MQ=255)
                                                                                                                cAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGt                                     >  8:319242/1‑149 (MQ=255)
                                                                                                                     tGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCAcc                                 >  6:104756/1‑148 (MQ=255)
                                                                                                                      gggCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGa                               >  2:122884/1‑149 (MQ=255)
                                                                                                                            cAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCg                         >  4:138780/1‑149 (MQ=255)
                                                                                                                                 ctggtgctggtgCTTGAGGGCTTCAAGCACCTGCTACGTGAGCGTCGCCGTTCGGTGTTAGCCGTACCTGGCGTACGCGTCGCTGTTCAACCTCCTCCACGGTGAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAg                     <  7:239591/148‑1 (MQ=255)
                                                                                                                                        tggtgctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCAACGACCATCGGCAGGCGGTCCa             >  7:58326/1‑149 (MQ=255)
                                                                                                                                        aggtgctTGAAGCCCTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTc               >  7:102770/2‑147 (MQ=255)
                                                                                                                                        aggtgctTGAAGCCCTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCACCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTc               >  6:51066/2‑147 (MQ=255)
                                                                                                                                          gtgctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGc           <  8:206513/149‑1 (MQ=255)
                                                                                                                                           tgctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCa          <  8:229813/149‑1 (MQ=255)
                                                                                                                                            gctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCAACCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGc           >  1:258558/1‑148 (MQ=255)
                                                                                                                                            gctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTc               >  4:24926/1‑143 (MQ=255)
                                                                                                                                            gctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTc               <  3:24926/143‑1 (MQ=255)
                                                                                                                                             ctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAg        >  5:370167/1‑149 (MQ=255)
                                                                                                                                              ttGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAg        >  4:357285/1‑148 (MQ=255)
                                                                                                                                                   ggCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAg                                               >  4:130896/1‑104 (MQ=255)
                                                                                                                                                   ggCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAg                                               <  3:130896/104‑1 (MQ=255)
                                                                                                                                                   ggCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCtcatc  <  3:103469/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCTCATC  >  NC_002947/3893337‑3893631

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: