Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 973,244 T→A 8.7% I512I (ATT→ATA secD → protein translocase subunit SecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947973,2440TA 8.7% 256.4 / 9.8 88I512I (ATT→ATAsecDprotein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base T (41/39);  new base A (4/4);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGTGGGGCATGCTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCC  >  NC_002947/973097‑973391
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                              >  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                                                                                      >  4:173779/1‑149 (MQ=255)
                                                              tGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTAtt                                                                                      >  5:329695/1‑149 (MQ=255)
                                                              tGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTAtt                                                                                      >  2:56417/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                                                                 <  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                                                  <  7:82714/130‑3 (MQ=255)
                                                                                                                     gctgGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGTGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAgcac                                                  >  8:82714/1‑128 (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)
                                                                                                                                                    |                                                                                                                                                    
GCTGTGGGGCATGCTGTTCGTCTCGCTGTTCATCATCGCCATCTACCGCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCC  >  NC_002947/973097‑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: