Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 604,827 C→T 6.4% T46I (ACC→ATC) ‡ thiL → thiamin monophosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947604,8270CT 6.4% 265.5 / 4.0 89T46I (ACC→ATC) ‡thiLthiamin monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (40/43);  new base T (4/2);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTGCCTGGCCATGGGTGAGTTCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGC  >  NC_002947/604678‑604967
                                                                                                                                                     |                                                                                                                                            
gcgcTGCCTGGCCATGGGTGAGTTCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTcgac                                                                                                                                              <  1:961980/150‑1 (MQ=255)
     gCCTGGCCATGGGTGAGTTCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGaca                                                                                                                                         <  1:2181339/150‑1 (MQ=255)
       cTGGCCATGGGTGAGTTCGAGCTGATCAACCATTACTTCGCCGCAGCGTGCAGTGGGCAGGGCGGCGAGGCCTTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGGTCGACAcc                                                                                                                                       <  2:573079/150‑1 (MQ=255)
        tGGCCATGGGTGAGTTCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACAccc                                                                                                                                      <  1:1765991/150‑1 (MQ=255)
          gCCATGGGTGAGTTCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTg                                                                                                                                    >  1:1297770/1‑150 (MQ=255)
                 gTGAGTTCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCg                                                                                                                             >  1:1345331/1‑150 (MQ=255)
                  tGAGTTCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTc                                                                                                                                 <  1:407648/145‑1 (MQ=255)
                  tGAGTTCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTc                                                                                                                                 >  2:407648/1‑145 (MQ=255)
                  tGAGTTCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCgg                                                                                                                            >  2:1798479/1‑150 (MQ=255)
                     gTTCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGTACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGt                                                                                                                         <  1:540211/150‑1 (MQ=255)
                     gTTCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGt                                                                                                                         <  1:1260084/150‑1 (MQ=255)
                       tCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCAGGTACGACCGACACCCTGGTCGCCGGGGTGc                                                                                                                       <  1:1268845/150‑1 (MQ=255)
                         gAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCAc                                                                                                                     >  2:43680/1‑150 (MQ=255)
                             tGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTc                                                                                                                 <  2:55601/150‑1 (MQ=255)
                             tGATCAACCATTACTTCGCCGCAGCGCACTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTc                                                                                                                 <  2:102828/150‑1 (MQ=255)
                              gATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCt                                                                                                                                     >  1:1177134/1‑129 (MQ=255)
                              gATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCt                                                                                                                                     <  2:1177134/129‑1 (MQ=255)
                                tCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCctg                                                                                                              <  1:2258583/150‑1 (MQ=255)
                                  aaCCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCctgct                                                                                                            <  1:837950/150‑1 (MQ=255)
                                  aaCCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCctgct                                                                                                            <  2:1345331/150‑1 (MQ=255)
                                  aaCCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCctgct                                                                                                            <  2:1447518/150‑1 (MQ=255)
                                    ccATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGt                                                                                                          >  1:2030817/1‑150 (MQ=255)
                                     cATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGtt                                                                                                         >  1:1136050/1‑150 (MQ=255)
                                           ttCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGAt                                                                                                   >  1:2053999/1‑150 (MQ=255)
                                           ttCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGAt                                                                                                   >  2:2395534/1‑150 (MQ=255)
                                               ccgcAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGt                                                                                               >  1:2146034/1‑150 (MQ=255)
                                                    ccgcCCTGCGCACAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTtgctg                                                                                          >  1:85148/2‑150 (MQ=255)
                                                        ccTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCttt                                                                                      <  2:1297770/150‑2 (MQ=255)
                                                         cTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGtt                                                                                              >  2:2403214/1‑141 (MQ=255)
                                                         cTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGtt                                                                                              <  1:2403214/141‑1 (MQ=255)
                                                         cTGTGCGCAGGACGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCACCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTgg                                                                                     >  1:2370684/1‑150 (MQ=255)
                                                          tgtgCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGc                                                                                    >  2:1135071/1‑150 (MQ=255)
                                                            tgCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCa                                                                                  <  2:790915/150‑1 (MQ=255)
                                                            tgCGCAGGGCGGCGAGGGCATGGCCCTGGCTATCAGCTACTACTGCGCCCTGCTGGACCTTCCCGGCGGCCGGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCa                                                                                  <  1:43680/150‑1 (MQ=255)
                                                               gcAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGcg                                                                               >  2:456437/1‑150 (MQ=255)
                                                                 aGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACCGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGACGGTGTCGACCGACCCCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGCTTCGTTGCTGCTTTGACAGcact                                                                             >  2:383714/1‑147 (MQ=255)
                                                                  gggCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGcgttc                                                                            <  1:456437/150‑1 (MQ=255)
                                                                  gggCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGcgttc                                                                            <  2:2146034/150‑1 (MQ=255)
                                                                  gggCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGcgttc                                                                            <  2:1136050/150‑1 (MQ=255)
                                                                    gcggcgAGGGCGTGTCGCTGGGTACCGGCGACGACTGCGCCCTTCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCAACACCTTGCTCGCCGGGGTGCCCTGGCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGcgttcgt                                                                          <  2:85148/150‑1 (MQ=255)
                                                                      ggcgAGGGCGTGGCTCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACAccc                                                                                                                                      <  1:163288/88‑1 (MQ=255)
                                                                      ggcgAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACAccc                                                                                                                                      >  2:163288/1‑88 (MQ=255)
                                                                            ggCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCtgtg                                                                  >  1:668302/1‑150 (MQ=255)
                                                                                tGGCGCTGGGTATCGGCGACGACTGCGCGCTGCTGGACCTTCCCCCCGCTGAGCAACTTGCTGTGTCGACCGACACCCTGGTCGCCGGCGTGCCCTTTCCTGCTGTTTGGGATCCGTTGCTTCTTGGCCAGCGTTCGTTGCCttcttcag                                                              <  2:1296930/150‑8 (MQ=255)
                                                                                        ggTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGAc                                                      <  2:1094094/150‑1 (MQ=255)
                                                                                         gTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTg                                                                                                                                    >  2:2164432/1‑71 (MQ=255)
                                                                                         gTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTg                                                                                                                                    <  1:2164432/71‑1 (MQ=255)
                                                                                            tCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGtt                                                                                              >  2:2048688/1‑106 (MQ=255)
                                                                                            tCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGtt                                                                                              <  1:2048688/106‑1 (MQ=255)
                                                                                            tCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTgg                                                  <  1:1798479/150‑1 (MQ=255)
                                                                                             cGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGc                                                 >  2:2282106/1‑150 (MQ=255)
                                                                                              ggCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTtg                                                                                             <  1:887877/105‑1 (MQ=255)
                                                                                              ggCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTtg                                                                                             >  2:887877/1‑105 (MQ=255)
                                                                                               gCGACGACTGCGCCCTGCTGGACCTTGGATCGGGGGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGTGTGCACTTTCCTGCTGTTTGCGATCCGTTGCGGCTTGGCCAGCGTTCGTTGGCTGGAGCGGCCAGTGACCTGGCAg                                               >  1:675986/1‑150 (MQ=255)
                                                                                                   acaacgcgcCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCAt                                           <  1:2394355/145‑1 (MQ=255)
                                                                                                    gacTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGc                                                                                                      <  1:1982589/90‑1 (MQ=255)
                                                                                                    gacTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGc                                                                                                      >  2:1982589/1‑90 (MQ=255)
                                                                                                     acTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATgg                                         <  1:383714/150‑1 (MQ=255)
                                                                                                       tGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTg                                                                                                                                    >  1:147302/1‑57 (MQ=255)
                                                                                                       tGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTg                                                                                                                                    <  2:147302/57‑1 (MQ=255)
                                                                                                       tGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGgc                                       >  1:2284138/1‑150 (MQ=255)
                                                                                                       tGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGgc                                       >  2:1748167/1‑150 (MQ=255)
                                                                                                            ccTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGgc                                       <  1:686084/145‑1 (MQ=255)
                                                                                                            ccTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGgc                                       >  2:686084/1‑145 (MQ=255)
                                                                                                            ccTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGAc                                  >  2:2345695/1‑150 (MQ=255)
                                                                                                            ccTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGAc                                  >  1:1206108/1‑150 (MQ=255)
                                                                                                                  ggACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAAt                            >  1:177563/1‑150 (MQ=255)
                                                                                                                     ccTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCgg                         >  1:1535838/1‑150 (MQ=255)
                                                                                                                       ttCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTTGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGGCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGAGCGACTGCAATCGGCt                       >  1:627914/1‑150 (MQ=255)
                                                                                                                           ccccGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCAc                   >  1:1047705/1‑150 (MQ=255)
                                                                                                                            cccGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTg                                                                                      <  2:1382220/82‑1 (MQ=255)
                                                                                                                            cccGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTg                                                                                      >  1:1382220/1‑82 (MQ=255)
                                                                                                                             ccGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCAccc                 <  1:1135071/150‑1 (MQ=255)
                                                                                                                             ccGGTGAGCAACTGGCGGTGTCGACCGACACACTGGTCGCCGGGGTATCGCGGGCTGGTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGTGACTGCAATCGGCTTCAccc                 <  2:675986/150‑1 (MQ=255)
                                                                                                                               ggTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTc               <  2:668302/150‑1 (MQ=255)
                                                                                                                                   aGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCAc           <  2:1535838/150‑1 (MQ=255)
                                                                                                                                    gCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTggc                                                                >  2:951333/1‑96 (MQ=255)
                                                                                                                                    gCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTggc                                                                <  1:951333/96‑1 (MQ=255)
                                                                                                                                      aaCTGGCGGTGTCAACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGa        >  1:1284872/1‑150 (MQ=255)
                                                                                                                                       gaccagggTGTCGGTCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATTGTCGCTACTGTATTCGGCTTCCACCTCGCACTGAc       >  2:998780/7‑150 (MQ=255)
                                                                                                                                       aCTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGGGGCCAGTGCCCAAGCAGTCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGAc       <  2:1667586/150‑1 (MQ=255)
                                                                                                                                         tGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCtt                                                                                       >  2:1285003/1‑68 (MQ=255)
                                                                                                                                         tGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCtt                                                                                       <  1:1285003/68‑1 (MQ=255)
                                                                                                                                          cagggTGTCGGTCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTg                                                                                      >  2:1991896/4‑68 (MQ=38)
                                                                                                                                          cagggTGTCGGTCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTg                                                                                      <  1:1991896/65‑1 (MQ=255)
                                                                                                                                          cagggTGTCGGTCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGTTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACcct    >  1:752439/4‑150 (MQ=255)
                                                                                                                                          cagggTGTCGGTCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACcct    >  1:1292110/4‑150 (MQ=255)
                                                                                                                                            cGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACcctgc  <  1:2395534/150‑1 (MQ=255)
                                                                                                                                              gtgtCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCa                    >  1:2245862/1‑130 (MQ=255)
                                                                                                                                              gtgtCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCa                    <  2:2245862/130‑1 (MQ=255)
                                                                                                                                                   gaccgacACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGAc       <  1:156844/138‑1 (MQ=255)
                                                                                                                                                   gaccgacACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGAc       >  2:156844/1‑138 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GCGCTGCCTGGCCATGGGTGAGTTCGAGCTGATCAACCATTACTTCGCCGCAGCGCCCTGTGCGCAGGGCGGCGAGGGCGTGGCGCTGGGTATCGGCGACGACTGCGCCCTGCTGGACCTTCCCCCCGGTGAGCAACTGGCGGTGTCGACCGACACCCTGGTCGCCGGGGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGC  >  NC_002947/604678‑604967

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: