Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 973,245 G→C 8.9% A513P (GCC→CCC)  secD → protein translocase subunit SecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947973,2450GC 8.9% 269.7 / 11.9 89A513P (GCC→CCC) secDprotein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base G (41/40);  new base C (4/4);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

CATGCTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCC  >  NC_002947/973106‑973391
                                                                                                                                            |                                                                                                                                                   
cATGCTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCagg                                                                                                                                             >  1:184280/1‑146 (MQ=255)
cATGCTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCagg                                                                                                                                             <  2:184280/148‑3 (MQ=255)
 aTGCTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCgt                                                                                                                                           <  5:93650/149‑1 (MQ=255)
    cTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGtt                                                                                                                                       <  4:124588/149‑1 (MQ=255)
     tGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGtt                                                                                                                                       >  8:34144/1‑148 (MQ=255)
     tGTTCGTCTCGCTGTTCATCATCGCCATCTACAGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGTGGCGAACCTGACCCTGCAGGGAATTGCCGGCATCGTGtt                                                                                                                                       >  6:282685/1‑148 (MQ=255)
       ttCGTCTCGCTGTTCATCATCGCCATCTACCGCGCGTTCGGTGTGCTGGACACCATCGGCCTGGCGGGGAACTTGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGAc                                                                                                                                    <  6:245233/149‑1 (MQ=255)
       ttCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGAc                                                                                                                                    <  7:316621/149‑1 (MQ=255)
          gTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCAt                                                                                                                                 <  4:30772/149‑1 (MQ=255)
             tcGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTAGACCCTGCCGGGTATTGCCGGAATCGTTTTTACCATgg                                                                                                                               >  1:127106/1‑149 (MQ=255)
             tcGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATggg                                                                                                                              >  3:20168/1‑149 (MQ=255)
             tcGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATggg                                                                                                                              >  4:27121/1‑149 (MQ=255)
              cGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCag                                                                                                                                              >  8:292341/1‑132 (MQ=255)
              cGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCag                                                                                                                                              <  7:292341/133‑2 (MQ=255)
               gCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTGGCTGCTGGGGCCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATgg                                                                                                                               <  1:266362/146‑1 (MQ=255)
                   ttCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGc                                                                                                                        >  5:147754/1‑149 (MQ=255)
                       tcatcGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTg                                                                                                                    <  3:27121/149‑1 (MQ=255)
                       tcatcGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTg                                                                                                                    >  2:180076/1‑149 (MQ=255)
                         atcGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCag                                                                                                                                              <  3:152488/122‑2 (MQ=255)
                         atcGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCag                                                                                                                                              >  4:152488/1‑121 (MQ=255)
                          tcGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACACTGACCCTGCCGGGTATTGACGGCATCGTGTTGACCATGGGTATGGCGGTGGAc                                                                                                                 >  5:347683/1‑149 (MQ=255)
                             ccATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGc                                                                                                               >  6:147547/1‑148 (MQ=255)
                                   aCCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTg                                                                                                        >  5:359339/1‑149 (MQ=255)
                                    ccGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTg                                                                                                        <  1:285258/148‑1 (MQ=255)
                                    ccGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGGGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGc                                                                                                       >  7:150988/1‑149 (MQ=255)
                                        gcCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGAt                                                                                                   <  1:180076/149‑1 (MQ=255)
                                              ggTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGGGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTtctc                                                                                             <  8:150988/149‑1 (MQ=255)
                                               gtgtGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTcg                                                                                            >  2:313259/1‑149 (MQ=255)
                                                     tGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGa                                                                                                    >  5:190354/1‑135 (MQ=255)
                                                     tGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGa                                                                                                    <  6:190354/135‑1 (MQ=255)
                                                     tGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTAtt                                                                                      >  5:329695/1‑149 (MQ=255)
                                                     tGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTAtt                                                                                      >  2:56417/1‑149 (MQ=255)
                                                     tGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTAtt                                                                                      >  4:173779/1‑149 (MQ=255)
                                                        ccaccaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcgc                                                                                   <  6:147754/149‑1 (MQ=255)
                                                        ccaccaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcgc                                                                                   >  5:65984/1‑149 (MQ=255)
                                                          accaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcgcg                                                                                  >  6:53416/1‑148 (MQ=255)
                                                          accaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGa                                                                                 <  7:323830/149‑1 (MQ=255)
                                                           ccaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGaa                                                                                <  8:286348/149‑1 (MQ=255)
                                                           ccaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGaa                                                                                <  2:93917/149‑1 (MQ=255)
                                                             aTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAaga                                                                              <  3:20707/149‑1 (MQ=255)
                                                             aTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAaga                                                                              <  3:175198/149‑1 (MQ=255)
                                                                        gggCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTg                                                                   <  6:65984/149‑1 (MQ=255)
                                                                        gggCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCagg                                                                                                                                             <  1:74849/76‑3 (MQ=255)
                                                                        gggCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCagg                                                                                                                                             >  2:74849/1‑74 (MQ=255)
                                                                          gCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGt                                                                 <  3:25058/149‑1 (MQ=255)
                                                                          gCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGt                                                                 <  1:56417/149‑1 (MQ=255)
                                                                             aCATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATg                                                              >  3:148812/1‑149 (MQ=255)
                                                                              cATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAaga                                                                              <  2:132345/132‑1 (MQ=255)
                                                                              cATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGt                                                             <  6:329695/149‑1 (MQ=255)
                                                                                   tgctgctgcTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTg                                                        <  2:127106/149‑1 (MQ=255)
                                                                                        ctgctgGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAgc                                                    <  5:53416/148‑1 (MQ=255)
                                                                                         tgctgGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCTACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAgcgc                                                  >  8:189296/1‑149 (MQ=255)
                                                                                            tgGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCg                                                                                                                                                  <  5:41057/51‑1 (MQ=255)
                                                                                            tgGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCg                                                                                                                                                  >  6:41057/1‑51 (MQ=255)
                                                                                             ggCGTTGATGTAGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGcc                                               >  6:60705/1‑148 (MQ=255)
                                                                                                gTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATc                                            >  4:215868/1‑148 (MQ=255)
                                                                                                     tGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGaa                                      <  7:34144/149‑1 (MQ=255)
                                                                                                      gTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGaa                                      <  5:60705/148‑1 (MQ=255)
                                                                                                      gTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAg                                     >  3:194417/1‑149 (MQ=255)
                                                                                                       tCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAg                                     <  1:313259/148‑1 (MQ=255)
                                                                                                         gctgctGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCt                                  >  3:355483/1‑149 (MQ=255)
                                                                                                          ctgctgGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCtt                                 <  3:239746/149‑1 (MQ=255)
                                                                                                            gctgGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGTGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAgcac                                                  >  8:82714/1‑128 (MQ=255)
                                                                                                            gctgGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGTGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAgcac                                                  <  7:82714/130‑3 (MQ=255)
                                                                                                              tgGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATg                                                                                                                          >  5:302713/1‑56 (MQ=255)
                                                                                                              tgGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATg                                                                                                                          <  6:302713/56‑1 (MQ=255)
                                                                                                              tgGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGcc                                               <  3:79587/131‑1 (MQ=255)
                                                                                                              tgGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGcc                                               >  4:79587/1‑131 (MQ=255)
                                                                                                               gggggCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATc                            >  1:251359/1‑149 (MQ=255)
                                                                                                               gggggCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATc                            >  5:205961/1‑149 (MQ=255)
                                                                                                                   gCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTtctc                                                                                             >  5:59412/1‑80 (MQ=255)
                                                                                                                   gCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTtctc                                                                                             <  6:59412/80‑1 (MQ=255)
                                                                                                                   gCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCgcgc                        <  4:148812/149‑1 (MQ=255)
                                                                                                                    cGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGAtcttct                                                                                              <  1:341629/78‑1 (MQ=255)
                                                                                                                    cGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGAtcttct                                                                                              >  2:341629/1‑78 (MQ=255)
                                                                                                                      accctgaccctgCCGGGTATTGCCGGCATCGTGGTAGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGcc                       >  1:46104/1‑148 (MQ=255)
                                                                                                                        cctgaccctgCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAaga                                                                              <  4:335201/90‑1 (MQ=255)
                                                                                                                        cctgaccctgCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAaga                                                                              >  3:335201/1‑90 (MQ=255)
                                                                                                                        cctgaccctgCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACt                   >  4:362143/1‑149 (MQ=255)
                                                                                                                         ctgaccctgCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTc                  >  8:280488/1‑149 (MQ=255)
                                                                                                                            gaccctgCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCt                                                                    >  4:16620/1‑97 (MQ=255)
                                                                                                                            gaccctgCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCt                                                                    <  3:16620/97‑1 (MQ=255)
                                                                                                                               cctgCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTTTGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCaa             >  7:295508/1‑149 (MQ=255)
                                                                                                                                 tgCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAAtc           >  4:266825/1‑149 (MQ=255)
                                                                                                                                  gCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCg                                                           >  2:103155/1‑100 (MQ=255)
                                                                                                                                   ccGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAAtcg          <  6:359339/148‑1 (MQ=255)
                                                                                                                                   ccGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAAtc           <  5:282685/147‑1 (MQ=255)
                                                                                                                                      ggTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGa      >  6:16004/1‑149 (MQ=255)
                                                                                                                                          ttGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGcc  <  4:269216/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CATGCTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCC  >  NC_002947/973106‑973391

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: