Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 296,372 T→G 5.4% E19A (GAG→GCG) ‡ PP_0242 ← TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947296,3720TG 5.4% 273.1 / 6.0 92E19A (GAG→GCG) ‡PP_0242TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (40/47);  new base G (2/3);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

CAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGA  >  NC_002947/296223‑296521
                                                                                                                                                       |                                                                                                                                                     
cAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCt                                                                                                                                                       <  2:168999/150‑1 (MQ=255)
cAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCt                                                                                                                                                       <  2:584917/150‑1 (MQ=255)
         ccTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCAACCCTTGTTCGGCCAACAGCTCACGCGCAGAC‑‑GCCAGGATCGCCTCGATcctgg                                                                                                                                              >  1:2299054/1‑150 (MQ=255)
          cTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATcctggc                                                                                                                                             >  1:2282087/1‑150 (MQ=255)
          cTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGAGGCGATcctggc                                                                                                                                             <  1:1342979/150‑1 (MQ=255)
           tCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAAGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCGAAGAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATcctggcc                                                                                                                                            >  1:633384/1‑150 (MQ=255)
              aGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATcctggcctgg                                                                                                                                         >  2:1472030/1‑150 (MQ=255)
              aGCAAGGCCGGCACGCTGACGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATcctggcctgg                                                                                                                                         <  1:2323320/150‑1 (MQ=255)
                   ggCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCCAGGCCAGGATCGAGGCGATcctggcggctgcg                                                                                                                                      <  2:99505/150‑8 (MQ=38)
                      cGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGgc                                                                                                                                 <  2:184876/150‑1 (MQ=255)
                       ggCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGgcg                                                                                                                                <  2:51679/150‑1 (MQ=255)
                                   aaaaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTgcgc                                                                                                                    <  1:100852/150‑1 (MQ=255)
                                   aaaaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTgcgc                                                                                                                    >  1:2395343/1‑150 (MQ=255)
                                    aaaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGcgcg                                                                                                                   <  2:1222833/150‑1 (MQ=255)
                                    aaaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGcgcg                                                                                                                   <  1:952632/150‑1 (MQ=255)
                                     aaaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                                  >  2:766293/1‑150 (MQ=255)
                                      aaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCggg                                                                                                                 <  1:2096317/150‑1 (MQ=255)
                                      aaGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCggg                                                                                                                 <  1:766293/150‑1 (MQ=255)
                                          gATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtgg                                                                                                             <  2:2282087/150‑1 (MQ=255)
                                           aTACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggt                                                                                                            <  2:2395343/150‑1 (MQ=255)
                                           aTACACCGACGAAGGCGGAATCTGCGCCCGCTCAGCCACGCTGTAGCTATAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggt                                                                                                            >  1:792276/1‑150 (MQ=255)
                                            tACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggtg                                                                                                           <  1:591501/150‑1 (MQ=255)
                                               acCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATGCTGGGCTGACTGTTGTCTCGTGGTTTGCGCGGGgtggtggtg                                                                                                        >  2:998468/1‑150 (MQ=255)
                                               acCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                                  <  1:1563278/140‑1 (MQ=255)
                                               acCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                                  >  2:1563278/1‑140 (MQ=255)
                                               acCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGgtggtggtg                                                                                                        >  2:1411737/1‑150 (MQ=255)
                                                 cgacgaAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCg                                                                                                      <  2:1083112/150‑1 (MQ=255)
                                                  gacgaAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCg                                                                                                      <  1:516477/149‑1 (MQ=255)
                                                  gacgaAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCg                                                                                                      >  2:516477/1‑149 (MQ=255)
                                                  gacgaAGGCGGAATCTGCGCCCGCTCGGCCACGCGGTAGATAGAACGTGCTGCCACCCCTTATTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGACCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCgg                                                                                                     >  1:425939/1‑150 (MQ=255)
                                                        ggCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatc                                                                                               <  1:868170/150‑1 (MQ=255)
                                                        ggCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatc                                                                                               >  1:741419/1‑150 (MQ=255)
                                                         gCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatca                                                                                              >  2:1941509/1‑150 (MQ=255)
                                                           ggAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAc                                                                                            <  2:393655/150‑1 (MQ=255)
                                                             aagcTGCGCCCGCTCGGCCACGCTGTAGATAGAGCCTCCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGACCCTGGCCTGGCTGCTGGCGCGTGGTGTGCGCGGGGCGGTGGTGCGGTTCATCAActg                                                                                          <  2:633384/147‑1 (MQ=255)
                                                              aTCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAActgc                                                                                         <  1:1533419/150‑1 (MQ=255)
                                                                   cgcCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGCCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGACCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGCGGTGGTGGTGCGGTTCATCAACTGCTGAtt                                                                                    <  2:425939/150‑1 (MQ=255)
                                                                       cGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCa                                                                                <  1:956996/150‑1 (MQ=255)
                                                                        gCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGGTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAg                                                                               <  1:924456/150‑1 (MQ=255)
                                                                          tCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatc                                                                                               <  2:821716/132‑1 (MQ=255)
                                                                          tCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcatc                                                                                               >  1:821716/1‑132 (MQ=255)
                                                                              ccACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCGATTTCTGATTGCCAGGATGCt                                                                         >  2:2361293/1‑150 (MQ=255)
                                                                              ccACGCTGGAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCATCTGCTGATTGCCAGGATGCt                                                                         <  2:754605/150‑1 (MQ=255)
                                                                               cACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTg                                                                        <  2:608153/150‑1 (MQ=255)
                                                                                 cGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGc                                                                      <  1:1472030/150‑1 (MQ=255)
                                                                                    tGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGtt                                                                                                   <  1:1211865/118‑1 (MQ=255)
                                                                                    tGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGtt                                                                                                   >  2:1211865/1‑118 (MQ=255)
                                                                                     gTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGAGGCGATcctggc                                                                                                                                             <  2:1125176/75‑1 (MQ=38)
                                                                                     gTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGAGGCGATcctggc                                                                                                                                             >  1:1125176/1‑75 (MQ=38)
                                                                                       agatagaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTa                                                                <  2:1504465/150‑1 (MQ=255)
                                                                                         atagaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcat                                                                                                >  2:736724/1‑116 (MQ=255)
                                                                                         atagaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTtcat                                                                                                <  1:736724/116‑1 (MQ=255)
                                                                                          tagaAAGTGCTGCCACGCTTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGgc                                                             <  2:792276/150‑1 (MQ=255)
                                                                                            gaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAct                                                                                           <  2:904468/118‑1 (MQ=255)
                                                                                            gaAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAct                                                                                           >  1:904468/1‑118 (MQ=255)
                                                                                             aaaGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGcc                                                          <  2:1778904/150‑1 (MQ=255)
                                                                                                gTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCcgac                                                       <  1:2059781/150‑1 (MQ=255)
                                                                                                    tgcCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAAc                                                   >  2:1498485/1‑150 (MQ=255)
                                                                                                        aCCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGt                                               >  2:1562816/1‑150 (MQ=255)
                                                                                                          cccTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCg                                             >  2:1887518/1‑150 (MQ=255)
                                                                                                            cTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggtg                                                                                                           >  2:1197414/1‑86 (MQ=255)
                                                                                                            cTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGggtggtg                                                                                                           <  1:1197414/86‑1 (MQ=255)
                                                                                                            cTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGaa                                           <  2:957091/150‑1 (MQ=255)
                                                                                                               gTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAc                                                                                            <  1:943738/98‑1 (MQ=255)
                                                                                                               gTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAAc                                                                                            >  2:943738/1‑98 (MQ=255)
                                                                                                                ttCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCa                                       >  1:2106263/1‑150 (MQ=255)
                                                                                                                   ggCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGa                                            >  1:1580699/1‑142 (MQ=255)
                                                                                                                   ggCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGa                                            <  2:1580699/142‑1 (MQ=255)
                                                                                                                   ggCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCg                                    >  2:278930/1‑150 (MQ=255)
                                                                                                                         cAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCg                              >  2:2085567/1‑150 (MQ=255)
                                                                                                                         cAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCg                              <  2:741419/150‑1 (MQ=255)
                                                                                                                           gCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCGGTGACTTCgac                            <  2:471893/150‑3 (MQ=255)
                                                                                                                           gCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGtt                            <  2:2106263/150‑1 (MQ=255)
                                                                                                                            cTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGttt                           <  1:1411737/150‑1 (MQ=255)
                                                                                                                             tCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGtttt                          >  1:1182031/1‑150 (MQ=255)
                                                                                                                             tCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGGTCATCAACTGCTGATTGCCATGAGGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGtttt                          >  1:1741818/1‑150 (MQ=255)
                                                                                                                               aCGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGttt                           <  2:2105572/147‑1 (MQ=255)
                                                                                                                               aCGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCCGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGttt                           >  1:2105572/1‑147 (MQ=255)
                                                                                                                                 gcgacgcc‑‑gctAGGATCGCTTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTc                      >  2:600258/6‑150 (MQ=255)
                                                                                                                                     aGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAg                                                 <  2:1235782/119‑1 (MQ=255)
                                                                                                                                     aGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAg                                                 >  1:1235782/1‑119 (MQ=255)
                                                                                                                                     aGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGc                  >  2:2015883/1‑150 (MQ=255)
                                                                                                                                       cc‑‑gccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTc                <  1:998468/150‑1 (MQ=255)
                                                                                                                                           gccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGaa                                                    <  1:1570321/112‑1 (MQ=255)
                                                                                                                                           gccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGaa                                                    >  2:1570321/1‑112 (MQ=255)
                                                                                                                                           gccAGGATCGAGGCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCgg              >  1:183914/1‑150 (MQ=255)
                                                                                                                                            ccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGa             <  1:1941509/150‑1 (MQ=255)
                                                                                                                                                gATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                                  >  1:1053330/1‑45 (MQ=255)
                                                                                                                                                gATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCgg                                                                                                                  <  2:1053330/45‑1 (MQ=255)
                                                                                                                                                gATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGc         >  1:1042964/1‑150 (MQ=255)
                                                                                                                                                      cTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGTGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgcg   >  2:767646/1‑150 (MQ=255)
                                                                                                                                                      cTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGc                        >  1:109345/1‑129 (MQ=255)
                                                                                                                                                      cTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGc                        <  2:109345/129‑1 (MQ=255)
                                                                                                                                                      cTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgcg   >  2:1764718/1‑150 (MQ=255)
                                                                                                                                                       tCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgcga  <  1:1562816/150‑1 (MQ=255)
                                                                                                                                                       tCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgcga  <  1:1887518/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CAGTAAGTGCCTCCAGCAAGGCCGGCACGCTGGCGAAAAAGTGATACACCGACGAAGGCGGAATCTGCGCCCGCTCGGCCACGCTGTAGATAGAAAGTGCTGCCACCCCTTGTTCGGCCAACAGCTCACGCGCAGCC‑‑GCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGA  >  NC_002947/296223‑296521

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