Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 954,436 T→G 5.2% L93R (CTC→CGC)  cyoE → protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947954,4360TG 5.2% 276.7 / 4.4 91L93R (CTC→CGC) cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base T (48/36);  new base G (2/3);  total (52/39)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

TACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGCCTGTC  >  NC_002947/954293‑954583
                                                                                                                                               |                                                                                                                                                   
tACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCgggg                                                                                                                                                <  3:428353/149‑1 (MQ=255)
   ttCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTgg                                                                                                                                             >  4:259961/1‑149 (MQ=255)
    tCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGc                                                                                                                                            >  5:436165/1‑149 (MQ=255)
      gCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCgg                                                                                                                                          <  1:275406/149‑1 (MQ=255)
      gCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCgg                                                                                                                                          <  1:339960/149‑1 (MQ=255)
      gCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCgg                                                                                                                                          <  3:259961/149‑1 (MQ=255)
      gCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCg                                                                                                                                           <  7:182449/148‑1 (MQ=255)
      gCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCg                                                                                                                                           <  4:137270/148‑1 (MQ=255)
      gCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCACCCCGAGCAGGGTGGCgt                                                                                                                                          <  4:399390/148‑2 (MQ=255)
       cTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCgg                                                                                                                                          <  2:318087/148‑1 (MQ=255)
         ggtggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGtt                                                                                                                                       >  4:267089/1‑149 (MQ=255)
          gtggtTGCGTCCGGATGCGTGTTCAACAACTGCAGCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGt                                                                                                                                        <  4:256846/147‑1 (MQ=255)
           tggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTc                                                                                                                                     <  2:13121/149‑1 (MQ=255)
            ggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGTTCTACGCCACCCTGCTCGGGGTGGCGGGttt                                                                                                                                      >  5:370129/1‑147 (MQ=255)
            ggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGttt                                                                                                                                      >  6:332166/1‑147 (MQ=255)
              ttGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGc                                                                                                                                  >  1:103057/1‑149 (MQ=255)
                  gTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCctgc                                                                                                                              >  7:21134/1‑149 (MQ=255)
                  gTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCctgc                                                                                                                              <  7:273183/149‑1 (MQ=255)
                   tCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCATCCCGAGCAGGGTGGCg                                                                                                                                           <  6:423852/135‑1 (MQ=255)
                   tCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCATCCCGAGCAGGGTGGCg                                                                                                                                           >  5:423852/1‑135 (MQ=255)
                     cGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGt                                                                                                                           <  1:244739/149‑1 (MQ=255)
                     cGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGt                                                                                                                           <  3:267089/149‑1 (MQ=255)
                      ggATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTa                                                                                                                          >  2:249929/1‑149 (MQ=255)
                             tgtTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAg                                                                                                                   <  2:103057/149‑1 (MQ=255)
                                 caacaaCTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGcc                                                                                                                >  7:171577/1‑148 (MQ=255)
                                  aacaacTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCaa                                                                                                              >  7:257078/1‑149 (MQ=255)
                                  aacaacTGCATCGACCGTGACATAGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCaa                                                                                                              >  5:264075/1‑149 (MQ=255)
                                   acaacTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAAc                                                                                                             >  2:307971/1‑149 (MQ=255)
                                       cTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGc                                                                                                         >  3:186941/1‑149 (MQ=255)
                                          caacGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTCTCATGCTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATATACGCCACCCTGCAAGGGGTGGCGGGTATCAGCCTGCTGTATGTCCATCCCAACCCGCTg                                                                                                       >  4:306587/4‑148 (MQ=255)
                                          cATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGt                                                                                                      <  6:292390/149‑1 (MQ=255)
                                           aTCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAgg                                                                                                                  >  6:190855/1‑136 (MQ=255)
                                           aTCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAgg                                                                                                                  <  5:190855/136‑1 (MQ=255)
                                               aCCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCg                                                                                                 >  7:154599/1‑149 (MQ=255)
                                                  gTGACATCGACCACAAGATGGAGCTCACGAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTc                                                                                              <  7:251248/149‑1 (MQ=255)
                                                  gTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCa                                                                                                                                    >  6:21263/1‑111 (MQ=255)
                                                  gTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCa                                                                                                                                    <  5:21263/111‑1 (MQ=255)
                                                   tGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCg                                                                                                 <  5:78644/145‑1 (MQ=255)
                                                   tGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCg                                                                                                 >  6:78644/1‑145 (MQ=255)
                                                   tGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCt                                                                                             >  2:21847/1‑149 (MQ=255)
                                                   tGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCACGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCt                                                                                             >  2:403476/1‑149 (MQ=255)
                                                      cATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTgcg                                                                                          >  3:1699/1‑149 (MQ=255)
                                                         cGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGTATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCAc                                                                                       >  2:430717/1‑149 (MQ=255)
                                                            ccACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCa                                                                                        >  7:241444/1‑145 (MQ=255)
                                                            ccACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCa                                                                                        <  8:241444/145‑1 (MQ=255)
                                                            ccACAAGATGGAGCGCACCAAGAACCGTGTCATCGTTCAGGGCGGCATGTCGCTGCCCCTCCCGCTCATCTACGCCACCCTGCTCGGGGTGCCGGGTTTCAGCCTGCTCTATGACCAGGCCAACCCGCTGACAGCGTTCTGCGCACTg                                                                                     >  8:295796/1‑148 (MQ=255)
                                                             cacaAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGAt                                                                                   >  4:369226/1‑149 (MQ=255)
                                                             cacaAGATGGAGCGCACCAACAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGAt                                                                                   >  4:163957/1‑149 (MQ=255)
                                                              acaAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGAtt                                                                                  <  5:332166/149‑1 (MQ=255)
                                                              acaAGATGGAGCGCACAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTg                                                                                 >  6:235169/1‑149 (MQ=255)
                                                              acaAGATGGAGCGCACAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTg                                                                                 >  6:229241/1‑149 (MQ=255)
                                                               caAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTg                                                                                 >  7:82409/1‑149 (MQ=255)
                                                               caAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTg                                                                                 >  6:393334/1‑149 (MQ=255)
                                                                aaGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTgg                                                                                <  6:370129/149‑1 (MQ=255)
                                                                      gAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCAGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTgt                                                                          >  5:29162/1‑149 (MQ=255)
                                                                          gcACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTgtcgtc                                                                      >  6:372205/1‑149 (MQ=255)
                                                                            aCCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTa                                                                    <  3:211092/149‑1 (MQ=255)
                                                                             ccAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTAc                                                                   <  4:186941/149‑1 (MQ=255)
                                                                              cAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACg                                                                  >  1:45784/1‑149 (MQ=255)
                                                                              cAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACg                                                                  >  4:376938/1‑149 (MQ=255)
                                                                                 gAACCGTCTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGATGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCGGCGCACTGATTGGAATTGTCGTCTacgacg                                                               >  8:339727/1‑149 (MQ=255)
                                                                                         tCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTa                                                        >  1:56224/1‑148 (MQ=255)
                                                                                          cATGGTTCAGGGCGGCAAGTCGCTGCCCCTCGCGCTGATCTACGCCACCATGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCCACCCGCTGTCTGCGGTCTGCGCACTGATTGGCTTTGTCGTGTACGTCGGTTTCTa                                                        >  6:305426/1‑147 (MQ=255)
                                                                                              gTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCt                                                  >  4:375499/1‑149 (MQ=255)
                                                                                              gTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCt                                                  >  3:61672/1‑149 (MQ=255)
                                                                                              gTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCt                                                  >  2:51749/1‑149 (MQ=255)
                                                                                                tCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTgtc                                                                         >  1:184983/1‑124 (MQ=255)
                                                                                                tCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTgtc                                                                         <  2:184983/124‑1 (MQ=255)
                                                                                                  aGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCtgtg                                               <  6:436165/148‑1 (MQ=255)
                                                                                                  aGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTgg                                              >  2:91213/1‑149 (MQ=255)
                                                                                                        gcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAg                                        <  1:21847/149‑1 (MQ=255)
                                                                                                         cATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGc                                       >  5:152055/1‑149 (MQ=255)
                                                                                                            gTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGt                                     >  5:331563/1‑148 (MQ=255)
                                                                                                               gCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAAt                                 >  6:444627/1‑149 (MQ=255)
                                                                                                                   ccccTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGt                             >  1:291670/1‑149 (MQ=255)
                                                                                                                   ccccTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTATTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCg                               <  1:307971/147‑1 (MQ=255)
                                                                                                                    cccTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTg                            <  1:403476/149‑1 (MQ=255)
                                                                                                                    cccTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCAGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGt                             >  2:254848/1‑148 (MQ=255)
                                                                                                                      cTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCa                          >  6:195391/1‑149 (MQ=255)
                                                                                                                         gcgcTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCAc                         <  8:154599/147‑1 (MQ=255)
                                                                                                                         gcgcTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACgg                       <  8:171577/149‑1 (MQ=255)
                                                                                                                               aTCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCt                  <  8:82409/148‑1 (MQ=255)
                                                                                                                                  cccGCCACCCCGAGCAGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGt              >  3:210310/3‑149 (MQ=255)
                                                                                                                                   aCGCCACACTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGAGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTc             >  8:249607/1‑149 (MQ=255)
                                                                                                                                    cGCCACCCTGCTCGGGGTGGCTGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTgg               <  3:163957/146‑1 (MQ=255)
                                                                                                                                    cGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTAGTc             <  4:1699/148‑1 (MQ=255)
                                                                                                                                      ccACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCgg           >  5:93148/1‑148 (MQ=255)
                                                                                                                                        aCCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAg        <  8:21134/149‑1 (MQ=255)
                                                                                                                                         cccTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGc       <  5:229241/149‑1 (MQ=255)
                                                                                                                                         cccTGCGCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGc       <  5:235169/149‑1 (MQ=255)
                                                                                                                                              cTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGCCTGTc  <  5:372205/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGCCTGTC  >  NC_002947/954293‑954583

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: