Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,893,486 T→C 6.0% F284L (TTC→CTC)  rarD → chloramphenicol‑sensitive protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,893,4860TC 6.0% 197.8 / 2.5 67F284L (TTC→CTC) rarDchloramphenicol‑sensitive protein
Reads supporting (aligned to +/- strand):  ref base T (36/27);  new base C (2/2);  total (38/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 = 7.95e-01

CCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCTCATCAC  >  NC_002947/3893339‑3893633
                                                                                                                                                    |                                                                                                                                                    
ccATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCTGTGGCGTTGCTGCTGGGCGAAACAATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCt                                                                                                                                                      >  1:283398/1‑149 (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                                                                                                   <  7:28780/149‑1 (MQ=255)
                                                 tAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGc                                                                                                   <  8:193678/149‑1 (MQ=255)
                                                        tggcgtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTAccc                                                                                             >  7:32995/1‑149 (MQ=255)
                                                            gtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                             >  7:178399/1‑98 (MQ=255)
                                                            gtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                             <  8:178399/98‑1 (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                                                                                                                                             >  5:357688/1‑74 (MQ=255)
                                                                                    cccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                             <  6:357688/74‑1 (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                 <  3:24926/143‑1 (MQ=255)
                                                                                                                                          gctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTc                 >  4:24926/1‑143 (MQ=255)
                                                                                                                                           ctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAg          >  5:370167/1‑149 (MQ=255)
                                                                                                                                            ttGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAg          >  4:357285/1‑148 (MQ=255)
                                                                                                                                                 ggCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAg                                                 <  3:130896/104‑1 (MQ=255)
                                                                                                                                                 ggCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAg                                                 >  4:130896/1‑104 (MQ=255)
                                                                                                                                                 ggCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCtcatc    <  3:103469/149‑1 (MQ=255)
                                                                                                                                                    ttCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCTCATCAc  >  4:264143/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCTCATCAC  >  NC_002947/3893339‑3893633

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: