Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 973,240 G→C 8.8% G511A (GGT→GCT) ‡ secD → protein translocase subunit SecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947973,2400GC 8.8% 272.6 / 12.3 91G511A (GGT→GCT) ‡secDprotein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base G (42/40);  new base C (4/4);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TCGCTGTGGGGCATGCTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGA  >  NC_002947/973095‑973387
                                                                                                                                                  |                                                                                                                                                    
tCGCTGTGGGGCATGCTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTAt                                                                                                                                                   >  4:25058/1‑149 (MQ=255)
  gCTGTTGCGAATTCTGTTCGTCTCGCTGTTCATCATAGACATCTACCGCGCATTCGTTGTGCTGGCCACCATCGCCCTGGCGGTCAACATGGTGATTCTGCTGGCGTTGATGTCGCTGCTGTGGTCGACCCTTACCCTGCCGTTTAtt                                                                                                                                                  <  4:21547/148‑1 (MQ=255)
           cATGCTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCagg                                                                                                                                         <  2:184280/148‑3 (MQ=255)
           cATGCTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCagg                                                                                                                                         >  1:184280/1‑146 (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                                                                                                                          >  4:27121/1‑149 (MQ=255)
                        tcGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATggg                                                                                                                          >  3:20168/1‑149 (MQ=255)
                         cGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCag                                                                                                                                          <  7:292341/133‑2 (MQ=255)
                         cGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCag                                                                                                                                          >  8:292341/1‑132 (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)
                                    atcGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTa                                                                                                                                                    >  2:337863/1‑112 (MQ=255)
                                    atcGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTa                                                                                                                                                    <  1:337863/112‑1 (MQ=255)
                                    atcGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCag                                                                                                                                          <  3:152488/122‑2 (MQ=255)
                                    atcGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGCAATACCCGGCag                                                                                                                                          >  4:152488/1‑121 (MQ=255)
                                     tcGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACACTGACCCTGCCGGGTATTGACGGCATCGTGTTGACCATGGGTATGGCGGTGGAc                                                                                                             >  5:347683/1‑149 (MQ=255)
                                        ccATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCggg                                                                                                                                                      >  2:375121/1‑106 (MQ=255)
                                        ccATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCggg                                                                                                                                                      <  1:375121/106‑1 (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                                                                                  >  2:56417/1‑149 (MQ=255)
                                                                tGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTAtt                                                                                  >  5:329695/1‑149 (MQ=255)
                                                                tGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTAtt                                                                                  >  4:173779/1‑149 (MQ=255)
                                                                   ccaccaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcgc                                                                               >  5:65984/1‑149 (MQ=255)
                                                                   ccaccaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcgc                                                                               <  6:147754/149‑1 (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:175198/149‑1 (MQ=255)
                                                                        aTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAaga                                                                          <  3:20707/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                                                             <  1:56417/149‑1 (MQ=255)
                                                                                     gCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGt                                                             <  3:25058/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                                                                                                                                              >  6:41057/1‑51 (MQ=255)
                                                                                                       tgGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCg                                                                                                                                              <  5:41057/51‑1 (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                                           >  4:79587/1‑131 (MQ=255)
                                                                                                                         tgGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGcc                                           <  3:79587/131‑1 (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                                                                                          >  2:341629/1‑78 (MQ=255)
                                                                                                                               cGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGAtcttct                                                                                          <  1:341629/78‑1 (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)
                                                                                                                                                  |                                                                                                                                                    
TCGCTGTGGGGCATGCTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGA  >  NC_002947/973095‑973387

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: